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The  data  base  previously  developed  for  multicomponent 
SIALON  Ceramic  phase  diagrams  has  been  expanded  to  cover 
Ce2°3'  Be0  and  ^2°3  Editions.  Isothermal  sections  in 
the  Mg0-Si3N4-Si02,  Y203~Si02-Si3N4  and  Ce203-Si02“Si3N4 
systems  near  2000°K  have  been  computed  and  compared  with 
limited  experimental  data.  The  trajectory  of  ordering 
temperatures  for  A2/B2  and  B2/DC>3  reactions  has  been 
computed  along  the  Fe3Si-Fe3Al  composition  path  in  the 
BCC  phase  of  the  Fe-Al-Si  system  and  compared  with  experi¬ 
ment.  The  two  phase  (fee  +  bcc)  fields  for  ordered  phases 
in  the  iron-aluminum-nickel, iron-aluminum-manganese,  and 
the  iron-nickel-manganese  system  between  700°C  and  1200°C 
were  calculated.  Construction  of  a  data  base  covering 
fluoride  systems  consisting  of  systems  containing  ZrF4 
which  are  employed  to  synthesize  fluoride  glasses  has  been 
initiated  and  used  to  calculate  the  composition  of  maximum 
liquid  stability  in  the  ZrF4-LaF3~BaF2  and  the  ZrF4-BaF2~NaF 
ternary  systems  where  glass  formation  has  been  observed. 

The  calculations  have  been  extended  to  consider  the  effects 
of  A1F3  additions  on  the  glass  compositions  with  good 
results.  An  analysis  of  the  titanium-carbon-nitrogen 
system  coupling  the  thermochemical  and  phase  diagram  data 
has  been  performed  to  permit  calculation  of  the  ternary 
phase  diagram  and  thermochemical  properties  over  a  range 
of  temperatures.  This  information  should  be  useful  in 
vapor  deposition  of  coatings  based  on  this  system  which 
are  used  in  metal  cutting  and  wear  resistance  applications. 
The  addition  of  Zr02  to  the  data  base  has  been  initiated 
and  applied  to  calculation  of  the  Zr02~Al03  and  Zr02-Y2C>3 
binary  system  as  well  as  isothermal  sections  in  the  Al203~ 
Zr02~Y203  ternary  system  between  1427°C  and  2427°C. 
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INTRODUCTION  AND  SUMMARY 


U.S.  Air  Force  requirements  for  an  expanding  variety  of 
materials  and  energy  sources  with  the  constraints  imposed 
by  ever  increasing  costs  and  budget  restrictions  has  magnified 
the  need  for  evaluating  current  methods  of  materials  processing 
and  performance.  The  most  rational  basis  from  which  to 
perform  such  an  evaluation  is  an  extensive  thermodynamic 
frame  work  for  describing  and  predicting  the  behavior  of 
materials  systems  over  a  broad  range  of  environmental  conditions. 
Development  of  such  a  frame  work  would  permit  interrelation  of 
phase  diagram  and  thermochemical  data  yielding  a  more  complete- 
description  of  a  given  system.  The  utility  of  computer  based 
methods  for  developing  and  utilizing  coupled  phase  diagram 
and  thermochemical  data  on  metallic  systems  has  been  well 
documented,  the  present  work  is  directed  toward  achieving  a 
similar  facility  in  dealing  with  ceranic,  glass  and  compound 
based  systems  used  in  electronic,  metal  cutting  and  high 
temperature  structural  applications. 

Section  II,  pages  5-15,  published  CALPHAD  (1981)  5  163-84) 
extends  the  present  data  base  to  BeO,  Y2C>3  and  Ce2C>3  ceramics. 
Previous  segments  of  this  base  cover  combinations  of  Cr203  • 

MgO,  A1203,  Fe203,  Fe3C>4,  "FeO” ,  Si02,  CaO,  Si2N4  and  AIN. 

Lattice  Stability,  Solution  and  Compound  Phase  Parameters 
were  derived  covering  the  liquid,  spinel,  corundum,  periclase, 
crystobalite,  tridymite,  quartz,  hexagonal  and  beta  prime  phases 
which  appear  in  the  binary  systems  composed  of  pairs  of  these 
compounds.  Compound  phases  formed  from  specific  binary 
combinations  of  these  compounds  (i.e.  Mg0*Cr203)  were  also 
characterized.  This  description  is  based  on  observed  thermo¬ 
chemistry  and  phase  diagrams  for  the  binary  systens  of  interest. 
Selected  ternary  systems  have  been  computed  based  on  the  fore¬ 
going  data  base  for  comparison  with  experimental  sections  in  order 
to  illustrate  the  usefulness  of  the  data  base. 

Section  III,  pages  16-28,  (published  in  High  Temperatures- 
High  Pressures  (1982)  14^  619-31)  illustrates  application  of  the 
data  base  to  the  processing  and  performance  of  SIALON  Ceramics 
by  calculation  of  quasiternary  sections  in  the  Si02-Mg0-Si3N4 , 


SiOj-Ce-p j-Si^N^  and  Si02-Y2^2_  si3N4  ^etween  1900  and  2100K 

which  agree  well  with  experimental  findings. 

In  order  to  illustrate  application  of  these  techinques 
to  alloy  semiconductors,  Section  IV,,  pages  29-58,  deals  with 
development  of  a  data  base  covering  the  binary  systems  composed 
of  Aluminum,  Gallium,  Indium,  Phosphous  ,  Arsenic  and  Antimony 
has  been  constructed  by  analyzing  the  fifteen  combinations  of 
these  elements  in  terms  of  lattice  stability,  solution  phase 
and  compound  parameters.  These  compound  semi-conductors  are 
important  for  a  wide  range  of  electro-optical  applications 
because  then  properties  can  be  tailored  by  controlling  the 
composition  of  the  solid,  which  is  grown  from  the  liquid 
or  deposited  from  the  vapor.  Although  crystal  growers 
have  become  very  adept  at  controlling  the  temperature  (and 
pressure)  of  the  parent  phase  in  order  to  obtain  the  desired 
characteristics  of  alloy  semiconductors,  phase  diagram  data 
is  indi spensi ble  in  order  to  deal  with  increasingly  complex 
systems..  Partial  isothermal  sections  in  the  P-In-As,  As-In-Sb, 
P-Ga-As,  Ga-Sb-In  and  Al-Sb-Ga  systems  were  then  calculated 
using  the  foregoing  data  base  for  comparison  with  experimental 
isothermal  sections  and  quasi-binary  III-V  phase  diagrams. 

It  was  found  that  ternary  liquid  and  III-V  compound  interaction 
parameters  were  required  to  attain  good  agreement  in  some  cases. 
Similar  calculations  were  performed  for  the  Te-Cd,  Hg-Cd,  and 
i'e-Hg  binary  systems  and  the  Cd-Te-Hg  ternary  systems  at  pressures 
up  to  74  atmospheres.  Comparison  of  the  calculated  results  with 
experimental  data  on  tie-line  compositions  between  Cd-Te-Hg 
liquid  and  quasi-binary  CdTe-HgTe  alloys  is  important  in  the 
liquid  phase  epitaxial  growth  of  controlled  band  gap  electro- 
optical  materials.  Pages  29-58  comprise  a  paper  published  in 
CALPHAD  jj  (1981)  217  covering  this  work. 

The  recent  discovery  of  a  new  family  of  non-oxide  glasses 
based  on  mixtures  of  ZrF^  or  HfF^  with  other  metallic 
fluorides  by  M.  Poulain  and  coworkers  offers  great  potential 
in  optical  fiber,  window  and  source/detector  applications. 

Due  to  the  limited  phase  diagram  data  available  for  the 
binary,  ternary  and  multicomponent  fluoride  systems  currently 
employed  to  synthesize  these  glasses  most  of  the  progress 
in  identifying  new  compositions  has  proceeded  along  empirical 
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lines.  In  order  to  remedy  this  situation,  the  development 
of  a  data  base  covering  fluoride  glasses  has  been  undertaken. 
Coupled  thermochemical  analysis  of  the  ZrF4~BaF2,  ZrF^-LaF^, 
ZrF4-PbF2,  ZrF^NaF ,  ZrF4~KF ,  ZrF4~RbF  and  ZrF^CsF,  LaI’3-BaF2 
and  BaF^-NaF  systems  was  completed.  The  results  were  employed 
to  calculate  the  composition  of  maximum  liquid  stability  in 
the  ZrF4-LaF3“BaF2  and  the  ZrF4-BaF2-NaF  ternary  systems  where 
glass  formation  has  been  observed.  Section  V,  pages  59-G9 
(published  in  CALPHAD  (1983)  2  71-82)  illustrates  these  results. 

Requirements  for  engine  components  in  advanced  air  craft 
engine  designs  have  stimulated  interest  in  iron  base  alloys 
exhibiting  ordered  phases  which  are  stable  at  temperatures 
between  900  and  1500K.  Section  VI,  pages  71-84,  presents 
a  calculation  of  ordering  temperatures  in  the  quasi-binary  section 
between  Fe2Al  and  Fe^Si.  Section  VII,  pages  85-103,  examines 
ordering  reactions  in  the  iron-aluminum  system  and  the  effects 
of  ordering  on  the  fcc/bcc  equilibria  in  iron  rich  alloys 
containing  aluminum,  manganese  and  nickel  between  973K  and  1473K. 

The  use  of  titanium  carbide,  titanium  nitride  and  titanium 
carbo  nitride  as  cutting  tool  materials  or  as  coatings  for 
cutting  tools  and  wear  parts  has  expanded  substantially  during 
the  past  decade.  Section  VIII,  pages  104-129,  provides  a 
thermodynamic  evaluation  of  the  Ti-C,  Ti-N  and  Ti-C-N  systems 
covering  the  range  of  processing  and  performance.  Ternary 
isothermal  sections  are  presented  between  1600K  and  2500K  and 
isoactivity  lines  for  the  elemental  components  are  presented 
at  1800K. 

Ceramic  composites  based  on  Zr02  or  in  which  Zr02  is  present 
as  a  second  pahse  have  wide  applications  in  high  temperature 
technology.  Recently  it  has  been  found  that  composites  which 
contain  Zr02  dispersions  can  exhibit  enhanced  fracture  toughness 
due  to  a  mechanism  which  is  based  on  the  titrogonal-monoclinic 
transformation  in  Zr02.  Such  enhancement  has  been  reported 
for  Al202,  S*3N4  and  SiC-A^O^-SiC^  composites.  Accordingly 
Section  IX,  pages  130-136,  describes  the  initial  step  in  adding 
Zr02  to  the  data  base  covering  the  Zr02-Y203  and  Zr02“Al203 
binary  systems  and  the  Al203-Zr02-Y203  ternary  systems. 

Section  X,  pages  137-38,  provides  a  calculation  of  the 
heat  of  fusion  of  eutectic  glasses  in  the  ZrF  .-BaF_ 


and  ZrF^-LaF^  systems  while  Section  XI,  pages  139-160, 
provides  a  calculation  of  the  effect  of  AIF^  additions  on 
glass  compositions  of  ternary  ZrF^-LaF^-BaJ  2  fluorides 
which  agree  with  observations.  Additions  of  AlF^  are  found 
to  facilitate  fabrication  of  large  fluoride  glass  windows. 
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II*  CALCULATION  of  quasi-binary  ceramic  systems 
Kaufman,  F.  Hayes  and  D.  Birnie,  CALPHAD  (1981)  5  163-184 


(This  paper  was  presented  at  CALPHAD  X,  Vienna,  Austria,  July  1981) 

ABSTRACT.  A  data  base  is  being  developed  for  calculation  of  quasi-binary  and  quasi-ternary 
phase  diagrams  of  ceramic  systems  (1-3).  Previous  segments  of  this  base  cover 
combinations  of  Cr^,  MgO,  Al^,  Fe^,  Fe^ ,  "FeO",  Si02>  CaO,  Si^  and  AIN. 

Lattice  Stability,  Solution  and  Compound  Phase  Parameters  were  derived  covering 
the  liquid,  spinel, corundum,  pcriclase,  crystobalite,  tridymitc,  quartz,  hexagonal 
and  beta  prime  phases  which  appear  in  the  binary  systems  composed  of  pairs  of 
these  compounds.  Compound  phases  formed  from  specific  binary  combinations  of 
these  compounds  (i.e.  MgO-Cr.O.)  were  also  characterized.  This  description  is 
based  on  observed  thermochemistry  and  phase  diagrams  for  the  binary  systems  of 
interest.  Selected  ternary  systems  have  been  computed  based  on  the  foregoing  data 
base  for  comparison  with  experimental  sections  in  order  to  illustrate  the  usefulness 
of  the  data  base.  The  present  paper  extends  the  data  base  to  cover  BeO,  Y-,0^  and 

Ce2°3  additions*  Moreover,  ternary  sections  in  the  SiO^-MgO-Si.^ ,  Si02-Y203-Si ^ 
and  Si02*Ce20j-Si were  calculated  between  1900K  and  2100K  for  comparison  with 
experiment. 

I.  Introduction 

The  utility  of  computer  based  methods  for  coupling  phase  diagrams  and  thermochemicai 
data  for  metallic  systems  has  been  well  documented  in  many  papers  published  in  this  journal. 

A  considerable  effort  is  being  applied  toward  developing  an  extensive  base  for  metallic 
systems.  Recently  similar  efforts  have  begun  in  order  to  provide  a  similar  facility  with 
ceramic  systems.  The  expanded  studies  of  S1AL0N  composites,  combining  silicon  and  aluminum 
nitrides  with  oxides  of  silicon,  aluminum,  magnesium,  beryllium,  cerium  and  yttrium  and 
other  metals  has  provided  additional  motivation  for  predicting  multicomponent  phase  diagrams 
of  ceramic  systems. 


2.  Description  of  the  Thermochemicai  System 
ployed  to  Characterize  Solution  and  Compound  Phases 


The  method  utilized  for  describing  solution  and  compound  phases  is  the  sane  as  that 
employed  earlier  (1-3)  comprising  some  symbolic  usage  which  facilitates  data  handling  as 
indicated  below.  The  free  energy,  G  ,  of  a  liquid  (solution)  phase,  L,  in  the  binary  system 
W8*0  4*  *ivw*  by  ER****10"  ( 1)  where  T  is  in  Kelvins,  R-6.314  J/g.  at*K,  while 

x  is  the  atomic  fraction  of  BO(i.e.  BO-ijBeO)  and(l-x)is  the  fraction  of  YO(i.e.  1/SY203). 

The  ■»**  basis  is  thus  one  mole  of  atoms  (i.e.  a  gram-atom) 


GL=>(l-x)*GY0*xbg0»RT(xfnx*(l-x)fn(l-xJ*x(l-x)llYOBO(l-x)»xLBOYO]  J/g. at  (1) 

Were  *GyQ  and  *GgQ  are  the  free  amrgies  of  one  gram  atom  of  pure  1  iquid  YOandBO.  Table  1  defines 


the  lattice  stabilities  of  the  liquid  and  solid  forms  of  YO,  CO  and  the  other  components  of 
current  interest.  These  data  derived  from  earlier  studies  and  compilations  of  thermochemicai 
and  phase  diagram  data,  (1-8),  when  combined  with  the  solution  and  compound  phase  parameters 
shown  in  Tables  2  and  3  permit  calculation  of  the  binary  phase  diagrams  shown  in  Figures  1-12, 
The  solution  parameters  LYOBO  and  LBOYO,  which  describe  the  liquid  Y  O.-BeO  solution  are 
listed  in  Table  2  (LYOBO* LBOYO* -209 2  Joules/g.at) .  Similar  parameters''for  the  liquid  ami 
solid  phases  are  listed  in  Table  2.  The  free  energy  of  the  solid  phases  are  described  in  a 
manner  similar  to  equation  (1).  Thus,  for  example,  the  free  energy  of  the  body  centered 
cubic  (Mn20j)Y  phase  in  the  YO-BO  system  is  given  by 


GY*(l-x)*Gy0*xtg0*RT(xfnx*(l-x)«nf l-x))*x(l-x)  (YY0B0(l-x)*xYB0Y0]  J/g. at.  (2) 
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TABLE  1 


SUMMARY  OF  LATTICE  STABILITY  PARAMETERS 
(All  units  ir.  Joules  per  gram  atom  (mole  of  atoms),  T  in  Kelvins) 

AO  =■  (l/SJAl^Oj,  SO  =  (l/J)Si02,  MO  =  (l/2)MgO,  AN  =  ( 1/2) AIN,  SN  =  (l/7)Si3N4 

BO  =  ( 1/2) BeO,  YO  =  (l/SJY^,  CE  =  (l/S)Ce203 

P  *  Poriclase,  C  =•  Corundum,  X  =  Crystobalite,  T  =  Tridymite,  II  =  a  quartz 

R  -  3  quartz,  B  =  beta  Si.^,  N  =  hexagonal  AIN  and  BeO,  Y  =  body  centered  cubic  (Mn^^Y^Oj 
and  Ce,03  structure 

BOBOLN*  =  (1/2) BeO( Liquid) - (1/2) BeO( hexagonal) 

YOYOLY  *  (l/S)Y203(Liquid)-(l/S)Y203  (body  centered  cubic) 


YOYOLY 

3 

22694  - 

8.368T 

CECELY 

a 

20334 

- 

8.368T 

YOYOLN 

3 

0  - 

14.142T 

CECELN 

a 

0 

14.142T 

YOYOLX 

ss 

0  - 

2.092T 

CECELX 

? 

0 

- 

2.092T 

YOYOLB 

3 

0  - 

12.S10T 

CECELB 

= 

0 

- 

12.S10T 

YOYOLC 

a 

0  - 

10.209T 

CECELC 

3 

0 

10.209T 

BOBOLN 

a 

40376  - 

14.142T 

AOAOLY 

s 

0 

8.-36ST 

BOBOLX 

a 

0  - 

1.674T 

ANANLY 

= 

0 

- 

8.368T 

BOBOLC 

a 

0  - 

10.209T 

SNSNLB 

= 

33949 

- 

12.510T 

BOBO LB 

= 

0  - 

12.S10T 

SNSNLY 

= 

0 

- 

8.368T 

BOBOLY 

a 

0  - 

8 . 368T 

SOSOLY 

= 

0 

- 

2.092T 

SOSOLN 

E 

0 

- 

2.092T 

SOSOU 

a 

3347 

- 

1.674T 

SOSOLR 

3 

5042 

- 

3.096T 

*  These  differences  specify  the  free  energy  of  one  phase  (i.e.  liquid)  mi  us  the  free  energy 
of  the  second  phase  (i.e.  hexagonal)  for  a  given  compound. 


TABLE  2 

SUMMARY  OF  SOLUTION  PHASE  PARAMETERS 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms),  T  in  Kelvins) 

L  =  Liquid,  N  *  hexagonal,  X  =  Crystobalite,  T  «  Tridymite,  Y  •  body  centered  cubic 


h  =  beta  Si3N4,  AO  =  (l/SJAl^j,  SO  =  (l/3)Si02,  AN  «  (1/2) AIN,  SN  -  (l/7)Si3N4,  BO  =  (l/2)BeO 


YO  =  (1/S)Y203, 

CE  = 

(l/5)Ce203 

LSOBO  * 

LBOSO  * 

41840 

LBOSN 

- 

LSNBO 

20920 

0<xcr<0.60* 

NSOBO  => 

NBOSO  = 

83680 

BBOSN 

= 

BSNBO 

— 

83680 

-  so— 

XSOBO  =• 

XBOSO  =■ 

83680 

NBOSN 

= 

NSNBO 

2 

83680 

LYOSO  ■-  LSOYO  =  -92048+36. 82T 

TSOBO  * 

TBOSO  = 

83680 

0.6<xso<1.0 

LYOSO  =  -2 16, 000+66. 94T 

LYOBO  * 

LBOYO  = 

-2092 

LYOAO 

3 

LAOYO 

3 

12552 

LSOYO  »  -9414+16. 74T 

iYuBO  = 

YBOYO  =■ 

-2092 

YYOAO 

3 

YAOYO 

3 

83680 

YYOSO  =  YSOYO  «  83680 

NYOBO  * 

NBOYO  = 

83680 

CYOAO 

= 

CAOYA 

3 

83680 

XYOSO  =  XSOYO  =  83680 

TYOSO  *  TSOYO  *  83680 

LAOBO  * 

LBOAO  « 

8368 

LYOAN 

E 

LANYO 

B 

20920 

CAOBO  * 

CBOAO  * 

83680 

YYOAN 

8 

YANYO 

3 

83680 

0<xcnl0-60 

NAOBO  * 

NBOAO  * 

83680 

NYOAN 

8 

NANYO 

8 

83680 

—  SO— 

LCESO  »  LSOCE  •  -92048  ♦  36.82T 

1.  BOAS’  = 

LANBO  » 

20920 

LSNSO 

a 

LSOSN 

X 

29288 

0.6<xs£><1.0 

NBOAN  * 

NAN BO  = 

83680 

XSNSO 

= 

XSOSN 

3 

125520 

LCESO  ♦  -216, 000+66. 94T 

TSNSO 

= 

TSOSN 

■ 

125520 

LSOCE  =  -9414+16. 74T 

LCESN  = 

LSNCE  = 

4180 

YCESO  =  YSOCE  =  83680 

BCESN  = 

BSNCE  = 

62760 

LYOSN 

3 

LSNYO 

a 

4180 

XCESO  =  XSOCE  =  83680 

YCESN  = 

YSNCE  * 

62760 

BYOSN 

= 

BSNYO 

■ 

83680 

TCESO  =  TSOCE  =>  83680 

YYOSN 

YSNYO 

3 

83680 

*x  Q  *  atomic  fraction  SO  in  YO- 

or  CE-SO 


TABLE  3 


S IBMARY  OF  COMPOUND  PARAMETERS  FOR  BINARY  SYSTEMS 


(All  units  in 

Joules  per  gram  atom 

(mole  of 

atoms)  T  in 

Kelvins) 

Compound 

Name 

Stoichiometry 

Stability 

Base 

Compound  Parameter 

(Joules/g. at . ) 

(l/7)(Si02-2BeO) 

G 

SO0.429BO0.S71 

stable 

N 

(>2132-7.  106T 

(l/17)(5Al2OjBeO) 

q 

A00.882B00.118 

stable 

C 

75312 

(1/7) (AljOj-BeO) 

i 

AO0.714BO0.286 

stable 

C 

7  (>567 

(l/ll)(Al203-3Be0) 

j 

A00.45SB0O.S4S 

stable 

N 

01923*9. 91  OT 

(1/1S)(2Y203-A1203) 

U 

YO0.667AO0.333 

stable 

Y 

70291 

(i/io)(Y2o3-ai2o3) 

V 

Y°0. S00A°0.  500 

stable 

c: 

48534*20.08 5T 

(1/40)(3Y203-SA1203) 

WA 

Y00. 375AO0. 625 

stable 

c 

82425 

(I/5)(Si2N20) 

R 

SN0.7SO0.3 

stable 

B 

62760-1 .25ST 

(1/12) (Y203-Si3N4) 

S 

Y00.417SN0.S83 

stable 

B 

44769*14. 644T 

(1/8)(Y203-Si02) 

M 

Y00.625S00.37S 

stable 

Y 

-15900*29.2881 

(l/19)(2Y203-3Si02) 

K 

YOO.S27SO0.473 

stable 

Y 

-30334*36.8191 

(l/ll)(Y203-2Si02) 

0 

Y00.45SS0O.S45 

stable 

Y 

-2218*23. 849T 

(l/12)(Ce203-Si3N4) 

S 

CE0.417SN0.S83 

stable 

B 

44769*14. b4  tl 

(l/8)(Ce203-Si02) 

M 

CEO.625S00.37S 

stable 

Y 

418*24. 058T 

(l/62)(7Ce203-9Si02) 

Z 

CEO.S6SS00.43S 

stable 

Y 

6276*24. 434T 

(l/ll)(Ce203-2Si02) 

0 

CE0.4SSS0O.S45 

Stable 

Y 

7531*25. 230T 

0.01(1 4Si jN4 ■ BeO) 

BP 

SN0.98BO0.02 

stable 

B 

83680*4. 184T 

(1/150) (21Si3N4-Si02) 

BP 

SN0.98S00.02 

stable 

B 

71128*4. 184T 

0.0040(35Si3N4-Y203) 

BP 

SN0.98Y°0.02 

stable 

B 

83680*4. 184T 

0.0040(3SSi3N4-Ce203) 

BP 

SN0.98CE0.02 

stable 

B 

83680*4. 180T 

solution  parameters 

for  the 

Y  phase  (YYOBO»YBOYO= 

-2092  J/g. 

at)  are  listed  in  Table  2. The  free 

. At  M  V 

energy  difference  between  the  L  and  Y  forms  of  YO  and  BO,  i.e.  YOYOLY*‘,C^p-”G^Q  and 

"Go-- *GY  =BOBOLY  are  listed  in  Table  1.  The  free  energy  of  a  compound  phase  such  as 
BO  BU 

C*(  1/7) (SiO2’2Be0)=SO0  <29BO0  571  *s  defined  on  tbe  basis  of  Table  3  in  terms  of  the  base 
phase,  the  compound  parameter,  C,  and  the  stoichiometry  by  equation  (3) 

GG«0.429g!!  *0.S71G^„*(0. 429) (0.571) [0.429LSOBO*O, 571LBOSO  -C]  J/g.at.  (3) 

jU  BU 

Reference  to  Table  2  shows  that  LSOBO=LBOSO=41840  J/g.at  whilo  Table  3  shows  that  C=621 32-7 . 406T 
Thus 

GC=0.429GgO*0.S71GgO-4971*1.813T  J/g.at.  (4) 

At  298K  reference  to  Figure  1  shows  that  the  stable  form  of  quartz  is  the  beta 
form  (designated  R) .  Table  1  shows  that 

Gso‘cso=SOSOLR'SOSOLNl5042'10°4T  J/s-at-  Cr’) 

Substitution  of  equation  (5)  into  eq  it  ion  (4)  yields  an  expression  for  the  free  energy  of 
formation  of  (1/7) (SiO^'BeO)  at  298K  from  the  stable  forms  of  it  s  compound  components  as, 

AGf[298K]»GG-0.429GgO-0.S71CgO*-2808»1.382T  J/g.at. 
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(6) 


TABLF.  4 


CALCULATED  FREE  ENERGY  OF  FORMATION  OF  COMPOUND  PHASES 
(All  Units  in  Joules  per  gran  atom  (mole  of  atoms),  T  in  Kelvins) 


Compound 

Name 

Stoichiometry 

Eree  Energy  of  Formation,  £G-[298K] 
from  Component  ComDOunds  f 

(l/7)(SiO,-2BeO) 

G 

SO0.429BO0.S-l 

-2808  ♦  1.382T/-2809*0.777T(8) 

( 1 / 1 7 ) ( 3A 1 203 • BeO) 

Q 

AO0.882BO0. 118 

-2203  -  0.464T 

(1/7) (Al^Oj- BeO) 

l 

A00. 714B00. 286 

-2381  -  1.12ST 

(l/llJlAl^Oj  • 3BeO) 

J 

A00.4SSB00.54S 

-2524  -  1.049T 

(1/15)(2Y203-A1203) 

u 

Y00.667A00.333 

-4950  -  0.61ST 

(1/io)(y2o3-ai2o3) 

V 

YO0.500AO0.S00 

-2351  -  4.100T 

U/40)(3Y203-SA1203) 

NA 

Y00.37SA00.625 

-7866  ♦  0.690T 

(l/S)(Si2N20) 

R 

SN0.700S00.300 

-5514  -  0.038T/-S781  ■*  0.40T(S) 

(l/12)(Y203-Si3N4) 

S 

YO0.417SN0.583 

-404  -  1 .833T 

(l/8)(Y203‘Si02) 

M 

YO0.62SSO0.375 

-15958  ♦  1.389T 

(l/19)(2Y203-3Si02) 

K 

YO0.S27SO0.473 

-13000  -  0.475T 

(l/ll)(Y203-2Si02) 

0 

YO0.4SSSO0.S4S 

-19S27  -  2.669T 

(l/12)(Ce203-Si3N4) 

S 

CE0.417SN0.S83 

-1389  -  1.828T 

(l/8j(Ce;03-Si02) 

M 

CE0.62SSO0.37S 

-19787  ♦  2.615T 

(l/62)(7Ce203-9Si02) 

2 

CEO.56SSO0.435 

-21974  ♦  2.606T 

(l/ll)(Ce203-2Si02) 

0 

CE0.4SSSO0.54S 

-21945  ♦  2.326T 

0.01(14Si3N4-BeO) 

BP 

SN0.980BO0.020 

-423  -  0.113T 

(l/lS0)(2lSi3N4-Si02) 

BP 

SN0.980S00.020 

-720  -  0.102T 

0. 004 ( 55Si ,N  • Y  0  ) 

3  4  2  3 

BP 

SN0.980Y00.020 

-1104 

0.004(35SiJN4-Ce10J 

BP 

SN0.980CE0.020 

-1158 

Reference  to  Table  4  shows  that  this  result  is  in  keeping  with  the  experimental (8) 
thermochemical  data  on  the  free  energy  of  formation  of  2Be0-Si02.  Table  4  also  displays  the 

calculated  free  energy  of  formation  for  the  remaining  compounds  listed  in  Table  3  and  shown  in 
figures  1-12.  It  should  be  noted  that  this  compound  parameter  for  SNq  ^SO^  3  currently  shown 

as  C  =C27oO- 1 . 25ST  J/g.at.  has  been  revised  from  the  previously  stated  value  C=115060-2S. 10T 
J/g..-!t.(j)fhe  current  value  is  based  on  the  assessment  of  Doemer  ct  al  (5)  and  leads  to  slight 
revisions  in  the  previously  calculated  SN-SO  and  SO-AO-SN  phase  diagrams  presented  earlier  (3). 
lhe  revised  versions,  which  differ  slightly  from  the  earlier  results  are  shown  in  Figures  8 
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w/o  BcO  t 


w/o  BeO-^- 


V2°3  3  7  11  16  22  29  39  53  71 


Figure  2.  Calculated  (l/5)Y2Oj-(l/2)BeO  Phase  Diagram. 
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L(Liquid) 


vs 


w/o  Si_N. - 

3  4 


Figure  11.  Calculated  1/5  Ce^O^  -  l/7Si^N^  Phase  Diagram 


w/o  SiO^— 


III.  CALCULATION  OF  QUASI -TERNARY  CERAMIC  SYSTEMS 


L. 


Kaufman,  F.  Hayes,  and 
High  Pressures 


D.  Birnie,  High  Temperatures- 
(1982)  14  619-31 


Abstract.  Phase  diagrams  lor  ihe  quasibinaiy  systems  Si3N4-Si02,  Y203-Si3N4.  Y203-Si()2, 
(e2()3-Si3N4,  and  Ce203-Si02  are  recalculated  t<>  illustrate  a  method  hv  which  these  and  mote 
complex  equilibria  can  be  elucidated  from  theimodynamic  data.  New  results  lor  the  quasiternary 
systems  Si02-Mg0-Si3N4,  Si02-Y203-Si3N4.  and  Si02-Ce203-Si3N4  aie  derived  from  these  phase 
diagrams  and  presented  as  component  pairwise  equilibria  and  complete  isothermal  sections  of  the 
appropriate  ternary  phase  diagrams.  These  sections  are  compared  with  experimental  data  from  the 
literature,  and  good  agreement  was  obtained  in  the  cases  of  Si02-Mg0-Si3N4  and  Si02-Ce203-Si3N4 

1  Introduction 

The  utility  of  computer  based  methods  for  coupling  phase  diagrams  and  tliermn- 
chcmical  data  for  metallic  systems  lias  been  well  documented  in  many  papers  published 
in  the  journal  CALPHAD.  A  considerable  effort  is  being  applied  toward  developing 
an  extensive  data  base  lor  metallic  systems.  Recently  similar  efforts  have  begun  in 
order  to  provide  a  similar  facility  with  ceramic  systems.  The  expanded  studies  of 
SIAL  )N  composites,  in  which  silicon  and  aluminum  nitrides  arc  combined  with 
oxides  of  silicon,  aluminum,  magnesium,  beryllium,  cerium,  yttrium,  and  oilier  metals 
has  provided  additional  motivation  for  predicting  multicomponent  phase  diagrams  of 
ce  anno  systems. 

Previous  work  in  this  laboratory  on  the  development  of  a  data  base  for  the 
akulation  of  qujMbmary  and  quasitemury  phase  diagrams  of  ceramic  systems 
(Kaufman  and  Nesor  I ‘>78:  Kaufman  ll>7‘>a,  1 979b)  has  covered  combinations  of 
<  r,U3.  MgO.  Al,0„  l  ejOj,  f  e,04.  Te(>\  SiU,.  CaO,  Si3N4.  and  AIN.  Lattice 
viability,  solution,  and  compound  phase  parameters  were  derived  covering  the  liquid, 
spinel,  corundum,  periclase,  crystohalite.  tridymite,  quartz,  hexagonal,  and  P'  phases 
which  appear  in  the  binary  systems  composed  of  pairs  of  these  compounds. 
Compound  phases  formed  from  specific  binary  combinations  of  these  compounds 
(eg  MgO  CrjO.l  were  also  characterized  from  experimental  observations  of  the  binary 
'Ictus  of  interest.  Selected  ternary  systems  have  been  computed  from  this  data  base 
tor  comparison  with  experimental  sections  in  order  to  illustrate  the  usefulness  of  the 
data  base  In  this  paper  the  data  base  is  extended  to  cover  Y203  and  C'e20, ;  ternary 
sections  m  the  SiO,- MgO-Si ,N4.  Si()2- Y2()3-Si3N4.  and  Si02-Ce203-Si3N4  systems 
were  also  calculated  between  1‘XIO  K  and  2100  K  for  comparison  with  experiment. 

2  Thermochcmio.il  sy  stem  for  characterization  of  phases 

The  method  utilized  for  describing  solution  and  compound  phases  is  t he  same  as  that 
employed  earlier  (Kaufman  and  Nesor  1978).  and  has  been  devised  to  facilitate  the 
manipulation  of  data.  The  free  energy,  Gl,  of  a  liquid  (solution)  phase,  L,  in  the 
binary  system  Si,N4-.Si02  is  given  by  equation  (I )  where  /'  is  the  absolute 


temperature,  R  is  the  gas  constant,  x  is  the  atom  fraction  of  SO  (ic  iSiO,),  and 
(1  —  jc)  is  the  atom  fraction  of  SN  (ic  4SijN4).  The  amount  of  substance  basis  is 
thus  one  mole  of  atoms  (ie  a  gram  atom): 

=  (l-x)0G&M+x°C|o  +  K7’I*ln*+(l  ~x)ln(  I  -*)]  +  *( I  -,v)|(l  -,r)CU 

+  xGsosn  ]  •  <  1 ) 


Table  1.  Expressions  for  lattice  stability  parameters  for  various  phases  of  Si02,  Si}N4,  YjOj,  and 
CejO*.  Sources  of  data:  Levin  et  al  (1964);  Kaufman  (1979b). 


°gVo-°Gyo  =  22694-8-3687- 
°Gyo  “  °G*o  =  0-2  0927- 
°Gyo - °Gyo  =  0-  12-5107- 

°GtK-°G^E  =  20334- 8 -638 r 
0G£e-°g£e  =  0-2-0927- 
°6ce-°®ce  =  0-  12-510T 


°Gsn - °Gsn  =  33949-  1 2  5 1 0 r 
°Gsn  “  °Gsn  =  0-8-3687- 

°Gso-°Gso  =  3347-  I  -674r 
°Gko  ~°G'so=  5042-  3-  09<>  T 
°Gso-°G?o-  0-2-09:7- 


Phases:  L,  liquid;  Y,  body-centered  cubic;  X.  crystobalile;  B,/3-Si,N4;  R,(?-quart/.  T,  Iridymiti-. 
Atom  bases:  YO.iYjOj,  SN,4Si,N4;  TE.ICcjO,;  SO.  iSiO,. 

G  =  G/l  (gatom)"1;  T=  T/K. 


"'Stfl./wl'-’ 


Si,N4  10  20  30  40  50  60  70  80  ‘8)  SiO, 


(b)  S  SN 

Figure  1.  Calculated  phase  diagrams  for  binary  systems  (see  tables  1  -3  for  explanation  of  symbols). 
(a)Si,N4-SiO,;  (b)  YjO,-Si,N4;  (c)  Y20,-Si02;  (d)  Ce,0,-Si,N4 .  (c)  CcjO,-SiOj. 


Here  and  °G so  are  the  free  energies  of  one  gram  atom  of  pure  liquid  SN  and 
SO.  Table  1  gives  the  lattice  stabilities  of  the  liquid  and  solid  forms  of  SN,  SO.  and 
the  other  components  of  current  interest.  These  data,  derived  from  earlier  studies 


Figure  I  (continued). 


and  compilations  of  thermochemical  data  and  phase  diagrams,  can  be  used  to 
calculate  the  binary  phase  diagrams  shown  in  figure  1  when  combined  with  the 
parameters  for  solution  and  compound  phases  shown  in  tables  2  and  3.  The  solution 
parameters  G^nso  and  Gsosn.  which  describe  the  liquid  SijN4-Si02  solution,  are 
listed  in  table  2,  together  with  similar  parameters  for  the  liquid  and  solid  phases. 

The  free  energy  of  the  solid  solution  phases  are  described  in  a  manner  similar  to 
equation  (I).  Thus,  for  example,  the  free  energy  of  the  body-centered  cubic 
(MnjOj)  Y  phase  in  the  YaOa-SiOa  system  is  given  by 

GY  =  ( I  -x)0<7yo  +  -*°Gso  +  RT\x  ln.v+  ( I  —  .v > In ( I  -  v)|+  v(  I  -.01(1  -aK-yos,, 


+  *Gsn 


The  free  energy  differences  between  the  L  and  Y  forms  of  YOand  SO,  ic  °Gy0-  Y/y, 
and  °Csb- °^»so.  are  listed  in  table  I.  The  free  energy  of  a  compound  phase  such  is 
R  |ie  i\j(SiaN4-SiOa)  or  SN0  7SN0.j1  is  defined  from  table  3  in  terms  of  the  base 

Table  2.  Summary  of  solution  parameters  for  quasibmary  systems.  Sources  of  data:  levin  el  al 

(1964):  Kaufman  (1979b). 


Gcesn  =  Gsncf.  =  4180  gU,s<)  =  G'sovo  =  -92048  (  36  827"  (0  <-  A'1  <  0  6) 

Gcesn  =  GgNCn  =  62760  G'yoso  =  -216000  + 66-94T  I  ,  ,  .  ... 

Gcfsn  =  G?Nce  =  62760  G'sovo  =  -0414+  16-747"  l(0'6  0> 

Gsnso  =  Gsosn  =  29288  G?oso  =  G^  OYo  =  83680 

Gsnso  =  Gsosn  =  125  520  G$OSo  =  G&>yo  *  83680 

Gsnso  =  g£osn  =  125520  g(OSo  =  gJoyo  =  83680 

Gyosn  =  GkNYo  =  4180  Gceso  =  Gsoce  =  -l»2048+  36-827"  (0  <  t  <  0  6) 

Gyosn  =  G|Nyo  =  83680  G£fso  -  -216000+66-947"  I  <  _ 

Gyosn  -  Gsnyo  =  83680  Gsoce  =  -9414+  16-747"  I  ^ 

Get  so  '  Gsoci  ~  83680 
Gceso  =  Gsoce  =  83680 
Gceso  =  Gsoce  =  83680 

*  x  is  the  atom  fraction  of  SO  present  (see  equation  (2)). 

For  explanation  of  other  symbols,  see  table  I . 


-  Gsncf.  = 


Gsnso  - 
Gsnso  = 


-u 

u  sncf 
Gsncf 
Gsosn 


;* <  i  0) 


(0-6  <  v  <  10) 


Table  3.  Summary  of  compound  parameters  for  binary  systems.  (All  compounds  listed  are  stable 
at  298  K.)  Sources  of  data:  Jack(!977);  Wills  (1974,  1975);  Morgan  (1976):  Levin  et  a)  (1964). 


Gyosn 

Gyosn 

n  Y 

Gyosn 


Gsosn  = 

Gsnyo 

Gsnyo 

Gsnyo 


29288 
125  520 
125520 

=  4180 
=  83680 
:  83680 


Compound 


Name  Stoichiometry 


C/S  (g  atom) 


i(SijNjO) 

R 

SN0.7SO0. 3 

B 

62760-  1  2557" 

A(YjO,Si,N,) 

S 

YOo.4nSN0.  jij 

B 

44769+  14-6447" 

*(Y,0,SiOa) 

M 

YO0.sj5SO0.j75 

Y 

-15900+  29 -2887" 

i(2Y,0,-3Si0,) 

K 

YOo-sa7SOo.47j 

Y 

-30334+36  8197" 

■ir(YiOj-2SiOj) 

0 

YO0.455SOo.M5 

Y 

-2218+23-8497" 

lll(CcjO|-SijN4) 

S 

CE0.4|7SNo.si] 

B 

44769+  14-6447" 

*(Ce,0,-Si0a) 

M 

CE0.8j5SO0.j75 

Y 

418+  24-0587" 

&(7CeaO,-9SiO,) 

Z 

CEo.5e5SOo.8j5 

Y 

6276+  24-4347" 

A(Ce,0,-2SiOa) 

0 

CE0.855SO0.545 

Y 

7531  +25-2307" 

Tfc(21SiJN4SiOa) 

BP 

SNo.9|S00.oi 

B 

71  128+4-1847" 

0  0040(35SijN4  Y,0,) 

BP 

SN0  9|YO0  01 

B 

83680+  4- 1847" 

0-0040(35SijN4-Ce,0,) 

BP 

SNo.$gCE<).<n 

B 

83680+  4  1807" 

*  See  table  1  for  explanation  of  these  and  other  symbols. 

phase,  the  compound  parameter,  C,  and  the  stoichiometry,  by  the  equation 

CR  =0-7oCgN+0-3oG§o  +  (0-7)(0-3)[0-7C^NSO+0  SCljosN-CI.  (3) 

Reference  to  table  2  shows  that  Gsnso  =  G^OSn  =  29  288  i  (g  atom)'1,  while  table  3 
shows  that  C  =  62760-  I -25ST  J  (gatomT1,  where  T  —  T/K.  Thus,  we  obtain 

GR  =  0-7°GfN+0-3oGfo-7029+0-26367\  (4) 

where  G  =  (i/S  (g  atom)"1.  At  298  K  reference  to  figure  1  shows  that  the  stable 
form  of  quartz  is  the  0  form  (designated  R).  Table  I  shows  that 

°Gs0~0Gso  =  5042-  1-0067".  (5) 

Substitution  of  equation  (5)  into  equation  (4)  yields  an  expression  for  the  free  energy 
of  formation  of  i^(Si3N4SiOj)  at  298  K  from  the  stable  forms  of  its  compound 
components  as 

AG,.„g  =  GR-0-7°Gj?N-0-3°G£o  = -5516- 0  038T.  (6) 

Reference  to  table  4  shows  that  the  result  given  in  equation  (6)  is  in  keeping  with 
the  experimental  thermochemical  data  on  the  free  energy  of  formation  of  Si3N4-Si03 
given  by  Doerner  et  al  ( 1 979). 

Table  4  also  displays  the  calculated  free  energy  of  formation  for  the  remaining 
compounds  listed  in  table  3  and  shown  in  figure  1.  It  should  be  noted  that  the 
above  compound  parameter  for  SN0  7SO0  3  has  been  revised  from  the  previously 
stated  value  C  =  115060-25  I0T  J  (g  atom)"1  (Kaufman  1979b);  the  current  value 
is  based  on  the  assessment  of  Doerner  et  al  (1979). 

Table  4.  Calculated  free  energies  of  formation  of  compound  phases  from  component  atom  bases*. 


Compound 

Name 

symbol 

Stoichiometry 

(g  atom) 

\(Si2NjO) 

R 

SNo.7floSOo-.1M 

-55l4-0  038Tb 

AtYjOj  Si,N4) 

S 

YOO.417SN0.5g5 

-404-  1-8337- 

itYjO,  StO,) 

M 

Y00-«:sSOo  375 

-15958+  1-389T 

A(2Y303-3Si0}) 

K 

YOo  jj7SOo-473 

-13000-0-4757- 

AtYjOj  2SiOj) 

O 

YOo  45sSOo  s45 

-19527-  2-669 T 

T^tCcjOj  SijN«) 

S 

CI-O.417SN0  5S3 

-1389-  I  -8287- 

1(  Ce20jSi02 ) 

M 

cto  sjsSOo.  J75 

-19787+  2-6157- 

A(7Ce,Oj  9S.O,) 

Z 

CEo  56sSOo  435 

-21974+  2-6067" 

A(Ce,0,-2SiO,) 

O 

t  Co-4SsSOo  s45 

-21  445+  2-3267" 

2 1  SijN4  SiOj  | 

BP 

SNo  osoSOo-ojo 

-720-0  1207" 

0  004(35SiJN4  YJ0J) 

BP 

SNo- 9*0 YOo  0:0 

-1104 

0  004(35SijN4  Ce3O,) 

BP 

SNophoCEo-ojo 

-1158 

a  See  table  1  for  explanation  of  symbols. 
h  -5781  +0-40T< Doerner  et  al  1979). 


3  Calculation  of  quasiternary  phase  diagrams 

The  free  energies  of  ternary  solution  phases  were  calculated  from  the  binary  solution 
phases  on  the  basis  of  the  Kohler  equation  as  in  previous  papers  (Kaufman  and  Nesor 
1978;  Kaufman  1979a,  1979b).  For  example,  the  free  energy  of  the  liquid  phase  in 
the  SO-MO-SN  system  (MO  =  iMgO),  where  x  is  the  atom  fraction  of  MO  and  y 


the  atom  fraction  of  SN,  is  given  by 

GL  =  (I  -x-y)°G^0  +  x°G +  R T[(  I  - x  -y ) ln(  I  -  x -  »■) 

+  x\nx+y  lnj’j  +  ( I  -  x  -  y)x(  I  ->')'*((  I  -x-  y)G}-OMO  +  vG'^oso) 

+  xy(x  +y T ' [xG ^OSN  +  vG|NMO|  +  (1  -  x - y );•(  I  - x Y 1 1 (I  - x -  v )f,^OSN 

+  >'^SNSol-  (7) 

The  solution  parameters  required  to  specify  Gl  were  given  previously  (Kaufman 
1979a,  1979b)  as  (7sosn  =  ^snso  =  f'MosN  =  ^snmo  =  29  288  J  (g  atom)  1  and 
Gsomo  =  Gmoso  =  - 1 06 274  +  42 •  0 1 T  for  0 -4  <  v  <  I  0  and  G !,OMO  =  -  24  26 7  + 

24-27 T,  Gmoso  =  “220283+  68  627"  for  0  0  <  x  <  0-40.  These  parameters  permit 
explicit  definition  of  GL  in  composition  ranges  where  (x/(  1  —  v  —  >•  1 1  >  0-667.  and 
ranges  where  |x/(l  —  x  — .+ ) |  <  0-667.  In  the  latter  range,  the  miscibility  gap  present 
in  the  SO-MO  binary  system  propagates  into  the  ternary  SO-MO-SN  as  shown  in 
figure  2a.  Here,  the  ternary  miscibility  gap  shows  the  tie  lines  connecting  the 
coexisting  liquid  compositions  L,  and  Lj. 

The  remaining  solution  phases  in  this  system,  ic  P,  X,  B,  etc,  are  defined  in  similar 
fashion  to  that  illustrated  by  equation  (7).  Tire  resulting  equilibria  between  the  L 
and  P  solutions  are  depicted  by  the  tie  lines  traversing  the  two-phase  L+P  field 
shown  in  figure  2b.  The  compound  phases  F,  R,  and  BP  which  appear  in  the  SO- 
MO,  SO-SN,  and  MO-SN  binary  systems  are  defined  along  the  lines  previously 


SN  MO 

(c)  (d) 

Figure  2.  Calculated  pairwise  equilibria  in  the  SiOj-MgO-Si,N4  system  at  2100  K  (see  tables  I  -5 
for  explanation  of  symbols)  (a)F-L;  (b)  L-P  (P.  periclase);  (c)  L-R;  <d)  BP-L. 
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established  (Kaufman  and  N'esor  1978;  Kaufman  1979a,  1979b)  and  are  specified  in 
table  3.  Thus,  for  example,  the  free  energy  of  the  SN0  ?SO0  J  phase  (R)  in  the  SO- 
MO-SN  system  is  defined  as: 

Gr  =  z0G§n+.v0G§o+»'°G&o  +  (i-j^ACa  +  ^AG,, 


+  RT\y  Iny  +  z  Inz  -  (1  —  x)ln(  I  -jc)]  +  AG'F  ,  (8 

where  x,  y,  and  z  are  the  atom  fractions  of  SO,  MO,  and  SN,  respectively.  Since  the 
compound  phase  runs  from  SN9  ,0SO0  J0  to  M00.7oSOq.jo,.v  =  0-3,  and  z  =  0  7-y. 
Also, 

AGa  =  <0-3)(0-7)(0  7GW>  +  0  3G|;osn-0  -  -7029  +  0  2(^T .  (9) 


and 

AGb  =  (0  3)(0-7)(0-7Gfc,OSo+0'3GfeOMo-O  =  -22318+  8-822 T, 

from  tables  3  and  5  and  the  previous  values  of  G^NSO,  Gl-osn-  G'moso  >  and  ^somo 
given  in  the  above  text.  The  excess  free  energy  of  mixing,  AGF,  for  the  compound 
is  given  by 


However,  the  parameter  CAU  is  taken  as  zero  in  all  the  cases  treated  here;  thus 
equation  (8)  can  be  written  explicitly  as  follows: 

GR  =  (0-70-.v)°GfN  +  0  30Gf0  +  y°G5o-  (l  ~  q^)(7°29  -  0-263D  (10 

^(22318-8-822n  +  8-3147Tvlnr-+(0-7-.v)ln(0-7-  +  )-0-7ht0  71 


Table  5.  Summary  of  counterphase  stoichiometry  and  parameters  employed  in  ternary 
calculations.  The  counterphase  parameter  is  taken  to  be  zero  in  all  cases3. 


System 

Stable  phase 
(symbol) 

Base 

Counterphase 

Base 

SO-MO-SN 

SO0  429MO0  S7l0‘  ) 

P 

SN0.439MO0.571 

P 

SOo- 90oS  N0  ■  XX)  ( R ) 

B 

SO0.300MO0.700 

B 

SN0  9wS00.<no(BP) 

B 

SO0.500MO0.500 

B 

SN  0 . 980  M  Oo  •  <no(  B  P) 

B 

MOq.  soqSOo  •  500 

B 

SO-YO-SN 

TOoojsSOo-  j7j(M ) 

Y 

YO0.6j5SN0.375 

Y 

YOfl.  5J7SO|).«7j(K  ) 

Y 

YO0.5j7SN0.473 

Y 

YOo«5jSOoj4s(0) 

Y 

Y00  455SN0.545 

Y 

SO0.,0oSN0  70o(R) 

B 

SOo- jooYOo-too 

B 

SN0.,„SOoo:o(BP) 

B 

SOosnoYOosoo 

B 

SN0  9*oY00  wo(BP) 

B 

YO0  500SO0  500 

B 

SO-CE-SN 

CEo415SOo37j(M) 

Y 

CE0.6J5SN0.  375 

Y 

CEo.S4!SOo.4,5(Z) 

Y 

CE0.545SN0.435 

Y 

CE0455SO0. 545(0) 

Y 

CE0.455SN0.545 

Y 

SOo  30oSN0.70o(R) 

B 

SO0300CE0700 

B 

SNo98oSOo.(no(BP) 

B 

SO0500CE0.500 

B 

SNo9ioCEo.ojo(BP) 

B 

CE0.500SO0.500 

B 

3  See  table  I  for  explanation  of  symbols. 
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Since  the  lattice  stabilities  of  SN,  SO.  and  MO  (ic  °(/'sn  -°(»sn-  “Cso  ~  ('so •  and 
('mo~°('mo)  have  already  been  specified  the  R-L  equilibria  (figure  2c)  can  he 
computed,  figures  2a-2d  illustrate  file  steps  which  are  taken  in  synthesizing  the 
ternary  system  SO-MO-SN.  (ie  ^SiOj-iMgO-4Si3N4)  from  the  component  binary 
systems.  Initially  each  of  the  equilibria  between  pairs  of  solution  and/or  compound 
phases  are  calculated  individually.  In  the  present  case,  these  pairs  consist  first  of  the 
liquid  miscibility  gap  (L,-L2  in  figure  2a),  followed  by  L-P  in  figure  2b  Nest.  I  he 
L-F  and  L-R  equilibria  which  are  shown  in  figures  2a  and  2c  are  calculated.  Super- 
imposition  of  the  Lj-Lj,  L-F.  L-P,  and  L-R  equilibria  shows  tlut  there  arc  no 
interactions  between  any  of  the  phase  pairs  except  in  the  case  of  I  -R  and  L,-l  2. 

The  latter,  which  is  shown  in  figure  2c  results  in  a  three-phase  field  between  R.  1.,. 
and  Lj.  The  boundaries  of  this  three-phase  field  are  the  lines  traversing  the  L,  +  L2 
miscibility  gap  illustrated  in  figure  2a  and  the  tie  lines  which  define  the  K-l  equilibria 
shown  in  figure  2c. 

The  calculation  of  the  isothermal  section  at  2100  K  can  be  concluded  by  adding 
the  L-BP  equilibria  as  shown  in  figure  2d.  In  this  case,  as  in  the  R-L  case,  there 
are  interactions  with  phase  pairs  which  have  been  considered  previously.  First  there 
is  an  interaction  between  L-BP  and  the  liquid  gap,  L,-L2.  as  illustrated  in  figure  2d. 
then  there  is  a  second  interaction  between  the  L-BP  and  L-P  pairs.  Comparison  of 
figures  2b  and  2c  with  2d  shows  that  while  the  L-R/L-BP  and  the  L-BP/L-P 
interaction  pairs  form  stable  three-phase  fields,  the  L-BP/L|-I  2  interaction  pair 
which  is  shown  in  figure  2d  is  metastable.  This  can  be  seen  by  comparing  figures  2c 
and  2d  and  noting  that  the  three-phase  BP  +  L)+  Lj  field  is  ‘covered’  by  the  R  +  L,+  L2 
and  R+L  fields.  Detailed  calculation  of  the  L-BP  equilibria  can  be  performed  by 
defining  the  free  energy  of  the  BP  phase  on  the  bases  of  tables  1  -4. 

For  the  case  of  the  SO-MO-SN  (BP)  phase,  where  the  compound  runs  from 
SN0  9*oSOo  om  lo  SO0.50M0.S0  and  from  SN0.9gMOQ.0j  to  MOq-soSOq-^q,  similar 
procedures  are  followed.  For  the  former  case, 

Gw  =  r°CsN  + jt°Gso +  >°C  mo+  (l  -  AGa  +  j^AC»+  AGk 

+  Irty  +  z  Inz-  (I  -  jc) In ( 1  —  x)J  .  (II) 

Here  x.  =0-50  and  x'.  =  0  02,  so  that  p  =  0-96,  x  =  x'.  +  yp  =  0  02  +  0 •  96>  =  vso. 
y  =  yMO>  and  ^  =  1  —  jc  — ^  =  0-98—  1  -96 y  =  zSN.  Moreover.  A CE  =  0  and  AGa 
and  A GB  can  be  explicitly  defined: 

AGa  =  (0  02)(0-98)(0- 98G^nso+0- 02GkosN-C)  = -820-0  082 T,  (12) 

AG„  =  (0-50)(0-50)(0-50G^,MO+0-50GUso-O  =  -26568+  10  503T:  (12) 

Gnp  =  (0-98-  I  96v)°GsN  +  (0  02  +  0-96v)oGso++°Gmo-  (l-yz^)(820 

+  0  0828T)-1~(26  568-  I0-503T)+ 8-3I4TL)  lny  +  (0-98 

-  1 -96y)ln(0  98-  1 -96i )- (0-98- 0-96r )ln(0-98- 0-96y  ))  .  (14) 

Thus  the  equilibrium  between  the  liquid  and  BP  phase  can  be  computed  as  in 
figure  2d.  Combination  of  the  pairwise  equilibria  in  figure  2  yields  the  computed 
isothermal  section  at  2100  K  shown  in  figure  3.  Isothermal  sections  similarly  derived 
at  2000  and  1900  K  arc  shown  in  figure  3.  along  with  an  observed  section  at  1500  K 
due  to  Midler  (1981 ).  The  latter  is  in  keeping  with  the  calculations  shown  at  1900  K. 
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Similar  calculations  were  performed  in  the  SO-YO-SN  and  SO-CE-SN  systems 
The  results,  which  are  complicated  by  the  existence  of  stable  quasiternary  compounds 
(listed  in  table  6),  are  displayed  in  figures  4-7.  The  SO-YO-SN  and  SO-CE-SN 
sections  were  computed  in  a  manner  similar  to  that  in  the  case  of  SO-MO-SN,  in 
that  the  systems  were  first  computed  as  if  the  ternary  phases  were  absent. 


<c)  (d) 

Figure  3.  Calculated  and  observed  isothermal  sections  in  the  SiOj-MgO-SijNa  system  (see  tables 
l-S  for  explanation  of  symbols),  (a)  2100  K;  (b)  2000  K;  (c)  1900  K;  (d)  observed,  1500  K, 
IV  =  SiOj-MgO  (SO0.»MO0.„)  (Muller  1981 ). 


Table  6.  Description  of  quasiternary  compounds  in  the  SiOj-YjOj-SijNt  and  Si0j-Ccj0j-SijN4 
systems*.  Sources  of  data:  GauVIer  et  al  (1980);  Lange  (1980). 


Compound  Base  Stoichiometry 

C  i(YSiOjN)  Y  SOo  isoYOo  sooSNo.jso 

D  i^YaSijC^Nj)  Y  SO0.100YO0t47SN0.j33 

E  VtlYjfSiOabN)  Y  SOojjjYOo  soSSN0  0sj 

Ci(CeSiOjN)  Y  SOo.1soOEosooSNo.350 

DA(Ce4Si,0,N,)  Y  SOo.100OEott7SNo.j3s 

E  jV(Ces(Si04)jN|  Y  SOo.37jCEo.59s  SN0  ■0B3 


C/J  (g  atom)"1  TJK  AG'(ii9»/J  (g  atom)' 1 

29096+  8-3687"  1953  -4393-2-4107" 

13004+16-7367"  2110  -678-4-1277" 
-2155  +  20-927"  1996  - 1 1  456- 0-8837" 

31  505+  8-3687"  1914  -5422- 2-4107" 

14769+  16-7367"  2020  -1259-4  1277" 
3975+  20-927"  2l09h -13585-0-8837" 


1  Tm,  melting  point;  see  tables  1 , 3, 4,  and  5  for  explanation  of  oilier  symbols. 
b  Decomposes  peritectically  into  Z  and  liquid  above  2 100  K. 


O 
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Subsequently  the  ternary  phases,  designated  as  C,  D,  and  F,  were  inserted.  In 
keeping  with  previous  practice,  the  free  energy  of  the  ternary  phases  C.  D.  and  K 
were  defined  by  choosing  a  base  phase  and  then  defining  a  compound  parameter  in 
conformity  with  experimental  observation.  In  the  present  case  the  Y  structure  was 
chosen  as  the  base  phase,  since  this  ternary  phase  composition  occurs  near  the  YO  and 
CE  corners  of  the  SO- YO-SN  and  SO-CE- YO  systems,  as  is  seen  in  figures  5  and  6. 

The  free  energy  of  the  ternary  compound  D  in  the  SO-YO-SN  system 
(D  =  ^  Y4Sij07N3  =  SO(j.  iooY00.<67SNo.jj3)  is  defined  as 

GD  =  O'  l0°Gso+0-667oGyO4- 0-23^°Gsn+  (0  - 1  OHO  667)  ‘<0  1  Gfeovo 

+  0-667GVoso-O+<010M0  233X0  333)  '((MGio.iv,  *  0-233GlsNS(,  <> 
+  (0-667)(0-233)(0  Wr1(0-667G,YosN  +  0-223G^NYo- <  >  (15 

Since  C  =  13004+  16-736 T,  for  this  phase, 

Gn  =  0-10°Gjo  +  0-667°G?o+0- 233°G£n-  9094-  3-0O2T  .  do 

Thus  the  free  energy  of  formation  of  f)  from  the  Y  form  of  SO.  YO.  and  SN  is  given 
by  -9094- 3-0627".  The  free  energy  of  formation  of  D  from  the  stable  forms  of 
SO,  YO,  and  SN  (ie  R,  Y,  and  B.  respectively)  can  be  computed  by  using  the  lattice 


The  lattice  stabilities  from  Ref.  9 


EGbcc,dis  =  xx  (-32345-37 .92T)X_  +  (-150252+51. 42T)X..  J/mol 

m  re  A1  Fe  A1 

W1  =  1588,  w2  -  794  (k  units  k-13.8‘10-24  J/atom),  K-0.68 

v  frr 

m  =  XFeXAl  (~25022“38 •01T>xFe  ♦  (-106841+66 . 70T)XA1  J/mol 

ECL  -  X_  XAl  (-62760-14. 63T)X-.  +  (-96232+32. 72T)X..  J/mol 

m  Fe  Al  Fe  A1 

°G®  ..  -  0.25°GpCC-0.75  °GA“  -  -35275+8. 28T  J/mol 

Fe0.25A10.75  Fe  A1 

°GG  ..  -  0.286  °G^C-0.714  °G^C  -  -41169+11 .40T  J/mol 

Fe0.286A10.714  Fe  A1 

d  F 

r,  °  bre  °  hrr 

Ve0.333A10.667  "  0,333  GFe  "°-667  Ga1  "  -37872+4. 73T  J/mol 

oji  ..  -  0.4  GA*c-0.6  °clCC  -  -23863-2. 77T  J/mol 

Fe0.40M0.60  Al  te 
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A  THERMODYNAMIC  EVALUATION  OF  THE 


TITANIUM-CARBON-NITROGEN  PHASE  DIAGRAM. 


John  Agren 


1 . 0  Introduction 

The  use  of  vapor  deposited  TiC,  TiN  and  Ti(C,N) 
coatings  for  improving  the  performance  of  cutting  tools  and  wear 
resistant  surfaces  enhances  the  value  of  information  on  the 
thermochemistry  and  phase  diagram  of  the  Ti-C-N  system. 

Since  the  crystal  structures  of  the  compounds  TiC  and 
TiN  are  the  same,  it  is  likely  that  these  compounds  form  a  solid 
solution.  However,  the  experimental  Ti-C  and  Ti-N  phase  diagrams 
reported  to  date  arc  not  very  well  established  and  very  little 
thermodynamic  data  is  available  for  the  binary  cases  or  the 
Ti-C-N  ternary.  The  purpose  of  the  present  work  is  to  evaluate 
all  the  information  available  in  terms  of  thermodynamic  models 
for  the  individual  phases  and  synthesize  the  binary  and  ternary 
phase  diagrams. 

2 . 0  Thermodynamic  Models 
2 . 1  Liquid 

The  simplest  models  that  cover  the  whole  concentration 
range  are  the  regular  or  subregular  solution  models.  The  formation 
of  the  solid  compounds  TiC  and  TiN  in  the  binaries  suggests  that 
some  ordering  or  formation  of  molecular  species  might  take  place 
in  the  liquid  around  the  stoichiometric  compositions.  However, 
a  high  degree  of  association  is  not  observed  in  the  gas  phase  and 
vaporization  of  titanium  carbide  yields  titanium  and  carbon  polymer 
The  Fe-S  system  has  recently  been  treated  by  assuming  the  presence 
of  molecular  species  (1)  and  by  assuming  ordering  (2).  Both 
treatments  yield  a  satisfactory  agreement  with  experimental  data, 
and  predict  a  very  rapid  change  in  activity  close  to  the 
stoichiometric  composition.  This  rapid  change  is  difficult  to 
describe  with  a  regular  or  subregular  solution  model,  and  generally 
many  terms  will  be  needed  in  the  series  expansion  of  the  excess 


F  (SiOj-MgO)  and  (MgSiNj)  above  446  K.  Thus  the  sections  shown  in  figure  3 
represent  metastable  equilibria. 


Figure  7.  SiO,-CCjOj-SijN«  system:  (a),  (b)  calculated  pairwise  equilibria  at  1900  K; 

(c)  calculated  isothermal  section  at  1900  K;  (d)  observed  (Lange  1980)  subsolidus  equilibria  at 
1923-2023  K  (S  and  D  phases  were  expected  but  not  observed). 
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IV  CALCULATION  OF  TERNARY  III  -V  and  II-IV  SYSTEMS 


L.  Kaufman,  J.  Nell,  K.  Taylor  and  F.  Hayes,  CALPIIAD  (1981)  5 


(This  paper  was  presented  at  Calphad  X,  Vienna,  Austria  .July  1  ‘,*8 1 1 

ABSTRACT.  A  data  base  covering  the  binary  systems  composed  of  Aluminum,  Gallium,  Indium, 
Phosphorus,  Arsenic  and  Antimony  has  been  constructed  by  analyzing  the  fifteen 
combinations  of  these  elements  in  terms  of  lattice  stability,  solution  phase  and 
compound  parameters.  Partial  isothermal  sections  in  the  P-In-As,  As-tn-Sb,  I’-Ga  As. 
Ga-Sb-In  and  Al-Sb-Ca  systems  were  then  calculated  using  the  foregoing  data  base 
for  comparison  with  experimental  isothermal  sections  and  quasi-binary  III -V  phase 
diagrams.  It  was  found  that  ternary  liquid  and  I1I-V  compound  interaction  para¬ 
meters  were  required  to  attain  good  agreement  in  some  cases.  Similar  calculations 
were  performed  for  the  Te-Cd,  Hg-Cd  and  Te-llg  binary  systems  and  the  Cd-le-Mg 
ternary  systems  at  pressures  up  to  74  atmospheres.  Comparison  of  the  calculated 
results  with  experimental  data  on  tie-line  compositions  between  Cd-Te-Hg  liquid 
and  quasi-binary  CdTe-HgTe  alloys  is  important  in  the  liquid  phase  epitaxial  growth 
of  controlled  band  gap  electro-optical  materials. 

1 .  Introduction 

Alloy  semiconductors  are  Important  for  a  wide  range  of  electro-optical  applications 
because  their  properties  can  be  tailored  by  controlling  the  composition  of  the  solid,  which  is 
grown  from  the  liquid  or  deposited  from  the  vapor.  Although  crystal  growers  have  become  very 
adept  at  controlling  the  temperature  (and  pressure)  of  the  parent  phase  in  order  to  obtain 
the  desired  characteristics  of  the  crystalline  alloy  semiconductors,  phase  diagram  data  is 
indispensable  in  order  to  deal  with  the  increasingly  complex  systems.  Stringfellow  and 
co-workers  (1-6)  have  pioneered  in  developing  simple,  useful  models  for  calculating  quasi - 
binary  systems  based  on  III -V  compounds.  Stringfel low’ s  pioneering  work  has  been  extended 
and  expanded  by  Ansara  and  co-workers  (7),  Osamura  and  Murakami  and  co-workers  (8-10)(  and 
Brebrlck  (11)  who  have  used  more  extensive  experimental  thermochemical  data  and  less  restrict¬ 
ive  models  than  the  regular  solution  (and  quasichemical)  model  employed  by  Stringfellow. 

The  current  paper  attempts  to  assemble  a  data  base  for  III -V  compounds  which  can  be  employed 
to  compute  multicomponent  phase  diagrams  over  a  wide  range  of  temperatures (and  pressure?  )  both 
outside  and  within  the  quasi-binary  plane.  In  addition,  the  Cd-Hg,  Cd-Te  and  Hg-Te  binary 
systems  are  analyzed  and  combined  to  calculate  the  Cd-Te-Hg  ternary  system  over  a  range  of 
temperatures  (and  pressures)  in  order  to  compute  tic-line  compositions  between  quasi-binary 
CdTc-HgTa  alloys  and  the  Cd-Te-Hg  liquid.  The  former  alloys  are  of  current  interest  in 
synthesizing  controlled  band  gap  electro-optical  materials'.  In  the  analyses  of  both  the 
III-V  and  I I -VI  (Cd-Te-Hg)  systems.the  methods  and  formalism  which  are  applied  are  those 
employed  in  dealing  with  metals  and  metallic  alloys  ( 12, 13) , sil icides  (14), and  oxynitride 
systems  (15). 

2.  Thermochemical  System  Employed  to  Characterize  Binary  111-V  Phase  Diagrams 

The  method  employed  for  describing  solution  and  compound  phases  is  the  same  as  that 
employed  to  describe  a  variety  of  metal,  metalloid,  and  oxide  systems  (13-1S)  incorporating 


some  sybolic  usage^which  facilitates  data  handling  as  indicated  below.  The  free  energy  of 
the  liquid  phase, G  ,  in  the  binary  system  Al-In  is  given  by  equation  (1)  where  T  is  in 
Kelvins,  R=8.314  J/g.at?K,  and  x  is  the  atomic  fraction  of  Indium: 

GL  -  (l-x)*G^j»x*Cjn*RT(xfnx*(l-x)fn(l-x))  *x(l-x)  [  (l-x)LALlN^xUNAL) .  (1) 

Similarly,  the  free  energy  of  the  fee  phase  in  the  Al-tn  is  defined  by  Equation  (2)  as 

GfCC  *  (l-x)*G^jC*x*Gj^C*RT(xfnxx(l -x)  £n( l -x)>x( 1 -x) ( ( 1-x) AALlN*xAlNAL)  (2) 


where  the  difference  in  free  energy  between  fee  and  liquid  aluminum, ALALIA,  and  the  difference 
in  free  energy  between  fee  and  liquid  indium,  ININLA,  are  listed  in  Table  I,  while  the 
solution  parameters  LALIN,  LINAL,  AALIN,  and  A1NAL  are  given  in  Table  2.  The  liquid  para¬ 
meters  for  the  In-As  and  Ga-As  systems  were  taken  directly  from  Brcbrick's  analysis  (11). 
Brebrick  also  provides  parameters  for  the  ln-Sb  and  Ga-Sb  systems,  namely  LINSB=-3431 - 19 . 782T, 
LSBIN* -27698*1 9. S8 IT,  LGASB-21652-32. I03T  and  LSBCA=-30020»28. 259T  J./g.at.  While  these 
parameters  provide  a  good  description  of  the  ln-Sb  and  Ga-Sb  phase  diagrams  and  the  excess 
free  energy  of  mixing  for  the  liquid  phases  in  these  systems, they  do  not  describe  the 
enthalpy  of  mixing  for  the  liquid  phase  in  the  ln-Sb  and  Ga-Sb  systems  as  measured  by  Ansara 
et  al..  (7).  Thus  in  the  case  of  the  ln-Sb  system,  the  enthalpy  of  mixing  defined  by  Brcbrick's 
LINSB  and  LSBIN  parameters  is  given  by  Equation  (3)as 

EHL  *  -x(l-x)  [(l-x)3431*27698x]  J/g.at.  (3) 

where  x*Xg^.  By  contrast,  the.  equation  for  EHE  provided  by  Ansara  et^  a_h  (7)  for  ln-Sb,  which 
is  based  upon  experimental  measurements  of  the  heat  of  mixing,  is  given  by 

EHL  *  -x(l-x)(U689«19774x-46661x2*24980x3)  J/g.at.  (4) 

Ansara  et  ahalso  provide  the  following  equation  for  the  excess  entropy  of  mixing  of  the 
1 iquid  as 

ESL  •  x(l-x)[4.602-7,799x+12.104x2-6.782x3]  J./g.at.°K  (S) 

Although  the  excess  free  energy  provided  by  Brebrick* s  parameters  and  those  generated  by  Eqs. 
(4)and(S)in  the  temperature  range  of  interest  (near  700'K)are  in  relatively  good  agreement, 
the  enthalpy  of  mixing  described  by  Equations (3)and  (4)differ  by  800-1600  J/g.at.  Consequently, 
Equations  (4)and  (5)were  employed  to  derive  LINSB*-1S380-3.607T  and  LSBIN=-1029V2.393T  shown 
in  Table  2.  These  expressions  provide  a  closer  fit  to  Equations  (4  land  (5 )than  Brebrick* s 
parameters  and  still  yield  a  satisfactory  agreement  in  the  calculated  and  observed  ln-Sb 
phase  diagram.  A  similar  comparison  can  be  made  for  the  Ga-Sb  system.  In  the  latter  case 
Brebrick’s  analysis  yielded  LGASB=21652-32. 108T  and  LSBGA=-30020«28. 2S9T  J./g.at.  By  contrast, 
Ansara  et  al_.  (7)  suggest  that 

EtlL  »  x(  1-x)  {531-21599x*3217Sx2-  15788x3}  J/g.at.  (6) 

where  x=x<-b  and 

ESL  *  x(l-x)[4. 849-14. 456x»24.S27x2-1S.786x3]  J/g.at.aK  (7) 

The  present  parameters,  LGASB*4962-3. 209T  and  LSBCA=-9715-0.456T  provide  a  better  description 
of  the  experimental  heat  of  mixing  for  Ga-Sb  alloys  embodied  (12)  by  Equation  (5)  than  do  the 
parameters  suggested  by  Brebrick.  The  parametric  descriptions  of  the  In-Al ,  Ga-Al,  and  ln-(.a 
systems  shown  in  Tables  1  and  2  were  derived  from  the  earlier  description  provided  by 
Ansara  et_  al.  (12).  The  remaining  111 -V  systems  were  characterized  as  shown  in  Tables  1-3 
based  on  the  available  thermochemical  and  phase  diagram  data  employing  the  standard  compila¬ 
tions  of  such  information  (16-20)  and  available  lattice  stability  data.  Table  4  compares 
the  calculated  and  observed  thermochemical  properties  of  III -V  compounds  based  on  this 
description  while  Figures  1-14  display  the  partial  phase  diagrams  computed  with  the  current 
description.  The  heat  of  mixing  of  liquid  At-  -Sb-  -  at  1400*K  is  given  as  -2477  J./g.at. 
based  on  experimental  data  (17)  and  -1126  J./g.at.  ofrthe  basis  of  Table  2. 


The  free  energy  of  the  ternary  liquid  phase  in  the  P- In-As  system  is  defined  on  the 
basis  of  Kohler’s  equation  as 

GL  =  (l-x-y)°GE+x°Gf  E»y°(',Es+P.T((l-x-y)£n  (l-x-y)»x£nx+y£ny) 

+(l-x-y)x(l-y)'1((l-x-y)LPPIN*xLlNPP]  +  (l-x-y)y (1 -x) _1 ( (1-x-y) LPPAS*yLASPP) 

♦  xy (x*y) ” * [xLINAS+yLASIN]  ♦  xy(l-x-y)  TRNL  J./g.at.  (8) 


where  x  is  the  atom  fraction  of  In,  y  is  the  atom  fraction  of  As  and  TRNL  is  the  ternary  inter¬ 
action  parameter.  Values  of  TRNL  for  five  ternary  systems  which  have  been  investigated  here 
are  listed  in  Table  5.  The  specific  value  for  the  P-In-As  system  is  12S52  J/g.at.  This  is  a 
rather  small  value  since  the  maximum  value  of  the  ternary  term  occurs  at  x*y« (l-x-y)=l/3.  Thus 
the  maximum  contribution  of  this  term  is  only  46S  JL/g.at.  If  attention  were  restricted  to  the 
quasi-binary  join  (i.e.,  where  x=0.5)  between  In_  .P.  .  And  In.  -As.  .  the  maximum  contribution 
is  390  J/g.at.  when  x=0.S,  y=0.2S  and  (l-x-y>0 . 25.  0-5  0,5  0,5 

The  free  energy  of  the  In.  CP  -In,.  .As„  -  zinc-blende  phase  (S)  is  defined  by  Equa¬ 
tion  (9)  as  follows  0,5  °*5  °*5  °-5 
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TABLE  1 


I 

I 


ALALLA  -  10711  -  11.S06T* 
ALALLO  «  0  -  18.451T 

GAGALO  -  5590  -  18. 45 IT 
GAGALA  =  -5021  -  8.368T 
GAGALT  ■  -1674  -  8.368T 

ININIT  *  3264  -  7.594T 
ININLA  -  2908  -  7.113T 
ININLO  «  0  -  18.4S1T 


*  ALALLA  -  c“quid  -  cf'C 
A1  A1 


A»fcc*,  T»tetragonal ,  O=orthorhombic, 
phosphorus  structure,  W*white  phosphorus  structure 


PPPPLX 

B 

17991  - 

21.506T 

PPPPWX 

S* 

18171  - 

19.539T 

PPPPLA 

= 

-18828  - 

8 . 3fa8T 

PPPPLR 

K 

4184  - 

23.012T 

SBSBLR 

s 

19874  - 

21.966T 

SBSBLA 

= 

-18828  - 

8.368T 

ASASLR 

- 

24874  - 

23.012T 

ASASLA 

B 

-18828  - 

B.368T 

ASASLX 

= 

21631  - 

21 .63  IT 

10711  -  11.506  T  J./g.at. 


SWWARY  OF  LATTICE  STABILITY  PARAMETERS 
(All  units  in  Joules  per  gram-atom  (mole  of  atoms),  T  in  Kelvins) 

L» liquid* 

R*rhoabohedral ,  X*red 


TABLE  2 

SUKWARY  OF  SOLUTION  PHASE  PARAMETERS* 

(All  units  in  Joules  per  gram-atom  (mole  of  atoms),  T  in  Kelvins) 


•for  Ai-In  system,  L«Iiquid,  A<=fcc,  x=*ln 
GL-(1-*)gJ^  ♦  xGjn  ♦  RT(jc£nx  ♦  (l-x)fn(l-x))  ♦  x(l-x)  ((l-x)LALIN  ♦  xLlNAL] 


-fee 

G  a 

C1-*)G*C  ♦  xG^c  ♦  RT(*tnx  ♦ 

(l-x)£n(l-x))  +  x(l- 

x)  [(l-x)AALIN  ♦  xAINAL] 

LALAS 

-16317  -  41.338T 

LALIN 

■ 

25522 

-  3.347T 

lasal 

-16317  -  41.338T 

LINAL 

e 

20502 

-  2.510T 

LALPP 

AALIN 

8 

60250 

-  3.347T 

-38911  -  26.275T 

AINAL 

8 

55229 

-  2.510T 

LPPAL 

-38911  -  26.275T 

LINPP 

-29288  ♦  5.188T 

LALSB 

8050 

9.706T 

LSBAL 

-17054  -  9.706T 

LPPIN 

-79496  ♦  S.188T 

RINPP 

8368 

LSBAS 

LASSB 

«  0 

RPPIN 

8368 

RSBAS 

RASSB 

=  8368 

LGAPP 

-90374  ♦  34.225T 

LSBPP 

LPPSB 

'  -15899 

LPPGA 

-58576  -  49.036T 

RSBPP 

RPPSB 

■  -20920-4. 184T 

XGAPP 

8368 

XSBPP 

XPPSB 

•  -15899 

XPPGA 

8368 

LASPP 

LPPAS 

-  -15899 

LGASB 

4962  -  3.209T 

RASPP 

RPPAS 

*  -20920-4 . 184T 

LSBGA 

-971S  -  0.4S6T 

XASPP 

8 

XPPAS 

*  -15899 

RSBGA 

RGASB  •  8368 

LGAAS 

-1845 

-  21.209T 

LALGA 

3280  -  2.021T 

LASGA 

-41551 

♦  7.192T 

LGAAL 

2071  -  0.74ST 

RGAAS 

RASGA 

*  8368 

AALGA 

-1046  ♦  5.439T 

AGAAL 

-20962  ♦  26.778T 

LINAS 

-243  - 

29.S22T 

OGAAL 

OALGA  ■  8368 

USIN 

-48124 

♦  17.472T 

TASIN 

TINAS 

•  8368 

RASIN 

RINAS 

*  8368 

TABLE  2  (CONTINUED) 


SUMMARY  OF  SOLUTION  PHASE  PARAMETERS 
(All  units  in  Joules  per  gran-atom  (mole  of  atoms)  T  in  Kelvins) 


LINSB  •  -1S380  -  3.607T 
LSBIN  .  -10293  -  2.393T 
TSB1N  •  TINSB  •  8368 
RSBIN  «  RINSB  «  8368 


LGAIN  =  4435  ♦  2.268T 
L1NGA  -  4435  ♦  0.9S4T 
TGAIN  *  -3431  ♦  21.907T 
TINGA  •  -3376  ♦  24.393T 
OCA IN  •  OINGA  =  12552 


TABLE  3 

COMPOUND  PARAMETERS  FOR  III-V  ZINC-BLENDE  COMPOUNDS 
q  fee  fee 

G  «  0.5  g‘  ♦  0.5  G*  ♦  0.25  [0.5  LI1JJ  ♦  0.5  LJJII  -  C]  Joules/g.at. 


Compound 


Compound  Parameter  (C)  Compound 
(Joules/g.at.  *K) 


Compound  Parameter  (C) 
(Joules/g.at.  T) 


^o.s^.s 

364 

,443  -  60.  51 7T 

^O.S^O.S 

179,343  - 

77.S71T 

^0.5AS0.S 

305 

,348  -  78.994T 

In0.5P0.5 

188,749  - 

64.099T 

^0.SSb0.5 

174.! 

992  -  SS.81ST 

Ino.sAso.s 

187,217  - 

73. 806T 

^o.s^.s 

280 

,328  -  75.647T 

,n0.5Sb0.S 

125,227  - 

63. OUT 

# 

GaQ  5As0  5  259. S25  -  83.496T 

TABLE  4 

COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  THERMOCHEMICAL  PROPERTIES  OF  III 

-V  COMPOUNDS 

Compound 

(J./mol 

AHf[298*K] 

.,  J./mol . *K) 

AS  [298'K) 

AS 

[fusion] 

# 

calculated 

exptl.  calculated  expt 1.(16) 

calculated 

exptl. 

AlP 

-164,849 

-164,431*2929(16) 

-3.97  -3.85*5.0 

61.66 

61.50*1.6(16) 

MAs 

-117,127 

-122,591*5021(16) 

-4.18  -3.76*4.8 

70.91 

70.42(43) 

MSb 

-51,045 

-5^208*418(16) 

-9.46  -8.86*0.8 

59.32 

61.67(43) 

c 

GaP 

-129,971 

-122^173*8368(16) 

-10.92  -11.51*  1.9 

66.10 

66. 53(43) 

GaAs 

-8^291 

-88,575(11) 

-13.51  -12.47*1 . 3 

70.015 

69.63(11) 

CaSb 

-41,547 

-41,840*1255(16) 

-14.19  -9.20*1.9 

67.061 

67.31(11) 

InP 

-84,391 

-75,312*8368(16) 

-21.02  -20.88*1.9 

59.07 

59.08(43) 

InAs 

-61,647 

-61^23(11) 

-18.74  -18.83*2.5 

63.92 

63.36*6.89(16) 

InSb 

-2^974 

-30,962(11) 

-15.93  -15.69*1.1 

58.53 

59.82(11) 

-32- 
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Figure  J.  Calculated  Al-Ca  Phase  Diagram 


U 


I 


Met ast able! 

_ Note:  As  valor¬ 
izes  at  87(i°K 
and  1  atmos¬ 
phere 


Figure  4. 


Calculated  In  -  As  Phase  Diagram 


Figure  7. 


Calculated  Sb-As  Phase  Diagram 


Figure  8.'  Calculated  Ga-Sb  Phase  Diagram 


•v 


V. 


GS  =  (0.5-y)°GpCC  ♦O.S°G^C*y°G^C  ♦(l-Zy)  (0.25)  (0.S  Ln* I N+ 0 . 5 L r NFP -C j ) 

♦2y(0.2S) (0.5LINAS»0.SLASIN-C2)  *  CAB  y(l-2y) 

♦RT(yiny*(0.S-y)tn(0.5-y)-0.Sln0.S)  J./g.at.  (9) 

where  Cj,  the  compound  parameter  for  ln^  gP^  g  =  188,749-64 . 099T  J^g.at.  and  C^,  the  compound 
parameter  for  In^  gAs^  g=187, 217  -  73.806T  J. /g .at., as  given  in  Table  3.  The  interaction  para¬ 
meter  for  mixing  lnQ  gP^  g  and  InQ  gAsQ  g,  CAB,  is  equal  to  zero  in  this  case,  as  shown  in 

Table  S.  Thus  this  compound  is  defined  as  an  ideal  solution  of  its  components.  Since  the  free 
energies  of  the  liquid  and  zinc-blende  phases  can  be  defined  explicitly,  the  partial  ternary 
phase  diagrams  covering  the  L-S  equilibria  can  be  computed  as  a  function  of  temperature.  This 
procedure  has  been  carried  out  for  the  ternary  systems  listed  in  Table  5.  The  results,  dis¬ 
played  as  isothermal  sections  or  as  quasi-binary  joins. arc  compared  with  experimental  findings 
and  Stringfel low' s  calculations  (1)  in  Figures  15-2S.  A  calculated  isothermal  section  due  to 
Ansara  et  al.  (7)  for  the  Sb-Ga-In  case  at  823°Kis  included  for  comparison  in  Figure  21.  The 
latter  results  are  in  good  agreement  with  the  current  findings  shown  in  Figure  2  3.  Although 
it  is  not  possible  to  make  an  exact  comparison  of  the  interaction  parameters  employed  by 
Stringfellow  in  his  quasi-binary  calculations  (1,6)  and  those  employed  here, an  approximate  com¬ 
parison  can  be  made  in  the  following  way.  In  Stringfcllow’s  study. the  quasi-binary  system 
InP-lnAs  was  treated  as  if  the  liquid  and  solid  phases  were  both  regular  solutions.  The 
interaction  parameter  for  the  liquid  was  taken  to  be  5238  J^mol.  Since  the  heat  of  mixing  is 
the  product  of  the  mol.  fraction  and  the  interaction  parameter,  the  maximum  value  of  the  heat  of 
mixing  is  (0. 25) (S238) =1309  J/mol=655  J/g.at.  for  liquid  O.SInQ  gPu  g-0..1n0  gAsQ  g.  This 

value  is  listed  in  Table  6  along  with  comparable  parameters  for  the  other  systems  of  current 
interest  cited  by  Stringfellow.  The  latter  results  are  compared  with  those  derived  from  the 
current  study  on  the  basis  of  Equations  (8)  and  (9)  and  Table  S.  As  indicated  above, the  maxi¬ 
mum  contribution  along  the  quasi-binary  join  is  equal  to  TRNL/32  and  CAB/8.  Table  6  lists 
these  values  for  comparison  with  Stringfel low' s  results  in  Table  6.  While  there  is  no  obvious 
correlation  between  the  current  results  and  Stringfel  low*  s  values,  aU_  of  the  values  listed  arc 
small!  The  analysis  of  the  Ga-Sb-In  system  carried  out  by  Ansara  £t  al.  (7)  utilized  a  ternary 
heat  of  mixing  equal  to 


EHL  *  xInxGaXSb  [23234X.  -24066xr  ♦S82xc.J.  (10) 

In  (j3  Sd 

The  present  analysis  sets 

V-*  -xInXGaxSb  8368  •  <1') 

At  x.  =x_a=x  =0.333, Equation  (10)  yields -3. 1  J/g.at . , whi le  Equation  (11)  yields  -309  J./g.at. 
At  compositions  near  the  edges  of  the  ternary  (P.1,0.4,9.5)  Eq.  (11)  yields  -167  J/g.at ., whi  le 
Eq-  (10)  yields  results  between  +188  and  -188  J.g./at.  All  of  these  contributions  are  small. 


TABLE  5 

COMPILATION  OF  TERNARY  PARAMETERS  FOR  III-V  LIQUID  AND  ZINC-BLENDE 

COMPOUND  PHASES 


(All  Units 

in  Joules  per  gram 

-atom  (mole  of  atoms) 

T  in  Kelvins) 

System 

TRNL 

(Joulcs/g.at.) 

Zinc-blondc  Phase 

CAB 

(Joulcs/g.  a 

P-ln-As 

12552 

In0.5^P,As^0.S 

0 

As-In-Sb 

0 

In0.5(As,Sb^0.5 

9623 

P-Ga-As 

-31380 

Ga0.5^P,AS^0.5 

2092 

Ga-Sb- In 

-8368 

(Ga,ln)0  gSb0  5 

4602 

Al-Ga-Sb 

-8368 

(Al.Ca)0  5Sb0i5 

-4184 

o 
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TABU  i) 

COMPARISON  OF  THE  MAXIMUM  QUASI-BINARY  HEAT  OF  MIXING  FOR  Till  LIQUID  AND  SOI  ID  COMPOSITIONS 
FROM  STRINGFELLOW  (1,6)  AND  THE  CURRENT  ANALYSIS  (. Ionics/*. at .  ) 


System 

LIQUID 

Stringfellow(l) 

Current 

SOLID 

St  ringfel low  Current 

In.  p  -In„  ,  As„  , 
0.5  0.5  0.5  0.5 

655 

392 

306 

0 

,nO.SAsO.S‘InO.SSbO.S 

244 

0 

1197 

\202 

GaO.SP0.5'GaO.SAsO.S 

655 

-980 

515 

20  0 

G'VsSb0.5'InO.SSb0.5 

1409 

-260 

965 

•T  1 

AVsSbo.s'Gao.ssbo.s 

664 

-260 

12 

-S2  t 

Thus  the  current  analysis  suggests  a  simple  method  for  calculating  isothermal 
sections  for  III-IV  ternary  systems.  In  particular,  if  3  given  section  is  required,  the  d  it  \ 
base  contained  in  Tables  1-3  can  be  employed  to  calculate  the  section  of  interest  with  1KM.« 
CAB=0.  The  process  can  then  he  repeated  employing  a  range  of  TRNI.  and  CAB  similar  to  those 
shown  in  Table  S  to  predict  the  spread  of  expected  results.  If  a  few  experimental  data  are 
available  at  a  given  temperature,  they  can  be  readily  employed  to  fix  TRNI.  and  CAB,  which  in 
turn  can  be  employed  to  predict  the  behavior  at  other  temperatures.  It  is  hoped  that  future 
studies  will  permit  improvement  of  the  current  data  base. 

3.  Thermochemical  System  Employed  to  Characterize  Binary  11-V1  Phase  Diagrams 

The  II- IV  compounds  of  Cd-  Te.  and  Hg  PTc  .  form  a  scries  of  cubic  sine- blende 
solid  solutions  which  offer  the  'attrActivi'  property  of  avariahlc  band  gap  that  is  a  function 
of  composition,  the  (Cd,Hg)0  5TeQ  ,-solid  solution  displays  band  gaps  which  run  in  a  nearly 
linear  fashion  between  the  wide  gap  semiconductor  CdQ  5TeQ  (E^=1.6  eV)  and  the  semi- 
metallic  compound  Hg0  ,.Te0  ,  which  can  be  considered  to  be ’a  semiconductor  with  a  negative 

band  gap  (37).  Alloys  of  tdQ  jHg0  4Te0  5  and  Cd0  2HgQ  jTe0  5>  exhibiting  band  gaps  ranging 
from  0.1  to  0.5  eV,  respectively,  are  of  particular  interest  for  detection  of  infrared  radia¬ 
tion  in  the  2  to  20  micron  range.  The  growth  of  such  alloy  compositions  from  the  melt  presents 
some  difficulties  due  to  the  "width"  of  the  two  phase  liquid  plus  solid  field  along  the 


^0.5^e0.5 "^®0  5^e0  5  <luasx-binary  join.  Liquid  phase  epitaxial  (LPE)  growth  methods  have 
been  employed  to  circumvent  this  problem  (37),  but  controlled  growth  requires  a  knowledge  of 
the  tie-lines  between  the  alloy  compound  and  the  liquid.  Such  information  can  be  generated 
experimentally  with  considerable  effort.  Alternatively,  It  is  possible  to  use  the  CAI.PIIAD 
method,  which  has  been  applied  to  a  very  wide  range  of  systems,  for  calculating  the  ternary 
Cd-Te-Hg  system  from  a  knowledge  of  the  component  binary  systems  in  order  to  provide  the 
required  information  on  tie-line  composition.  The  description  provided  below  provides  an 
account  of  such  calculations. 

4.  Description  of  the  Vapor,  Liquid  and  Zinc-blende  Phases 

Table  7  summarizes  the  current  descriptions  of  the  liquid(L),  vapor(V)  and  zinc- 
blende,  ),  phases  in  the  cadmium-tellurium-mercury  system  as  a  function  of  temperature  (T,°K), 
pressure  (P,  atmospheres),  atomic  fraction  tellurium(x) ,  and  atomic  fraction  of  mercury(y). 
These  equations,  along  with  the  associated  lattice  stability  values  relating  the  free  energy 
of  the  vapor,  liquid,  fee,  and  stable  forms  of  cadmium,  tellurium  and  mercury, have  been  assem¬ 
bled  by  employing  existing  (38,14-17)  thermochemical  and  phase  diagram  data  for  the  pure 
elements  and  binary  systems.  Figures  26-31  show  the  partial  binary  phase  diagrams  computed 
for  the  Te-Cd,  Hg-Cd  and  Te-Hg  systems  at  P=1  atmosphere  and  P*16  atmospheres.  These  phase 
diagrams  follow  directly  from  the  ternary  equations  listed  in  Table  7  when  they  are  reduced 
to  binary  systems  for  each  of  the  edge  binary  systems  in  question.  The  small  ternary  terms 
assigned  to  the  liquid  phase  (-33472xy(l-x-y))  and  the  zinc-blende  phase  (-4l80y  (t-2y)) 
Joules/g.at.  were  chosen  in  conformity  with  the  experimental  results  along  the  Cd^  sTe()  j- 

Hgjj  jTSq  5  quasi-binary  join  shown  in  Figures  32  and  33.  Since  the  maximum  value  of  these 


terms  is  attained  at  x*0.5  and  y= 0.25,  it  is  apparent  that  the  ternary  correction  term  is  at 
most  about  -1000  Joules/g.at.  for  the  liquid  phase  and  -500  Joules/g.at.  for  the  zinc-blende 


phase  at  the  x*0.5,  y«0.25  composition.  Table  8  compares  calculated  and  observed  thcrmochcmi- 
cal  values  for  Cdg^Tegs  and  Hgg  jTeg.s  at  298  K. 


Figures  34-36  display  calculated  isothermal  sections  in  the  Cd-Te-Hg  systems  between 
I  and  74  atmospheres  at  temperatures  ranging  from  773°K  to  12I3°K.  These  calculated  iso¬ 
thermal  sections  show  the  tie-lines  between  the  zinc-blcnde  and  liquid  phases  in  addition  to 
the  three  phase  vapor  /liquid/zinc-blende  fields  which  arc  especially  important  in  crystal 
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•  Present  Calculation 

Stringfcllow  Calculation  (1) 
O  ref.  (22) 

O  ref.  (23) 

♦  ref.  (24) 
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Figure  16. 
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Miscibility  Gap  (calculated) 


In0.5As0.5 


Calculated  and  Observed  equilibria  Along  the 
*n0.5Sb0.5  "  In0  5As0  5  Quasi_binary  Join 
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Experimental  Points:  •  ref.  (JO);  q  ref.  (28)  .  A 


growth  operations.  The  calculations  of  tie-lines  can  be  compared  directly  with  the  recent 
experimental  results  of  Mroczkowski  and  Vydyanath,  who  determined  tie-line  compositions  by 
means  of  "closed  tube  tipping  experiments".  In  these  experiments,  crystals  of  the  zinc-blende 
phase  of  fixed  composition  are  equilibrated  with  the  liquid  phase  in  a  graphite  boat.  After 
equilibration,  the  temperature  was  decreased  by  S  to  7 °C  and  the  liquid  was  decanted.  Subse¬ 
quent  electron  microprobe  analysis  of  the  concentration  profiles  along  the  length  of  the 
crystal  (from  the  zinc-blende  single  crystal  phase  into  the  polycrystalline  phase,  which  forms 
when  the  contact  liquid  freezes)  permits  determination  of  the  tie-line  compositions  (37).  At 
818°k,  the  following  comparison  can  be  made  between  the  calculated  and  observed  tic  lines. 


Zinc 

-blende 

Liquid 

Atomic  Percent 

Hi 

Cd 

tin 

Cd 

Observed  (37) 

32.5 

17.5 

20.3 

2.3 

Calculated 

32,5 

17.5 

19. S 

2.0 

This  excellent  comparison  provides  an  independent  check  of  the  description  of  the 
liquid  and  zinc-blende  free  energies  shown  in  Table  7.  These  equations  can  be  used  to  com¬ 
pute  the  liquid/vapor/zinc-blende  equilibria  over  a  wide  range  of  temperatures  and  pressures, 
which,  in  turn,  can  be  used  to  define  the  conditions  for  stable  conventional  and  LPE  growth  of 

Te-rich  and  Te-poor  (Cd,Hg).  cTe„  .  .  , 

'  1  s  0.5  0.5  compositions. 

5.  Conclusions 

The  examples  of  calculated  III - V  and  1 1 -VI  multi -component  systems  presented  here 
illustrates  the  scope  and  utility  of  the  method  presented.  The  range  of  systems  to  which  such 
computations  can  be  applied  will  expand  as  additional  binary  systems  are  added  to  the  data  base. 
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TABLE  7 


SUMMARY  OF  FREE  ENERGY  EQUATIONS  DESCRIBING  TOE  VAPOR 
LIQUID, AND  ZINC-3LENDE  PHASES  IN  TOE  CADMIUM- TELLURIUM-MERCURY  SYSTEM 
(All  units  in  Joules  per  gran-atom  (mole  of  atoms),  °K) 


L-liquid  phase,  V=vapor  phase,  S=tinc-blendc  phase,  A=fcc 
R=rhombohedral  (stable  Hg) ,  T’trigonal  (stable  Te) ,  E=hcp  (stable  Cd) 
x=atom  fraction  Te,  y=atom  fraction  Hg,  P=pressure  (atmospheres) 
R-8.314  J/g.at'K 


Lattice  Stability  Values 


GL.  -  G?  = 
Hg  Hg 

pA  R 

TOg  '  TOg  " 


C-  -  Sg 


4 
4  -  4  1 
4  -  4  ■ 
K2-  4 


2292-9. 832T;  T=234CK 
2297-1.464T 
59174-93. 303T+RT£nP 
17489-24. 184T;  T=  723°K 
17489-15. 816T 
58409-45 .982T+0 . SRT£nP 


GCd  ‘  GCd  '  6J92-10.418T;  T=594K 

C*.  -  g!!  ,  =  1924-2. 050T 
Cd  Cd 

GCd  *  GCd  =  ,520-98 ,6S9T+RT£nP 


Ternary  Phases 

GL  »  (l**-y)Gpd  ♦  xGGe  ♦  yG^g+  RT  [ (l-x-y)£n(l-x-y)+x£nx+y£ny] 

-xyU+y)’1  [(30962+5. 02T)x+(3514+16.74T)y]  ♦  y(l-x-y) [-10878+4 . 184T] 
♦x(I-x-y) (1-y)*1 [ (1-x-y) ( 112968 -159T)  -x(75312+20.92T)] -  33472  xy(l-x-y) 


GV  »  (l-X-y)GJd  +  0.5x  G^  +  y  <{g 

♦  RT((l-x-y)£n(l-x-y)+x£nx+y£ny) 


GS  *  (0.5-y)  G*d  +  0.5  G^e  +  y  G^g  ♦  (l-2y) [-61191 +13. 075T] +2y[-26798+7 . 11T] 

♦  RT[yfny  ♦  (0.5-y)£n(0.5-y)-0.s£n0.5]  -4180y(l-2y) 

0  1  y  <  0.5 


TABLE  8 

COMPARISON  OF  CALCULATED  AND  OBSERVED  (16) FREE  ENERGY  OF  FORMATION  AT  298*K 


Compound 

AH  [2981 

.  (J./g.at.) 

AS  [2981  (J./g.at.0 

Cd0.5Te0.5 

-50900  i 

400  (obs.) 

-4.18  t  1.25  (obs.) 

-51500 

(calc.) 

-4.14  (calc.) 

Hg0.STe0.5 

-15900  ♦ 

2100  (obs) 

-6.28  ±  1.26  (obs) 

-18050 

(calc.) 

-3.39  (calc.) 

Figure  26.  Calculated  Partial  Te-Cd  Phase  Diagram  at  1  Atmosphere, 


Hg 

Figure  28.  Calculated  Partial  Hg-Cd  Phase  Diagram  at  1  Atmosphere. 


Figure  29.  Calculated  Partial  Hg-Cd  Phase  Diagram  at  16  Atmospheres 


Molt  Fraction  x 


HfTt 


•  -  Steininger  (39) 

o  LiqwOut  Oau  -  Sian  and  Mtwnham  (40) 

V  Sol»duaOau  -  Ha'<nan  and  Sarauit  (41) 

□  SoMuxOtta  -  Schm.t and  SpttricNitidtr  (42) 

Figure  32.  Liquidus  and  Solidus  Curves  for  Pseudo- 
binary  Hg^^Cd^Te  Compositions  (39) 


Figure  33.  Composite  P-T-x  Diagram  for  the  Pseudo- 
binary  HgTe-CdTe  System  (39) 


Figure  34.  Calculated  Isothermal  Sections  in  the  Cd-Te-llg  System 


V  CALCULATION  OF  TERNARY  FLUORIDE  GLASS  COMPOSITIONS 


L.  Kaufman,  J.  Agren,  J.  Nell  and  F.  Hayes,  CALPHAD  (1983)  ]_  71~ 


(This  paper  was  preaented  at  CALPHAD  XI,  Argonne,  Illinois  May  1982 

ABSTRACT.  The  recent  discovery  of  a  new  family  of  non-oxide  glasses  based  on 

mixtures  of  ZrF^  or  HfF^  with  other  metallic  fluorides  by  M.  Poulain 

and  coworkers  offers  great  potential  in  optical  fiber,  window  and 
sour ce /de t ec t or  application.  Due  to  the  limited  phase  diagram  data 
available  for  the  binary,  ternary  and  multicomponent  fluoride 
systems  currently  employed  to  synthesize  these  glasses  most  of  the 
progress  in  identifying  new  compositions  has  proceeded  along 
empirical  lines.  In  order  to  remedy  this  situation,  the  CALPHAD 
method  for  coupling  phase  diagram  and  thermochemical  data  has  been 
applied  to  develop  a  data  base  covering  metallic  fluorides.  The 
objective  is  to  permit  computation  of  multicomponent  phase  diagrams 
which  can  be  used  to  identify  the  composition  range  where  the  liquid 
is  most  stable.  The  latter  offers  opportunities  for  glass  formation 
a a  demonstrated  by  predictions  of  new  metallic  glasses.  Currently 
the  data  base  covers  combinations  of  0.2  ZrF^(ZF),  0.25  LaF  (LF), 

0.333  BaF2(BF),  0.333  PbFjCPF),  0.5  NaF(NF),  0.5RbF(RF),  0.5CsF(CF)  and 

0.5  KF(KF)  which  have  been  developed  along  the  lines  described 
earlier  fox  Itl-V,  II-VI  snd  SIAL0N  systems.  These  results  were 
used  to  generate  computed  liquidus  contours  in  LF-ZF-BF  and  BF-ZF-NF 
to  disclose  the  range  of  composition  in  which  the  liquid  has  the 
greatest  stability.  These  compositions  agree  well  with  those  in 
which  Poulain  and  coworkers  have  discovered  glass  formation. 

I  .  Introduct  ion 

The  utility  of  computer  based  methods  for  coupling  phase  diagrams  and 
thermochemical  data  for  metallic  systems  has  been  well  documented  in  many 
papers  published  in  this  journal.  A  considerable  effort  is  being  applied 
toward  developing  an  extensive  bese  for  metallic  systems.  Recently  similar 
efforts  have  begun  in  order  to  provide  a  similar  facility  with  ceramic  systems 
(I)  and  semiconducting  I1I-V  compound  systems  (2).  The  prediction  of 
compositions  of  high  liquid  stability  in  the  t i t an ium-bc r y 1 1 ium-z i r con ium 
in  1973  (3)  was  followed  by  development  of  commercial  metallic  glasses  in  1975 
at  Allied  Chemical  (A).  In  view  of  the  recent  discoveries  by  M.  Poulain  and 
coworkers  (5,6)  that  zirconium  fluoride  glasses  offer  potential  as  optical 
materials  it  seems  appropriate  to  consider  development  of  a  data  base  capable 
of  computing  composition  regions  of  liquid  stability  where  glass  formation 
could  be  facilitated. 


2  ■  Description  of  the  Thermochemical  System 
Employed  to  Characterize  Solution  and  Compound  Phases 

The  method  utilized  for  describing  solution  and  compound  phases  is  the  same 
as  that  employed  earlier  (1,2)  comprising  some  symbolic  usage  which  facilitates 
data  handling  as  indicated  below.  The  free  energy,  G  ,  of  a  liquid  (solution) 
phase,  L,  in  the  binary  system  ZrF  -RbF  is  given  by  Equation  (1)  where  T  is  in 
Kelvins,  R-8.314  J/g.  at  °K  ,  while  x  is  the  atomic  fraction  of  RF  (  i  .  e  .  RF  =  -j  Rb  F  ) 
and  (l-x)  is  the  fraction  of  ZF  (i.e.  l/5ZrF,).  The  mass  basis  is  thus  one  mole 

of  atoms  (i.e.  a  gram-atom) 

CL-( l-x)°c£  *x°g!t  +RT(xlnx+(  I  -  x ) In (  l-x) )*x(  l-x)[LZFRF(  l-x)+xLRFZF  ]  J/g  .  at  (I) 


where  °G^_  and  °CpE.  are  the  free  energies  of  one  gran  atom  of  pure  liquid  ZF  and 
Z  F  R  F 

RF .  Table  t  defines  the  lattice  stabilities  of  the  liquid  and  solid  forms  of 
ZF,  RF  and  the  other  components  of  current  interest.  These  data  derived  from 
various  compilations  of  t he rmochemic a  1  and  phase  diagram  data,  (7—15),  when 
combined  with  the  solution  and  compound  phase  parameters  shown  in  Tables  2  and  3 
permit  calculation  of  the  binary  phase  diagrams  shown  in  Figures  1—10.  The 
solution  par ame t er sLZFRF  and  LRFZF,  which  describe  the  liquid  ZrF.-RbF  solution 
are  listed  in  Table  2  (  LZFRF.-707  I  0  +  29 . 2  88T  J/g.at,  LRFZF=-3  1  380  +  29 . 288T 
Joules / g .  at )  .  Similar  parameters  for  the  liquid  and  solid  phases  are  listed  in 
Table  2.  The  free  energy  of  the  solid  phases  are  described  in  a  manner  similar 
to  equation  (1).  Thus,  for  example,  the  free  energy  of  the  soldium  chloride 
type  solid  S  phase  in  the  ZF-RF  system  is  given  by 


CS-(  l-x)0G2F+x°C®f.  +  RT(xtnx  +  (  l-x)8.n<  1 -x)  ) +x(  1-x  )  [  SZFRF  (  1 -x  ) +  XSKFZF  ]  J/g.  at.  (21 


The  solution  parameters  for  the  S  phase  ( SZFRF=SRFZF=  1 6  7 36  J/g.at)  are  listed  in 

Table  2.  The  free  energy  difference  between  the  L  and  S  forms  of  ZF  and  RF, 

i.e.  ZFZFLS“°G,_-°G,  and  °G"'  -°G„„=RFRFLS  are  listed  in  Table  I.  The  free 
Ct  L  r  Kr  Hr 

energy  of  a  compound  phase  such  as  X  =  (  I / 7  )  ( Zr F^ • RbF )  =  ZFQ  7|4R1'o  286  *s  'JclF»eJ 

on  the  basis  of  Table  3  in  terms  of  the  base  phase,  the  compound  parameter,  C, 
and  the  stoichiometry  by  equation  (3) 

GX-0.7I4°g“f+0.286°g“f+(0.7I4)(0.286)[0.7I4LZFRF+0.286  LRFZF-C]  J/g.at.  (3) 

Reference  to  Table  2  shows  that  LZFRF«- 707  1 0  +  29 . 288T  J/g.at  and  LRFZF  = 

-31380  +  29. 288T  J/g.at.  while  Table  3  shows  that  C--20669+4  1  .84T  J/g.at.  Thus 

CX-0.7 l4°cj?F+0.286°c”F-79 I2-2.56T  J/g.at.  (4) 

At  298K  reference  to  Figure  I  shows  that  the  stable  form  of  RbF  is  Lhe 
sodium  chloride  S  form.  Table  I  shows  that 

°g“f-°G®f-RFRFLS-RFRFLM-i2866-1.381T  J/g.at.  (5) 

Substitution  of  equation  (5)  into  equation  (4)  yields  an  expression  for  the  free 
energy  of  formation  of  (  I / 7 )  ( Zr F,  •  RbF )  at  298K  from  the  stable  forms  of  its 
compound  components  as, 

AGf[298K]-GX-0.7  14  g“f-0.286  G*f— 4234  +  2 .95T  J/g.at.  (6) 

The  latter  value  is  shown  in  Table  3. 

Experimental  information  concerning  the  thermochemical  properties  of 
fluoride  solutions  and  mixed  f  luoride  compounds  is  scarce.  The  only  experimental 
data  located  (IS)  is  shown  in  Table  3  where  experimental  measurements  of  the 

heat  of  formation  of  liquid  ZrF^-NaF  solutions  at  I303K  are  compared  with 

calculated  values  on  the  basis  of  Equation  (I)  and  Table  2.  The  table  lists 
the  measured  molar  heat  of  formation  (15)  for  specific  molar  concentrations 
along  with  the  calculated  values.  The  latter  are  obtained  by  converting  the 
molar  composition  into  appropriate  fractions  of  NF  (i.e.  l/2NaF)  and  ZF 
(i.e.  (l/5)ZrF^).  Thus  0.75NaF-0.25  ZrF^  is  equal  to  I.5NF-I.25ZF  which  is  in 

turn  equal  to  2.75  gram  atoms  of  ZFQ  ^j^NF^  545-  The  calculated  heat  of 

mixing  of  the  liquid  at  'this  composition  from  Equation  (1)  and  Table  2  at  a 
composition  xN  ■0.545  is  -10090  Joules/g.at.  Hence  the  calculated  heat  of 
formation  of  a  0.75NaF-0.25  ZrF,  molar  liquid  is  2 . 75x(-  10090  )— 2  7 760 
Joules/mole.  The  calculated  values  in  Table  4  are  in  agreement  with  the 
measured  results. 


& 


Calculation  of  Quasi-Ternary  Phase  Diagrams 


The  free  energy  of  ternary  solution  nhases  was  synthesized  from  the 
binary  solution  phases  on  the  basis  of  Kl  .ler’s  equation  as  in  the  previous 
papers  (1,2).  On  this  basis  the  free  energy  of  the  liquid  phase  in  the 
tF-ZF-BF  system  in  Figure  II  where  x  is  the  atom  fraction  ZF,  y  is  the  atom 
fraction  BF  and  1-x-y  is  the  atom  fraction  LF  is  given  by  Equation  (7) 


o 
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TABLE  I 

SUMMARY  OF  LATTICE  STABILITY  PARAMETERS 
(All  Units  in  Joules  per  gram  atom  (mole  of  atoms),  T  in  Kelvins) 

ZF  =  (  l/5)ZrF4,  RF-(l/2)RbF,  NF»(l/2)NaF,  CF«(l/2)CsF,  KF=(1/2)KF,  BK  =  (  I / 3  )  Ba F 2 

LF-(  l/4)L«F3,  PF-(t/3)PbF2 

S-Sodium  Chloride  S t rue t ure,  S t ab le  form  of  NF,  KF ,  RF  and  CF 

L-Liquid,  M-Stable  form  of  ZF,  K-Stable  form  of  BF,  J=Stable  form  of  LF , 

C-Stable  form  of  PF 


NFNFLS  =  (  l / 2 ) NaF (Liquid )-(  1  / 2  ) NaF ( S  ,  sodium  chloride) 


ZFZFLM-I2845-I0.669T 
ZFZFLS-  0-11. 966T 
ZFZFLK-  0-  6.069T 
ZFZFLJ-  0-  7  .  1  13T 
ZFZFLC-  0-  5.3I4T 

RFRFLS-12866-I2.050T 
RFRFLM-  0-I0.669T 

NFNFLS- 16799- 13. 263T 
NFNFLM-  0-10. 669T 
NFNFLK-  0-  6.067T 

CFCFLS- 10857- 1 1 . I29T 
CFCFLM-  0-10. 669T 


KFKFLS- 14121-12 .463T 
KFKFLM-  0-10.669T 

BFBFLK-  9481-  6.067T 
BFBFLM-  0-10. 669T 
BFBFLS-  0-13. 263T 
BFBFLJ-  8786-  7 . I 13T 

LFLFLJ- 12560-  7 . I 13T 
LFLFLM-  0-10.669T 
LFLFLK-  8673-  6.069T 

PFPFLC-  5803-  5.314T 
PFPFLM-  0-10. 669T 


•These  differences  specify  the  free  energy  of  one  phase(i.e.  liquid)  minus  the 
free  energy  of  the  second  phase  (i.e.  sodium  chloride)  for  a  given  compound. 


TABLE  2 

SUMMARY  OF  SOLUTION  PHASE  PARAMETERS 
(All  Units  in  Joules  per  gram  atom  (mole  of  atoms),  T  in  Kelvins) 


LZFRF  — 707  10  +  29. 288T, 
LRFZF--3 I  380  +  29. 288T, 
MZFRF-MRFZF- 16736  , 
SZFRF-SRFZF- 16736  , 
LZFNF  — 54392  +  29. 288T, 
LNFZF--29288+29.288T, 
MZFNF-MNFZF- 16736  , 
SZFNF-SNFZF- 16736  , 


LZFCF--75312+29.288T, 
LCFZF--33472+29.288T, 
MZFCF-MCFZF- 16736, 
SZFCF-SCFZF-  16736  , 
LZFKF--753I2+29.288T, 
LKFZF--33472+29.288T, 
MZFKF-MKFZF-  16736  , 
SZFKK-SKFZF-  16736  , 


LZFBF--  10154  .LPFZF-LZFPF- I  2552 
LBFZF-  18272  ,  CPFZF-CZFPF-4  I  84 0 
MZFBF-MBFZF -4  I  84 0 , MPFZF-MZF PF-4 1840 
KZFBF-KBFZF-4 1840, 

LZFLF— 11088  .LRFNF-  12552 
LLFZF-  22595  ,LNFRF- 1 t . OT 
MZFLF-MLFZF  =  4  I  84 0  ,  SRFNF- 3 7 6 5 6 
JZFLF-JLFZF-4 I  84 0 , SN FRF-4 05 84 


LNFBF-  12552,  LBFNF-16736,  SNFBF- SBFNF-KNFBF-KBFNF-4 I  84 0 

LBFLF-LLFBF-2  1757-6 .82T ,  JBFLF- JLFBF-33  1  3 7- I  2 . 552T ,  KBFLF -KLFBF -2  5104-12. 552T 


"w  .'•  ."*  «> 


'»  .  ‘  s‘" 


•>av  v'v:%vv,'j>r*:vNJ 
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TABLE  3 

SUMMARY  OF  COMPOUND  PARAMETERS  FOR  BINARY  SYSTEMS 
(All  units  in  Joules  per  gras  atom  (mole  of  atoms),  T  in  Kelvins) 


■ 

ft 

1 

Compound  Name 

S  toich iome  try 

Base 

Compound  Parameter  Free  Energy  of 
Formation  (from 
stable  components) 

C 

(  Joules  /  g  •  at  ) 

ag£ 

(Joules/e.  at) 

.1 

• 

(  1 / 7 ) ( ZrF^ • RbF ) 

X 

ZF0.7I4RF0.286 

M 

-20669+4 1 . 84  T 

-4234-2. 95T 

j 

( l/30)(4ZrF, • 5RbF ) 

4 

R 

ZF0.667RF0.333 

M 

-19644+45. 52T 

-4  142-4. 07T 

j* 

(  1  /9 ) ( ZrF, *2RbF) 

4 

Q 

ZF0 . 556RF0 .444 

M 

-  15297  +  5 1 . 09T 

-366  1-6 .00T 

(  1/11) (ZrF, • 3RbF ) 

P 

ZF0 . 4  55RF0 . 54  5 

S 

-  13389  +  58 .58T 

-3054-6. 67T 

(  l/7)(ZrF4-NaF) 

X 

ZF0.7 14NF0.286 

M 

-  17575  +  4  1 .2  IT 

-  125  1-3  .  18T 

1 

(  1/1  l)(2rF,  ONaF) 

p 

ZF0.455NF0.545 

S 

-  1  1456  +  5  1 .09T 

-  1406-4  .82T 

p 

(  1  / 7  )  ( ZrF .  •  CsF ) 

4 

X 

ZF0 . 7 14CF0 . 286 

M 

-17322  +  4  1  . 84  T 

-6293-2. 68T 

j 

(  l/9)(ZrF4  -2CsF) 

q 

ZF0.556CF0.444 

M 

-  16928  +  5  1  .97T 

-5012-5. 79T 

J 

i 

(1/1  1  )(ZrF4  -3CSF) 

p 

ZF0.455CF0.545 

S 

-12259+58. 58T 

-4  I34-6.67T 

(  I/7)(ZrF4  ‘KF) 

X 

ZF0.7 14KF0.286 

M 

-  19874  +  4  ! . 84T 

-4828-3. 08T 

• 

(  1  / 9 ) ( ZrF  • 2KF ) 

4 

Q 

ZF0.556KF0.444 

M 

-13974+49. 62T 

-4289-5. 82T 

(1/1  1  )(ZrF4  -3KF) 

p 

ZF0.455KF0.545 

S 

-14351+58. 58T 

-36  19-6. 67T 

c 


TABLE  A 

COMPARISON  OF  CALCULATED  AND  OBSERVED(IS)  HEATS  OF 
FORMATION  OF  ZrF^-NaF  LIQUID  SOLUTIONS  (T-I303°K) 


Molar  Composition 


0.95NaF-0.05ZrF. 
0.90NiF-0.  lOZrF, 
0.85NaF-0.  l5ZrF4 
0.80NaF-0.20ZrF4 
0.75NaF-Q.25ZrF4 


AH  £ (Observed ) ( 15 ) 
Joule  a/mol 

AH£(Calculated)* 

Joules/mole 

XNF 

-6485*2510 

-7095 

0.884 

-12135*2510 

-13570 

0.783 

-17990*2510 

-19230 

0.694 

-24895*2510 

-23970 

0.615 

-29290*25  10 

-27780 

0.545 

XNF  gram  atoms 

mole _ 


EHL  ( Joules/g  .at . ) 
Heat  of  Mixing  of  Liquid 


Q 


0.884  2.15  ZFq  i16NF0  "I .9  NF-0.25ZF  -3300 

0.783  2.30  zr“''  Ilir  -1.8  NF-0.50ZF  -5900 

0.694  2.45  ZFq  3Q6NFo'fiqi*l  NF-0.75ZF  -7850 

0.615  2.60  ZFn  ,BSNFns  it*1  -6  NF-I.00ZF  -9220 

0.545  2.75  ZFo’4S5NF0  545"1  ,5  NF-I-25ZF  -10100 


*AH£«(number  of  gram  atoms/mole)  times 


ehl 


© 


C 
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Figure  4.  Calculated  Zirconium  Fluoride  -  Potassium  Fluoride 

Phase  Diagram 
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CL-( l-x-y)0G^F**°G2F  +  y°CgF  +  RT[ ( I  -x-y )  in{  l-x-y )  *x£.nx  +  y  iny  ] 

♦  (  l-x-y )x(  1  -y  )  “  '  [ <  l-x-y)LLFZF+xLZFLF]*xy(x+y)“l[xLZFBF+yLBFZF] 

♦  ( l-x-y)y(  l-x)- (  I-x-y)LLFBFeyLBFLF]  J/g.at  (7) 

All  of  Che  solution  parameters  required  to  specify  Equation  (7)  are  given  in 
Table  2.  The  remaining  solution  phases  in  the  ternary  system  i.e.  S,M,J,K  etc. 
are  defined  in  similar  fashion.  Such  definition  permits  calculation  of  the 
LF-ZF-BF  ternary  shown  in  Figure  11.  The  liquidus  boundaries  calculated 
between  1473K  and  773K  define  the  compositions  of  maximum  liquid  stability. 
Reference  to  Figure  11  shows  that  these  results  agree  well  with  the  glass 
forming  range  established  by  Lecoq  and  Poulain  (5). 

In  order  to  consider  a  case  like  BF-ZF-NF  in  Figure  12  where  compound 
phases  (like  X  and  P)  are  present  the  usual  definition  (I)  is  employed  in 
order  to  define  the  free  energy  of  the  compound  phase  NFq  j^jZFq  in  the 

ternary  system  BF-ZF-NF  where  x  is  the  fraction  of  ZF,  y  is  the  fraction  of 
BF  and  s  is  the  fraction  of  NF.  In  this  case  the  free  energy  of  P  is 

GP  -  *oG®F+x0G|F*y°G*F+(l-(y/(l-x)))AGA^(y/(l-x))AGB 

+RT(yZny+sin*-( 1 -x)in ( l-x ) )  J/g.at.  (8) 

Comparison  of  Equations  (8)  and  (3)  show  that  with  the  compound  running  from 
the  stable  NF0  545ZF0 . 4 55  t0  the  unatable  BF0 . 545ZF0 .455 ’  x“0,455  and  *  ia 
equal  to  0.545-y  where  y  runs  from  2ero  to  0.545.  Thus 

AGa-(0.545)(0.4  55)[0.545LNFZF+0.4  55LZFNF+1  14  56-5  1 .09T]  J/g.at  (9) 

The  free  energy  of  the  unstable  compound  BFQ  545^0  455  was  defined  by 
Equation  10  as 

AC  -(0.545)(0.455)[0.545LBFZF+0.455LZFBF-30543]  J/g.at.  (10) 

a 


in  order  to  complete  the  calculation  of  the  liquidus  contours  to  define  the 
tange  of  liquid  stability  in  Figure  12. 

These  results  illustrate  how  the  development  of  a  data  base  for  fluoride 
systems  can  be  employed  to  usefully  predict  the  most  appropriate  compositions 
to  investigate  in  the  search  for  tirconium  fluoride  base  glasses.  Needless  to 
say  the  method  has  even  greater  value  for  multi-component  systems  containing 
more  than  three  compounds. 

4  .  References 

1.  L.  Kaufman,  F.  Hayes  and  D.  Birnie  CALPHAD  5^  163  (1981). 

2.  L.  Kaufman,  J.  Nell,  K.  Taylor  and  F.  Hayes  CALPHAD  £  185  (1981). 

j  .  L.  Kaufman  and  H.  Nesor,  Titanium  Science  and  Technology  R.I.  Jaffee  and 
H.  Burte  Editors,  Plenum  Press,  New  York  2  773  (  1973)  . 

4.  L.  Tanner  and  R.  Ray,  Scientific  American  233  15  (1975). 

5.  A.  Lecoq  and  K.  Poulain,  J.  Non-Crystalline  Solids  34.  101  (1979). 

6.  K.  Poulain,  M.  Poulain,  and  J.  Lucas,  Rev.  de  Chem.  Mineral  J_6  267  (1979). 

7.  J.  Lumsden,  Thermodynamics  of  Molten  Salt  Mixtures  Academic  Press,  London 
(1966). 

8.  0.  Kubaschewski  and  C.B.  Alcock,  Metallurgical  Thermochemistry,  Pergamon 

Press,  Oxford  (Fifth  Edit  ion) (  19  79).  ~~ 

9.  JANAF  Thermochemical  Tables,  Dow  Chemical  Co.,  Midland,  Michigan. 

10.  R.E.  Thoma,  Advances  in  Molten  Salt  Chemistry.  J.  Braunstein,  G.  Manatov 
and  G.P.  Smith ,  eds.  Plenum  Press,  N.Y.  (1975)  p.  302. 

11.  E.  Levin,  C.R.  Robbins,  and  H.  McMurdie,  Phase  Diagrams  for  Ceramists, 

American  Ceramics  Society,  Columbus,  Ohio  (  196  4).  — ™~ 

12.  Ibid.,  First  Supplement  (1969). 

13.  E.  Levin  and  H.  McMurdie,  Second  Supplement  (1975). 

14.  R.S.  Roth,  T.  Negas  and  L.P.  Cook,  Ibid.,  Third  Supplement  (1981). 

15.  A.  Fontana  and  K.  Winand,  J.  Nuclear  Materials  44  305  (1972). 


VI  CALCULATION  OF  ORDERING  TEMPERATURES  FOR  THE 
QUASI-BINARY  SECTION  Fe3Al-Fe3Si 


1.  Introductioh 


Where  limited  experimental  data  is  available,  the  method  adopted  by  Inden  (1-7) 
is  attractive  because  it  requires  only  two  ordering  parameters  for  each  of  the  three  component 
binary  systems,  provided  the  temperatures  of  interest  are  above  any  magnetic  transition 
temperatures.  Those  ordering  parameters  are  the  first  and  second  nearest  neighbor  (chemical) 

(1) 

interaction  parameters  denoted  W  and  Wv  .  These  can  be  related  to  heats  of  formation  as 
well  as  critical  ordering  temperatures,  as  will  be  shown  in  the  following  section.  The  Fe-Si 
binary  system  (Fig.  1)  has  already  been  characterized  in  some  detail  (3,4)  and  the  previously 

determined  values  of  W^_  and  W^_  have  been  used  unaltered  in  the  present  work. 

re-ui  Fe-si 

Table  I).  Their  validity  has  been  confirmed  by  the  successful  prediction  of  ordering 
equilibria  in  ternary  Fe-Si-Co  alloys  (8).  The  Fe-Al  system  is  inherently  more  connlicated 
(9,10)  (Fig.  2)  and  several  additional  effects  have  been  recently  added  by  KBster  et  al  (111. 
Nonetheless  the  major  features  have  been  reproduced  by  the  simplified  Inden  model(12;.  Only 

a  minor  change  has  been  made  in  the  value  of  ^  used  previously  by  Inden  (6)  in  the 

light  of  the  available  experimental  data  for  Fe-Si-Al  alloys  (Table  II).  The  t.jn  vation  of 

suitable  values  of  suitable  antl  W^si-\1  v:llues  Poses  3  problem  in  so  far  that  all 

experimentally  accessible  properties  for  these  alloys  refer  to  either  FCC  Aluminum  or  diamond 

cubic  Silicon.  One  has  therefore  to  extract  the  relevant  values  of  wU)  and 

7.  bi-.-u 

W1"'..  ..  from  whatever  data  is  available  for  BCC  ternary  allovs. 

Si-Al 

2.  Experimental  Information  on  Fe-Al-Si  Alloys 

Lihl,  Burger,  Sturm  and  Ebel  (14)  have  given  a  general  indication  of  the  location 
of  the  ordering  region  in  Fe-Al-Si  alloys  (Fig.  3).  However  there  is  considerable  divergence 
between  this  estimate  and  the  investigation  of  Katsnel'son  and  Polishchuk  (13).  More  weight 
has  been  given  to  the  latter  work  as  it  utilized  relative  intensities  of  superlattice  lines 
(Fig.  4)  whereas  the  work  by  Ebel  et  al  merely  interpreted  changes  in  d-spacing  as  ordering 
effects.  An  NMR  study  of  Fe3Si^  ^Al^  alloys  (15)  suggests  that  Aluminum  and  Silicon 

substitute  interchangeably  on  one  specific  sublattice,  which  confirms  that  the  interaction 
energy  between  Silicon  and  Aluminum  is  relatively  low. 


Fig,  4  Critical  Temperatures  suggested  by  the  temperature 
dependence  of  superstructure  lines  from  Katsnel'son  (13) 
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Ordering  Model 


(a)  Definition  of  order  parameters 

The  main  feature  of  the  3h'G  model  is  tne  description  of  atomic  distributions  inside 
specifically  defined  sublattices.  Thechoice  follows  from  the  kind  of  interaction  between 
atoms.  Any  BCC  lattice  may  be  divided  into  four  FCC  sublattices,  I,  II,  III,  and  IV,  with 

the  atomic  configurations  described  by  occupation  probabilities  'p.1'1  of  sublattice  site  ' L* 
by  the  .component  * i * .  (Fig.  S)  1 


Silas  of  sub-lattice  I  — o 


N 


O 


Fig.  5  Relative  position  of  sublattices  in  the  DO.  structure. 

The  atomic  configuration  in  a  binary  alloy  can  then  be  described  with  the  aid  of  the 

following  three  parameters,  which  are  a  combination  of  the  independent  occupation  probabiliti 


x 


y 


z 


i 

1 

» 

* 


r  .  n  rr  n  rrr  _  rv, 

Cpa  Pa  "  pA  '  pA  5 

(la) 

CPA”1  *  pAIV> 

c r 

CPA'  -  PA") 

(lc) 

The  greatest  value  of  these  parameters  corresponds  to  the  highest  degree  of  order,  and  the 
concentrations  of  A  and  B,  (C.  and  Cg)  set  limits  on  the  values  of  x,  y  and  z.  For  the 
binary  case  these  are  given  by: 


o 


,X. 


:  ».*  %•  v 

■ 


-V.-w*. 


*>■  *w 


-?4- 


Ca) 

0S<  x  <  CB  for  0  <  CB  <  0.5 

(2a) 

(M 

0  s  y  f  (CA  -  x,  CB  ♦  x) 

(2b) 

(c) 

0  j  z  <  min  (CA  ♦  x,  Cg  -  x) 

(2c) 

Case  1 

The  random  BCC  lattice  (A2)  has  an  equal  distribution  of  the  atoms  on  the  four  sublattices 
and  corresponds  to  x  «  y  =  z  =  0. 

Case  2 

The  state  x  j*  0,  y  ®  z  «  0  indicates  'A'  preferentially  occupying  sublattices  I  and  II  (x  >  0), 
or  III  and  IV  (x  0).  It  follows  that  the  number  of  AB  neighbors  in  nearest  neighbor  (n.n) 
positions  is  increased,  i.e.  this  corresponds  to  the  B2  structure. 

Case  3 

It  follows  that  y  ^  0,  describes  a  surplus  of  A  atoms  in  next  nearest  neighbor  (n.n.n) 
positions,  which  corresponds  to  the  DO  structure.  As  far  as  the  present  problem  is  con¬ 
cerned  we  need  not  consider  z.  (See  Refs.  1,2) 

(b)  Free  Energy  Expressions  and  Interaction  Parameters 

The  most  stable  atomic  configuration  then  follows  from  the  minimum  value  of  the 

BCC 

configurational  free  energy/mole  AG  which  is  given  by: 


AGBCC 

BMARY 


•  U°  -  NCxCb  (4w(l)  ♦  3w(2))  -  5 

♦  3W(2)  (y2  ♦  z2) 


(8  WO)  -  6W(2))  x2 


NkT  Z 


L,  L  '  L,  L 
PA  ln  P.A  *  pB  ln  pB 


(3) 


U°  is  the  internal  energy  of  pure  components;  the  second  term  is  the  energy  of  mixing;  the 
third  term  gives  the  energy  contribution  due  to  ordering,  and  the  fourth  term  is  the  entropy 
contribution  of  the  configuration  which  is  also  expressed  in  terms  of  x,  y,  s.  V  is  Avagadro1 

m  o\ 


number,  and  k  the  Boltzmann  constant.  VP  and  VP2'  are  defined  below  and  for  eqn  (3)  to  be 

-24  -24 

valid,  these  must  be  expressed  in  k  units  where  1  k  unit  •  3.3  x  10  cal  or  13.3  x  10  J. 
The  values  of  and  are  related  to  the  nearest  and  next  nearest  neighbor  bond  energies 

v  <k>  as  follows: 


wO)  -  -  zifoW  ♦  vMC>)  ♦  VggCO 

(4a) 

m  -  -  afoW  ♦  vM(2)  ♦  vbb(2) 

a* 

According  to  this  definition,  positive  values  for  and  W  ^  correspond  to  greater 

affinity  between  unlike  atoms.  All  interchange  energies  are  presumed  to  be  independent  of 
temperature  and  of  environment  of  other  atom  sites  where  the  interchange  process  is  considered 
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i1.  .ui.uii U|  I LIJ II.  »i 


PJ  •_*  >1  ii  j  ■  II  »I  n  ■! 


■  w  ?■„'  ,rj  r 


in  ■.  ji  u  j  i .  n 


fn  (2)  BCC 

If  the  values  of  Hr  and  Hr  are  known,  AG  can  be  calculated  with  respect  to  the  ordering 
parameters  x,  y,  2  using  the  necessary  equilibrium  conditions  including  particularly: 


Of 


m 


1# 


K:;, 


3G  3G  3G 

—  ■  ■  —  *  --  a  _ 

3x  3y  a: 

and  other  subsidiary  conditions. 


(S) 


(c)  Expressions  for  critical  ordering  temperatures  in  binary  alloys 

The  critical  temperatures  Tx(A2/B2),  T^(B2/DO.),  define  the  regions  in  which 

minimum  free  energy  is  obtained  by  the  A2,  B2  and  DO.  configurations  respectively.  For  the 
binary  case  these  are  given  by: 


kT  •  ^2/B2  «  (SWCO  -  6IV(2})  .  C,C„ 
x  (BWG)  BWG  BSVG  A  “ 


(6a) 


j^j-82/ED  3 

(BWG) 


6W(2)  •  (C, 
BWG  A 


X)  (CB  *  X) 


(6b) 


.here  x  vs  the  appropriate  value  of  the  ordering  coefficient  at  the  temperature  Ty. 

Since  are  energy  parameters,  they  are  most  easily  determined  from  energy  measurements 

like  enthalpies  of  mixing  of  random  alloys,  or  enthalpies  of  formation  AH  (C,T)  of  alloys 
with  given  atomic  configurations.  These  entities  can  immediately  be  expressed  by  means  of 

for  a  binary  alloy  Aj  ^.B^  with  reference  to  the  pure  components  in  the  same  crystal 

structure  as  the  alloy  formed.  For  BCC  alloys: 


AHA2 

•  -  N  C(1  -  C)  £  4W(»)  ♦  5W(2) J 

....  (7) 

ah82 

•  -  N  £c(4»0)  ♦  3i(2h-  6C2  W(2)J 

for  C  <  O.S 

ah003 

-  -  M  C  (4WC»)  ♦  31V(z)) 

for  0  <  C  <  0.25 

OR 

-  -  N  £l.SW(2)  ♦  C (4W C 1 )  -  3W(2))  J 

for  0.25  <  C  <  0.5 

AH  832 

-  -  N  £  C  (4W(»)  ♦  3W(2))  -  4C2  !t'(0  J 

for  C  <  0.5 

1  h>  values  obtained  from  an  analysis  of  critical  temperatures  (e.g.  via  equations  6a,  ob), 

QtoG  /  IQ 

can  be  related  to  the  Wv  1  values  in  thermochemical  equations  (e.g.  equ.  7V  by  a  scale  factor 
x  (listed  in  Table  I).  See  References  (152).  (The  differences  between  and 

are  essentially  due  to  the  omission  of  short  range  order  in  the  BWG  model.) 


o 
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4.  Expressions  for  the  Free  Energy  and  critical  ordering  temperature  in 
ternary  alloys  as  used  in  the  present  work 

In  the  ternary  case  there  is  a  slight  difference  in  the  definition  of  the  order 
parameters  (S).  If  the  total  number  of  atoms  is  N,  distributed  over  4N  lattice  sites,  the 

occupation  probabilities  p^  (i  *  A,BfC,;  L  *  I,  II,  III,  IV)  are  combined  to  give  the  order 

parameters  as  follows: 

(a)  x.  -  J  (p.1  ♦  Pj11  -  Pj111  -  PjIV)  C8a) 

(b)  /j  -  1  (Pj111  -  PjIV)  (8b) 

(c)  Zj  •  1  (Pj1  "  Pj111)  (8c) 

with  j  ■  A,  B 

A  similar  condition  to  Equation  (5)  occurs  except  now  there  are  two  x  parameters,  two  y 
parameters  and  two  z  parameters,  corresponding  to  the  degree  of  order  of  two  solutes  and  not 
just  one.  For  the  purpose  of  determining  T  and  T  temperatures  it  is  only  necessary  to 

x  y 

consider  and  Xg.  The  limiting  values  of  x,,Xg  in  the  ternary  case  are: 

fa)  The  maximum  value  of  x^  is  the  smaller  of  a  or  (1-a).  (9a) 

(b)  The  maximum  value  of  xg  is  the  smaller  of  b  or  (c  -  x^)  (9b) 


In  the  specific  case  ot  Fe^Al  ^^Si^  alloys,  where  a  =  Fe,  b  =  A1  and  c  =  Si,  this  leac  > 
to  a  situation  where  (a)  is  always  0.75  and  therefore  xA(raax)  =  0-25.  Likewise  since 

3  -b.  The  configurational  free  energy/mole  AG  of  the  three-component 


(1-a-b) ,  x 


B(max) 

solid  solution  is  then  given  by: 


4C(BCC1  -  u0  -  n  (4  (%«<>)„,  .  cAccwC  1>AC  *  cBcc.vC‘^c  , 

♦  Wt!,«  *  *  CBCc“(2,BC  )  } 

*  2  {  ERAC  *A  *  ^BC  xB  *  (ERAC  *  eR0C  '  xa  kb 

.  tyAz  *  4/)  ♦  (yB2  *  O  * 

-  »(2,ab)  OVb  *  Vb)}  •*![  Pia)  ln  Pi<U 


...  do) 


where  E*jj  ■ 


M°ij  -  «(2,ij) 


•“W)  .  O.S  {.I  CjCjEy  •  Jy :  qcjEi^  2  -  cAc8cc  ceachbc, 


T  [EAC  *  E5C  ■  EAB 


...  (11) 


where  ij  -  AB,  AC,  BC  and  -  f 8wC 1 )  -  6w(2)  ] 

J  1  ij  ij  ‘BWG 

k  T^lM/up  .  j(a.ynj  -  .»WK' 


"  ]bwG  «*  “  xfi)  Kg  kB^  Kq  *  XA  *  XB^ }  ' ' 

where  a  -  (kW^J^  <cb  *  CC  *  XA>  «*  -  *A) 

*  f"(!)BC)  ma  «A  *  =0  *  V  tC,  -  Xj) 

-  ("(2)ac  ♦  "U)BC  -  wC!)abJbwg  Kb  -  V  Ka  -  xA)i 

and  xA  and  xB  are  the  values  of  the  order  parameter  in  nearest  neighbour 
positions  at  T  *  Ty. 

S.  Calculation  of  T„  and  T  for  the  Fe-Al-Si  System 

The  following  procedure  was  adopted  to  obtain  optimum  values  of  and 

(BCC  Al-Si).  Examination  of  various  pairs  of  solutes  shows  that  for  BCC  alloys  W^/W^=0.5, 
as  is  indeed  the  case  for  the  other  pairs  concerned  in  the  present  case  (Table  I).  It  was 

therefore  assumed  that  W^/W^  BCC  (Al-Si)  »  0.5.  Inspection  of  the  superlattice  intensities 
Fig.  (5)  show  that  a  considerable  degree  of  order  exists  in  the  B2  lattice  prior  to  DOj 
formation;  the  maximum  theoretical  values  of  x.  and  xg  can  therefore  be  used  in  preliminary 
calculations.  ... 

Trial  values  of  Al-Si^8^  were  *ben  inserted  into  equations  until  a  reasonable 

fit  was  obtained  for  both  T  and  T  for  a  selected  alloy  composition. 

*  /  mi  (2) 

As  shown  in  (Fig.  6)  there  is  a  limit  beyond  which  further  changes  in  W  (and  W  ) 
are  counterproductive,  which  reflects  the  assumption  that  the  real  values  of  and  Xg  are 

lower  than  the  maximum  values  assumed  in  this  preliminary  calculation.  Varying  and  Xg 

in  an  iterative  fashion  can  then  be  used  to  establish  more  realistic  values  of  x^  and  xQ 

as  shown  in  (Table  III).  This  result  can  be  cross-checked  for  other  trial  compositions  and 
a  pattern  established  for  the  variation  of  x^  and  xQ  with  composition.  This  appears  to  be  a 

very  simple  variation,  with  xA  constant  at  0.18  (compared  to  the  initially  assumed 

maximum  value  of  0.25),  and  x„  equal  to  a  constant  fraction  of  CA.  (the  constant  being  the 


ratio  (xA/xA(max))  or  0.72).  These  values  of  Wv  ' ,  Wv  ,  x and  Xg  when  inserted  into 

equations  (11  8  12)  yield  the  results  shown  in  (Table  IV),  and  a  more  direct  comparison  with 
experimental  results  is  given  in  (Table  V),  and  fFig.  7). 
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Fig.  6  Effect  of  varying  Al-Si  on  calculated 

T  and  T  temperatures  using  the  assumption  of  full  A2/B2  ordering  at  T  . 
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6.  Discussion 

The  degree  of  agreement  obtained  is  very  satisfactory,  and  further  calculations  in 
this  system  based  on  these  interaction  parameters  should  yield  equally  reliable  results- 
There  is  some  indication  that  the  equilibrium  values  of  x,  and  xn  for  the  binary  alloys  Fe.Al 

and  FejSi  are  in  fact  slightly  higher  than  would  be  indicated  by  the  simplifying  assumptions 

just  outlined  (see  Table  II)-  Whether  it  is  necessary  to  calculate  equilibrium  values  of  x^ 

and  Xg  vigorously  for  every  composition  instead  of  using  the  simplifying  formula  depends  on 

the  accuracy  with  which  the  ordering  temperatures  need  to  be  known  and  also  to  some  extent  on 

the  cooling  rates  required  to  attain  equilibrium  under  practical  conditions.  Table  VI  gives 

some  idea  of  the  effect  of  small  variations  in  ordering  parameters.  Additionally,  the  values 

of  T  and  T  shown  in  (Fig.  7)  refer  to  the  so-called  constitutive  phase  boundaries  and, 
x  y 

analogous  to  the  T  lines  of  more  common  phase  diagram  calculations,  assume  no  two  phase 

regions  exist.  To  establish  the  full  phase  diagram  a  more  complex  minimization  program  must 
be  utilized  (1,  2,  8)  which  is  rather  lengthy  but  requires  no  further  new  input  data.  It 
should  also  be  noted  that  all  calculations  have  been  made  on  the  basis  that  the  temperatures 
of  interest  are  above  the  magnetic  transition  temperatures  in  the  system;  if  it  is  intended 
to  ir*'estigate  ordering  below  the  Curie  Temperature,  additional  terms  have  to  be  introduced 
into  equations  (10*12).  In  either  circumstance  the  quoted  values  of  wU)  and  can  be 
used  to  derive  heats  of  formation  for  solid  Al-SL  BCC  alloys  which  can  be  used  to  refine  the 
calculation  oi  the  iiquidus  and  solidus  line  in  these  systems,  (17,18). 
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VII  A  SIMPLIFIED  TREATMENT  OF  THE  THERMODYNAMICS  OF  THE  ORDER- 
DISORDER  REACTIONS  IN  THE  Fe-Al ,  Fe-Al-Hn  and  Fe-Al-H i  Systems 

The  Gibbs  energy  difference  between  ordered  and  dis¬ 
ordered  bcc  phase  in  the  Fe  -  A1  system  is  presented  as  an  explicit 
expression  of  composition.  A  power  series  expansion  based  on  the 
shape  of  the  Cp-curve  is  used.  The  thermodynamics  of  the  Fe  -  A1 
phase  diagram  is  evaluated  and  the  diagram  is  calculated  with  the 
new  n  idel  for  the  bcc  phase.  The  Fe-Al  -  Mn  and  Fe-Al  -  Ni 
ternaries  are  calculated  applying  a  subregular  approximation  of 
the  new  series-expansion.  The  deviations  from  previous  calculations 
where  no  distinction  between  ordered  and  disordered  bcc  was  made, 
are  found  to  be  very  small. 

1.  Introduction 

It  is  well  known  that  different  types  of  ordering  (i.e.  magnetic 
or  chemical)  contribute  considerably  to  the  thermodynamic  properties 
of  a  solid  solution.  Consequently  the  corresponding  phase  diagram  might 
be  strongly  affected  by  the  occurance  of  ordering.  Inden  (1)  has 
presented  extensive  theoretical  analyses  of  the  order/disorder  reactions 
in  binary  alloys.  In  a  disordered,  random  solution  the  Bragg-Williams 
approximation  and  different  regular  solution  models  yield  the  Gibbs 
energy  as  a  simple  algebraic  function  of  the  composition  of  the  solution. 

In  a  partially  ordered  phase  the  simple  models  also  predict  a  dependence 
of  the  degree  of  order  on  composition  and  temperature.  The  most  stable 
state  of  ordering  for  a  certain  composition  and  temperature  is  then  the 
one  that  yields  the  lowest  Gibbs  energy.  Thus  it  is  not  possible  to 
write  the  Gibbs  energy  simply  as  an  explicit  function  of  the  composi¬ 
tion.  In  many  applications  it  is  desirable  to  express  the  Gibbs  energy 
as  a  explicit  function  of  composition.  Inden  (2)  has  presented  a  method 
where  the  ordering  contribution  to  Gibbs  energy  is  expanded  in  a  power 
series  of  the  temperature.  The  coefficients  are  functions  of  composition 
and  adjusted  to  fulfill  certain  criteria.  In  the  present  work  a  slightly 
different  method  will  be  applied  to  the  Fe-Al  system.  The  method  is 
based  on  the  work  by  Inden  (3)  and  has  recently  been  applied  to  the 
magnetic  ordering  in  alloys  (4).  The  purpose  of  the  present  work  is  to 
show  how  the  order/disorder  reactions  can  be  incorporated  in  a  Calphad 
framework  without  an  undue  increase  in  complexity. 


A 


2.  The  Model 

The  integral  molar  Gibbs  energy,  G  ,  of  a  phase  can  be  separated  into 
two  parts. 


,dis  . 


G  =  GdlS  +  (G  -GdlS) 
mm  mm 


The  first  term  G^  is  the  Gibbs  energy  for  the  completely  disordered 
state  and  is  obtained  from  a  regular  solution  type  of  model,  i.e.  _ 


XFe  GFe  +  XAl  °A1  +  RT  (  XFeln  XFe  +  XA1  ln  XAl)  +  XFeXAl  LFe 


where  LpeA^  is  a  constant  for  a  regular  solution  and  exhibits  a  linear 
concentration  dependency  in  the  subregular  case.  °^Fe  an<*  are  the 

Gibbs  energies  for  the  pur>e  el'ements  for  the  structure  considered. 

The  second  term  will  now  be  evaluated  by  a  method  which  has  been 
applied  to  the  magnetic  ordering  in  bcc-Fe  (4,5,6),  and  is  a  modification 
of  a  method  initially  proposed  by  Inden(3).  Inden  found  that  the  follow  ig 
expressions  gave  an  adequate  description  of  the  contribution  from  order/dis¬ 
order  reactions  to  the  specific  heat  c  , 

P 


(1+t3) 

k  Rln  (  mF) 


.  .8  .  8m  _  (tS  +  1  ) 

Acp  =  k  Rln  (^5Tf) 


T  <  1 


T  >  1 


where!  =  T/T^,  T^  being  the  critical  temperature  where  all  long  range  order 

vanishes .The  superscripts  a  and  8  denote  the  states  with  long  range  order 

and  short  range  order  respectively,  and,  ka  and  k3  are  constants.  By  applying 

t he  first  three  terms  in  the  Me  Laurin  series  expansion  of  Eq.  3,  and 

3  cl  xs 

integrating  from  T=  °°,  where  G  =  G  ,  to  obtain  the  corresponding  enthalpy 

and  entropy,  we  can  write  down  the  expression  for  the  Gibbs  energy.  It  is 
.  8  ot  . 

convenient  to  express  k  and  k  in  terms  of  the  enthalpy  difference  between 
the  ordered  and  the  completely  disordered  state  at  0  K ,  because  this  quantity 
is  easy  to  obtain  by  counting  the  number  of  different  bonds.  By  introducing 
f  as  the  fraction  of  total  enthalpy  that  is  absorbed  above  Tp  (i.e.,  f=0 

v 

when  there  is  no  short  range  ordering  above  Tc  and  the  total  enthalpy  is  then  due 
to  long  range  ordering  only).  We  obtain: 

-RTC  k3=  J40  f  AHTbJ  0rd  (4a) 

-RTC  k“-  J20  <  1-0  (4b) 
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Inden  (7)  suggests  that  f=0. 4  should  be  used  for  bcc-phases  in 

general.  By  inserting  Eq.  4  in  our  expression  for  G  -G^ls  and  setting 

mm 

f=0.4  we  obtain: 

G  -  G 4Hdi’*ord  I(t) 
m  m  T=u 


Where  0.7089  (t'4/10+t_,4/315+t'24/600) 

I(T)  = ■  T  >  1 

1-1. 1046T+1. 0141  (tV+t'V-.+t16/,.-) 

6  135  600 

T  <  1 


For  B2  ordering  Inden  obtained  7a  and  7b  by  counting  the  number  of 
different  /-bonds; 

,.  .  -RX2  (4w’-3W2) 

...discord  A1 

AH  .  0<Ksli0.5 

T=0 

-R( 1-XA])2(4W'-3W2) 

0.5<Xa1  <  1 


and  the  Bragg-Williams  approximation  yields; 

TC  -  XA1<1~XAI)  <8W!  -  6W2)  (8) 

1  2 

W  and  W  are  the  so  called"exchange  energies,,(expressed  in  k-units)  for  the 
nearest  and  next  nearest  neighbours.  The  Bragg-Williams  approximation  does 
not  account  for  short  range  order  and  consequently  it  predicts  a  critical 
temperature  T^,  which  is  too  high.  Inden  found  that  one  could  introduce  a 
correction  factor  k  to  the  value  calculated  by  means  of  Eq.  8  and  obtain 
good  agreement  with  calculations  based  on  the  cluster-variation  method  which 
do  account  for  short  range  order. 

3.  Application  to  the  Fe-Al  Phase  Diagram 

For  Fe-Al  Inden  (2)  reports  w'=I478  and  W2=794  (k-units,  k=13.8.10  24 
J/atom)  and  k*0.68.  By  applying  these  values  it  is  possible  to  calculate 
the  critical  temperatures  for  B2  and  DO^  ordering.  This  was  done  by 
Miodownik  (8)  who  also  made  a  comparison  with  experimental  information  and 
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found  the  calculated  temperatures  too  low.  To  obtain  agreement  with  the 
experimental  temperatures  at  X^=0.25  he  chose  w'=1588  but  kept  the  other 
values  unchanged.  A  further  investigation  reveals  that  there  are  still 
large  discrepancies  when  XA^<0.25  or  XA^>0.25.  To  obtain  a  better  agreement 
over  the  whole  composition  range  and  must  be  allowed  to  vary  with 
composition.  This  can  be  done  with  a  reasonable  increase  in  effort  and 
good  agreement  with  the  low-temeprature  data  is  then  obtained.  However,  it 
is  found  that  this  procedure  overestimates  the  stability  of  the  bcc-phase 
relative  to  the  liquid  at  high  Al-content  (XA^>0.4)  where  no  information  on 
the  critical  temperatures  is  available.  To  overcome  this  problem  it  is 
necessary  also  to  introduce  temperature-dependencies  m  W  and  W  .  It  is  a 
substantial  increase  in  complexity  to  evaluate  both  a  concentration  and 
temperature  dependency  and  it  was  thus  decided  to  accept  Miodownik's  parameters 
to  avoid  these  difficulties. 

Eqs.  5-10  completely  determine  the  term  G  -G013  in  Eq.  1.  The  remain- 

mm 

ing  thermodynamic  quantities  are  now  to  be  determined  to  comply  with  the  available 
information. 

Kaufman  and  Nesor  (9)  analyzed  the  Fe-Al  phase  diagram  without 
distinguishing  between  the  ordered  and  disordered  bcc-phase.  The  present 
analysis  follows  their  evaluation  except  that  the  ordering  has  been  accounted 
for.  The  total  Gibbs  energy  of  the  bcc-phase  was  obtained  by  adding  the 
ordering  contribution  given  by  Eq.  5  to  the  Gibbs  energy  of  the  disordered 
solution  given  by  Eq.  2.  It  was  then  necessary  to  adjust  all  parameters 
given  by  Kaufman  and  Nesor  and  the  new  parameters  are  summarized  in 
Appendix  A. 

4 .  Calculated  Quantities 

The  recalculated  phase-diagram  is  shown  in  Figs.  1-3.  Fig.  2  shows  a 
magnified  view  of  the  compound  equilibria  and  Figs.  3a  and  3b  show  the 
equilibria  with  the  fcc-phase.  In  Figs.  4a  and  4b  comparisons  are  made  with 
thermochemical  data  from  Hultgren  ( 10)  the  calculations  according  to  Kaufman 
and  Nesor  (9)  and  the  present  analysis.  While  the  agreement  between  the 
enthalpy  of  formation  as  reported  by  the  different  sources  is  satisfactory, 

there  are  larger  discrepancies  for  the  excess  free  energy  of  Al.  However, 
the  experimental  scatter  was  large  in  the  latter  quantity  and  the  new  values 
fall  within  the  error  limits  given  by  Hultgren. 


AD-8136  085  COMPUTER-BASED  METHODS  FOR  THERMODVNAHIC  ANALYSIS  OF 
MATERIALS  PROCESSING(U)  MANLABS  INC  CAMBRIDGE  MASS 
L  KAUFMAN  30  NOV  83  AFOSR-TR-83-1099  F49620-80-C-0020 

F/G  9/2 


UNCLASSIFIED 


(Liquid) 


Calcul 


Molar  Enthalpy  of  Mixing  and  Heat  of  Formation  (kJ) 


Figure  4a. 


Molar  Enthalpy  of  Mixing  of  Iron-Aluminum  Alloys 
and  Heat  of  Formation  of  Compounds  as  a  function 
of  Composition  at  298  K.  The  Crosses  denote  the 
values  given  by  Hultgren  (Reference  10),  the  • 
dashed  heavy  curve  for  the  bcc  phase  is  due  to 
Kaufman  and  Nesor  (Reference  9)  and  the  open 
circles  ior  the  compound  phases  is  due  to  Kaufman 
and  Nesor  (Reference  9).  The  remaining  curves 
and  the  closed  circles  are  from  the  present 
analysis . 
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5.  Subregular  Approximation 

It  is  interesting  to  see  to  what  extent  the  new  formalism  can  be 
approximated  by  a  subregular  model.  The  following  approximations  were 
obtained  for  XAj<0.5  by  plotting  the  quantity  AG/X^X^  and  fitting  a  suit¬ 
able  straight  line  to  the  resulting  curve.  The  following  results  were  obtained 

for  T  >  T  and  XA->0.2 
cr  A1 

AG/XAlXFe  "  2-0,6-l°  T  (  3,75XFe~'  1,49XAl)  (9a) 

XAt  <0.2,  AG  -  0  (9b) 

A1 


For  T  <T 

cr 

AG/XpeXAl  -  R(4W]-3W2)  (o.4444Xpe  -  2.4I26XA1  J 

-R  1.1046  T  (  0.9789Xp,e  -  8.817IXA1  ^ 


-TV-™2)3  (  68-,9Xf.  *  I42xai) 


.  E.bcc  E_bcc 

where  AG  -  G  -  G,..  . 

(dxs) 


(10) 


with  the  actual  values  for  Fe-Al  inserted  we  obtain  (with  A2«disordered  bcc) 
T>Tcr-*A1>0-2: 

G°Cd-^2  -*FeXAl  (  7'455  22.844  I0,5xj  l‘4 


(11) 


T  <  T 


cr 


G°rd  -G*2  *=  X^X^  £  (14669-6. 6 IT- 1.09  10"9  T4)^ 

+  (  -79635  +  59.54T  -  2.27  10~9  T4  )  XAJ 


(12) 
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Both  above  and  below  the  critical  temperature  a  strong  asymmetric  term  is 
obtained.  The  temperature  dependencies  T  ^  and  are  also  strong. 

In  Figs.  Sa,  b,  c  and  d  these  approximation  formulas  have  been  com¬ 
pared  with  Eqs.  5-8.  The  agreement  is  quite  good  except  in  the  transition 
region  between  disordered  and  ordered  state. 

6.  Calculation  of  Isothermal  Sections  in  Fe-Al-X,  X=Mn,Ni 


Once  the  thermodynamics  of  the  binary  solutions  have  been  estabi  ishetl 
the  properties  of  a  ternary  phase  can  be  evaluated  directly,  assumin'' 
that  there  are  no  specific  ternary  effects.  In  a  rigorous  treatment  the 
order  disorder  contribution  is  still  given  by  eqs.  5-6,  but  AH  and  T 
are  now  to  be  evaluated  for  the  actual  ternary  composition.  It  is  a 
straight  forward  procedure  to  derive  expressions  for  AH  and  Tc  in  the 
ternary  case  but  the  final  expressions  are  rather  complicated.  In  the 
present  case  a  simpler  approach  will  be  taken.  The  Kohler  method  has 
been  applied  frequently  to  synthesize  the  properties  of  a  ternary  solution 
when  the  binaries  are  known,  see  for  example  ref. 11  .  Admittedly  there 
is  no  physical  justification  for  this  method  in  the  present  case  and  its 
general  features  have  recently  been  questioned^).  However ,  the  use  of  the 
method  is  widespread  and,  for  the  sake  of  simplicity,  it  was  applied  here 
as  follows:  First  a  ternary  calculation  was  made  using  the  subregular 
model  assuming  a  completely  disordered  bcc-phase  described  by  the  Kohler 
method.  Thereafter  a  new  calculation  was  performed  with  the  order-disorder 


contribution  added.  When  X, 


<  0.20  the  latter  contribution 


<W 

is  approximately  zero,  as  can  be  seen  from  Fig.  5,  and  the  phase  boundaries 
follow  the  first  set  of  curves.  When  Al^^Fe^Al^  °*30  the  subregular 
solution  approximation  is  a  good  approximation  of  the  order-disorder 
contribution  and  the  phase  boundaries  thus  follow  the  second  set  of  curves. 
The  curves  in  the  transition  region  were  drawn  by  hand  assuming  a, gradual 
shift  from  one  set  of  curves  to  the  other  set. 

All  parameters,  except  those  for  Al-Fe,  were  taken  from  the 
ManLabs  data  bank  (14).  The  results  are  presented  in  Figs.  6-7 


k..  - 
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Figure  5.  Calculated  regular  solution  parameter  for  bcc 

Iron-Aluminum  alloys  as  a  function  of  composition 
at  various  temperatures.  The  thick  curve  represents 
the  present  analysis  while  the  thin  segments 
represent  the  subregular  approximation. 
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Figure  6.  Calculated  fcc/bcc  Equilibria  in  the  Iron-Manganese 
Aluminum  System  as  a  function  of  temperature  between 
ordered  bcc  and  fee  solid  solutions.  The  dashed  curve 
shows  the  phase  boundary  for  the  disordered  bcc  phase. 
No  difference  was  noted  above  1173°K. 
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Figure  7 


Calculated  fcc/bcc  Equilibria  in  the  Iron-Nickel- 
Aluminum  System  as  a  function  of  temperature  between 
ordered  bcc  and  fee  solid  solutions.  The  dashed 
curve  shows  the  phase  boundary  for  the  disordered 
ase . 


It  can  be  seen  that  the  order-disorder  contribution  to  the  Gibbs 
energy  generally  yields  a  very  small  shift  of  the  calculated  phase  boundaries. 
Furthermore,  a  comparison  with  the  calculations  by  Kaufman,  (13)  who  did  not 
distinguish  between  ordered  and  disordered  state,  shows  that  there  are  no 
significant  discrepancies.  The  curves  fall  very  close  to  each  other. 

It  will  now  be  argued  that  the  present  procedure  probably  over¬ 
estimates  the  effect  of  ordering  on  the  bcc/fcc  equilibrium.  For  a  given 
temperature  the  composition  range  within  which  ordering  occurs  is  approximately 
given  by  the  parabola:  X^  Xpe  =  constant  which  starts  from  the  binary  Al-Fe-sid* 
and  loops  into  the  ternary.  The  exact  relation  is  complicated  (8)  and  requires 
detailed  knowledge  about  the  order-disordering  reactions  in  the  other  two 
binaries.  The  parabola  is  a  reasonable  approximation  when  the  maximum  ordering 
temperatures  are  low  (but  still  above  0°K)  in  the  other  systems.  In  Fig.  8 
the  parabola  has  been  plotted  at  different  XA1  Xpe  values.  When  these  plots 
are  superimposed  on  the  calculated  bcc/fcc  equilibria  it  is  found  that  all 
temperatures  the  bcc/fcc  phase  boundary  falls  well  outside  the  field  of  ordered 
bcc.  It  is  thus  concluded  that  the  true  phase  boundaries  must  follow  the  ones 
calculated  with  a  disordered  bcc  phase. 

As  a  final  conclusion,  it  would  be  worthwhile  to  examine  the  degree 
of  order  at  each  composition  in  the  bcc  phase  at  the  two  phase  bcc  +  fee  boundary 
For  consistency  this  degree  of  order  should  coincide  with  the  degree  of  order  in 
the  iron-aluminum  binary  at  the  same  X^/X^  ratio  which  characterizes  the 
ternary  composition.  If  this  is  not  the  case,  then  the  ternary  free  energy 
synthesized  from  the  binary  components  incorrectly  reflects  the  wrong  degree 
of  order.  In  the  present  case  the  errors  introduced  by  such  an  inconsistency 
do  not  appear  to  be  very  large.  However,  it  is  possible  that  such  errors  ^could 
be  substantial  in  other  cases  and  should  be  guarded  against. 

Note  in  Proof 

The  prior  discussion  of  Section  3  drew  attention  to  the  discrepancy 

between  the  experimental  ordering  temperature  and  the  calculated  value  for 

0.25  <  X.,  <  0.50.  It  was  noted  that  compos  it ionally  dependent  values  of 
—  A1 
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V  and  W  would  be  required  to  eliminate  this  discrepancy.  It  was  also  noted 
that  such  an  alteration  would  lead  to  a  stabilization  of  the  bcc  phase  at 
high  temperature  and  a  resulting  excessive  elevation  of  the  melting  points. 
Notwithstanding  these  difficulties  an  attempt  was  made  to  evaluate  W  and  W4 
as  a  function  of  composition  in  order  to  produce  an  ordering  temperature  vs. 
composition  curve  for  0.25  1  XA1  1  0.50  (8)  corresponding  to  1260K  at  X  = 
0.30  and  1570K  at  XAj“  0.40.  The  resulting  expression  for  W '  is  given  by 
Equations  13  and  14  as: 

W1  -  1480  +  625  XA1  ( 13) 

and 

IT  -  846  -  372  XA1  ( 14) 

1  2 

These  new  values  compare  with  the  previous  values  of  W  =  1478  and  W  =  794 

listed  in  Section  3.  These  expressions  for  W  and  W  as  provided  by 

Equations  (13)  and  (14)  were  substituted  directly  into  Equations  (5)  through 

(8)  to  reformulate  the  thermochemical  properties  of  the  iron-aluminum  system. 

No  further  changes  in  the  remaining  parameters  for  the  iron-aluminum  system 

were  made  at  this  stage.  This  procedure  did  indeed  produce  an  excessive 

stabilization  of  the  bcc  phase  with  respect  to  the  liquid.  Nevertheless, 

this  new  formulation  of  the  bcc  phase  was  employed  directly  to  derive  an 

effective  subregular  description  as  in  Section  5.  The  latter  step  was  accom- 

plished  by  evaluating  W  and  W  at  XAJ  -  0.5  from  Equation  13  and  14  since 

at  XA1  "  °‘5  the  largest  difference  between  this  approximation  and  the 

former  calculation  exists.  Thus  for  this  case  w'  *»  1793  and  W2  «»  660. 

Employing  these  values  in  Equations  (9)  and  (10)  permits  redefinition  of  the 

excess  free  energy  of  iron-aluminum  alloys  in  the  subregular  approximation. 

Subsequent  recalculation  of  the  Fe-Ni-Al  fcc/bcc  equilibrium  at  1273  with  the 

modified  results  did  not  yield  results  significantly  different  from  those 

shown  in  7c.  It  is  therefore  concluded  that  the  results  described  in  .Section 

6  would  not  be  materially  altered  by  introducing  compositionally  dependent 
12 

values  of  W  and  W  for  the  purpose  of  obtaining  a  better  description  of  the 
ordering  temperature  at  compositions  containing  more  than  25  atoms  percent 
aluminum. 
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Appendix  A:  Summary  of  Parameters  for  the  Fe-Al  System 
The  lattice  stabilities  from  Ref.  9. 

E  bcc.dis 


m 


*FeXAl  ^“32345-37 .92T)Xpe  +  (-150252+51 .42T)XA1  J/mol 


1  2  -24 

IT  =  1588,  VT  =  794  (k  units  k=13.8-10  J/atom) ,  K=0.68 


EGfcc  =  X„  X  -  (-25022-38. 01T)X_  +  (-106841+66 . 70T)XA . 

m  re  A1  Fe  A1 

EGL  =>  X,,  X..  (-62760-14 .63T)X  +  (-96232+32. 72T)X41 

m  Fe  A1  Fe  A1 


J/mol 

J/mol 


0  °  bcc  °  bcc 

5Fe  Al  ~  °‘25  GFe  -°*75  GA1  =  “35275+8. 28T 
0.25  0.75  A1 


J/mol 


*  n  °  bcc  °  bcc 

G  '  A1  -  0.286  G°c c-0. 714  G°,c  =  -41169+11 .40T 

Fe0.286A10.714  Fe  A1 


J/mol 


Fe0.40A10.60 


,bcc 

JA1 


Gp^C-0.667 

GA1C  ■  -37872+4. 73T 

J/mol 

0  her* 

-0.6  G” 

=  -23863-2. 77T 

J/mol 

r ' 
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VIII  A  THERMODYNAMIC  EVALUATION  OF  THE  TI-C-N  SYSTEM 


1 . 0  Introduction 

The  use  of  vapor  deposited  .LC,  TiN  and  Ti(C,N) 
coatings  for  improving  the  performance  of  cutting  tools  and  wear 
resistant  surfaces  enhances  the  value  of  information  on  the 
thermochemistry  and  phase  diagram  of  the  Ti-C-N  system. 

Since  the  crystal  structures  of  the  compounds  TiC  and 
TiN  are  the  same,  it  is  likely  that  these  compounds  form  a  solid 
solution.  However,  the  experimental  Ti-C  and  Ti-N  phase  diagrams 
reported  to  date  are  not  very  well  established  and  very  little 
thermodynamic  data  is  available  for  the  binary  cases  or  the 
Ti-C-N  ternary.  The  purpose  of  the  present  work  is  to  evaluate 
all  the  information  available  in  terms  of  thermodynamic  models 
for  the  individual  phases  and  synthesize  the  binary  and  ternary 
phase  diagrams. 

2 . 0  Thermodynamic  Models 
2 . 1  Liquid 

The  simplest  models  that  cover  the  whole  concentration 
range  are  the  regular  or  subregular  solution  models.  The  formation 
of  the  solid  compounds  TiC  and  TiN  in  the  binaries  suggests  that 
some  ordering  or  formation  of  molecular  species  might  take  place 
in  the  liquid  around  the  stoichiometric  compositions.  However, 
a  high  degree  of  association  is  not  observed  in  the  gas  phase  and 
vaporization  of  titanium  carbide  yields  titanium  and  carbon  polymers 
The  Fe-S  system  has  recently  been  treated  by  assuming  the  presence 
of  molecular  species  (1)  and  by  assuming  ordering  (2).  Botk 
treatments  yield  a  satisfactory  agreement  with  experimental  data, 
and  predict  a  very  rapid  change  in  activity  close  to  the 
stoichiometric  composition.  This  rapid  change  is  difficult  to 
describe  with  a  regular  or  subregular  solution  model,  and  generally 
many  terms  will  be  needed  in  the  series  expansion  of  the  excess 
energy. 


The  abnormal  behavior  of  the  activity,  compared  to  what 
the  regular  solution  model  predicts,  reveals  itself  in  the  phase- 
diagram  by  giving  the  compound  a  liquidus  curve  formed  more  like 
a  wedge  instead  of  a  parabola. 

In  the  present  case  the  experimental  data  on  the  L/TiC 
(liquid/TiC)  and  L/TiN  equilibia  is  too  meager  to  reveal  the  shape 
of  the  liquidus  curves  accurately,  and  no  thermodynamic  data  for 
the  liquid  phase  is  available.  For  the  sake  of  simplicity  the  liquid 
was  thus  described  with  the  subregular  solution  model.  The  cor¬ 
responding  equations  are  summarized  in  Appendix  b. 

2 . 2  Solid  Phases 

The  Tie  and  TiN  phases  both  have  a  Nacl-structure .  However, 
a  rather  wide  range  of  homogenuity  is  observed:  Ti^  X)CX  with 
0.32<  x  <  0.49.  First  an  attempt  was  made  to  describe  TiC  as  a 
subregular  solution.  This  attempt  was  not  very  successful  and  a 
new  approach  was  tried.  From  the  crystallographic  point  of  view 
TiC  and  TiN  consists  of  two  sublattices  and  the  variation  in 
composition  is  due  to  a  varying  number  of  vacant  carbon  or  nitrogen 
sites.  (3).  It  was  thus  decided  to  apply  the  sublattice  model 
which  was  developed  by  Hillert  and  Staffansson  (4)  and  later  gen¬ 
eralized  by  Sundman  and  Agren  (5)  and  recently  applied  to  the  Fe-S 
system  (2).  This  approach  was  more  successful  and  yielded  a 
satisfactory  description  of  all  the  available  information. 

Hexagonal  titanium  exhibits  a  rather  large  solubility  for 
N  and  the  choice  of  thermodynamic  model  might  be  important.  Recently 
(6)  the  hexagonal  solution  of  N  in  Fe  was  treated  by  means  of  the 
sublattice  model.  In  agreement  with  that  treatment  we  chose  the 
upper  limit  of  N-solubility  in  hcp-Ti  as  1/3  corresponding  to  the 
compound  T^N.  For  the  bcc-phase  we  also  applied  the  sublattice 
model  with  3  interstials  per  Ti-atom  as  in  the  Fe-C-case  (see  Ref. 6) 

The  equations  for  the  sublattice  model  are  summarized  in 
Appendix  B . 


3.0  Evaluation  of  Parameters  for  Phases  in  the  Titanium-Carbon 


The  free  energy  of  formation  of  stoichiometric  TiC  has 
been  tabulated  by  Hultgren  (3).  Figure  1  shows  the  data  given 
by  Hultgren  as  symbos  and  the  solid  line  is  a  least  square  fit 
(assuming  a  linear  temperature  dependency)  to  the  data.  The  dashed 
line  is  calculated  according  to  the  parameters  given  by  Kubaschcwski 
and  Alcock  (7).  The  result  of  the  least  square  fit  is: 

°nTic  ’  °GTiC'  °Gcr  =  -45012  +  3.428T  (cal/mol)  fl) 

and  was  accepted  in  the  present  study.  Storms  (S)  has  compiled  data 
on  the  Ti  partial  pressure  over  TiC  at  equilibrium  with  graphite  at 
different  temperatures  (See  Fig.  2)  Kubaschewski  and  Alcock  (7)  give 
the  vapor-pressure  of  pure  Ti  as: 

log  P°Ti  =  -  24400/T  -  0.91  log  T+10.299  f2) 

From  Eq?  we  can  evaluate  the  corresponding  free  energy  difference; 

°G|?S»  latm  -  °G^C  =  11164S  -  47.124T  +  1.808  TlnT  (cal/mol) (3) 
By  combining  the  information  from  Storms  (8)  with  Eq.3  and 

°G^C  -  °G^9c  *  -240  +  0.9T  (cal/mol)  (4). 

from  Kaufman  (9)  and  take  the  composition  of  Tie  in  equilibrium 
with  graphite  as  xr  =  0.49  according  to  Rudy's  phase  diagram  (10) 
we  can  calculate  the  partial  excess  quantities  GT.  and  G,  as 
functions  of  temperature  for  this  particular  composition.  The  phase¬ 
boundary  Ti/TiC,  as  reported  by  Rudy  (10),  can  be  used  to  evaluate 
GTi  at  different  temperatures  along  the  phase-boundary  when  the 
solubility  of  C  in  Ti  is  neglected.  Assuming  linear  temperbture- 
dependency  it  was  then  possible  to  evaluate  °h£ygCC  and 
from  the  partial  quantities.  The  following  expressions  were  obtained: 

°LCVaCC=  ~16099  +  0.36T  (cal/mol)  (5a) 

lLCVaCC  s  “38803L  +  H-OST  (cal/mol)  (5b) 

The  solid  line  in  Figure  2  is  calculated  by  the  application  of  F.qs. 

5a  and  5b.  The  agreement  with  the  experimental  information  is 
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The  Free  Energy  of  Formation  of  TiC  from  bcc  Ti 
and  Graphite  as  a  function  of  Temperature.  Solid 
Curve  designates  the  present  analysis  while  the 
Dashed  Curve  is  after  Kubaschewski  and  Alcock  (.Re! 
and  the  points  are  from  Hultgren  (Ref.  3) 


Figure  2. 


Figure  3. 


3.8  4.2  4.6  ( 1 0 , GQO/T°K) 


Pressure  of  Titanium  over  TiC  +  Graphite  as  a  function 
of  Temperature.  Solid  curve  is  from  the  present 
analysis,  points  are  from  Ref.  7. 


Yc 


Pressure  of  Titanium  over  TiC  at  190Q°K.  Solid 
Curve  is  from  Present  Analysis.  Points  are  from 


very  good.  In  Figure  3  the  Ti  partial  pressure  as  a  function  of 
composition  at  1900  K  has  been  plotted.  The  symbols  are  experimuit 
reported  by  Storms  (8)  and  the  solid  line  is  calculated  from  the 
present  analysis.  The  agreement  is  sat isfactory . 


The  position  of  the  melting-point  of  TiC  according  to  Rudy 
(10)  is  *£  =  0.44  and  =  3340k.  By  applying  the  new  evaluation  nl 
TiC  and  the  lattice  stabilities  reported  by  Kaufman  (9)  one  can 
evaluate  the  subregular  solution  parameters  that  yield  this  melMny- 
point.  Furthermore,  temperature  dependencies  in  these  paranete  n- 
can  be  evaluated  from  the  condition  that  the  L/TiC/Gr  eutectic 
occurs  at  3049K  and  the  L/TiC/Ti  eutectic  at  1923K  (from  Ref.  K' .  : 
The  following  values  are  then  obtained: 


°LCTi  =  _23185  "  10. ST 
lLLCTi  =  30516  "  10*5T 


( ca  1/mo  1 )  ( (> a  i 
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( r.h  i 


The  sublattice  model  for  the  bcc-phase  yields  a  reasonable  descript  i 
with  the  following  parameters: 


LCVaCC  =  “50000 


(ca 1/mol )  ( 7a i 


o~bcc  o^hep  <7  o,,gr  .  . ^ ^ 

GTic  “  r’Ti  “  3  (,c  =  ~44996  *  0 . 49T 


The  hep-phase  has  a  very  low  solubility  of  carbon  and  it  is  thus 
difficult  to  evaluate  the  interaction  parameter.  Arbitrarily, 
this  was  chosen  the  same  as  in  bcc  i.e. 


.hep  .bcc 

LCVa  LCVa 


-50000 


(cal /mol)  ( .4  a 


Subsequently 

oJicp  _  onhcp  ,  /0o..Gr  . 

'Vl^N  -  2  CTi  1/2  CC  "5574 

was  chosen  to  yield  a  three-phase  equilibrium  fcc/bcc/hcp  at  *1193K. 

All  parameters  are  summarized  in  Appendix  B. 

4 .  Evaluation  of  Parameters  for  Phases  in  the  Titanium-Nitrogen 
System. 

In  general  the  data  for  the  TiN-System  is  meager.  The 
phase  diagram  given  by  Hultgrcn  (3)  is  not  very  well  established 
and  is  basically  the  same  as  the  one  in  Hansen’s  compiliation  (11). 


■  f  - .  -t  - .  -  *.  . 

.  .  .  VstV.N-.vV.y.^. 
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The  free  energy  of  formation  for  stoichiometric  TiN  however,  seems 
well  known.  A  least  square  fit  to  the  values  given  by  Hultgren  '5 

yields : 

orfcc  orhcp  . /7  oreas 

’TiN  "  °Ti  *  L/ L  “N 2  =  -"’9666  ♦  21.42T  rCal/mol  '9 

This  expression  yields  values  in  good  agreement  with  those  calculated 
for  Kubaschewski  and  Alcock  (7).  In  addition  some  information  on  the 
enthalpy  of  mixing  in  the  solid  phases  is  available  from  Hultgren. 
Hultgren  also  gives  the  free  energy  of  mixing  for  0.45  <  x^.O.S  at  1930K. 

A  solution  parameter  *  -25000  +  4.4T  was  found  to  describe  the 

information  on  the  fee  phase  satisfactory.  In  Fig.  4  the  calculated 
values  have  been  plotted  together  with  Hultgren’s  values.  The  parameter.- 
for  the  hep-phase  t;as  subsequently  evaluated  from  the  fcc/hcp  phase 
boundaries  as  reported  by  Hultgren.  The  following  values  were  obtained. 

^’TijN  “  2  °GTiP  ‘  1/2  °r,§aS  ■  -95490  ♦  24. 7  IT  U0a 


LNVaCP  "  "998S  -3*44T 


r  iob  i 


The  enthalpy  of  mixing  in  the  hep  phase  can  now  be  calculated  as  a 
function  of  composition  and  the  result  is  presented  in  Fig.  6  together 
with  two  experimental  values  given  by  Hultgren.  The  agreement  is 
satisfactory. 

In  Hultgren's  phase  diagram  a  tetragonal  phase  Ti2N  with  a 
very  narrow  composition  range  is  shown  at  low  temperatures.  This 
phase  was  treated  as  a  stoichiometric  compound  in  the  present  evaluation 
and  it's  free  energy  of  formation  was  evaluated  from  the  reported 
phase  boundary  with  hep  at  1310K  and  1200K,  yielding  the  following 
expression : 

Gt  i  »i  —  2  G— .  —  1/2  °Ggas  »  -120692  ♦  43.06T  (cal/mol1 

Ti,N  1 l  N, 


The  highest  temperature  where  this  phase  is  stable  is  then  calculated 
as  1321K. 

It  is  very  difficult  to  evaluate  two  parameters  for  the  bcc- 
phase  because  the  N-solubility  is  rather  low.  Quite  arbitarily  it 
was  then  chosen  to  put: 

i T ibcc  n  c  i  7  i 

LNVa  "  0  fl2) 

The  remaining  parameter  was  then  adjusted  to  yield  a  hcp/bcc  equilibrium 

in  reasonable  agreement  with  Hultgren's  phase  diagram. 


The  following  value  was  obtained: 


o-bcc  orhcp 
nTiN3  "  GT i 


3/2 


#«s" 


-  104  155  *  31.22T 


( 13) 


No  thermodynamic  data  for  the  liquid  phase  is  available.  The  two 
subregular  solution  parameters  were  adjusted  to  yield  an  fcc/L 
equilibrium  in  reasonable  agreement  with  Hultgrens  phase  diagram. 
The  following  parameters  were  found: 


°^NTi  *  "64482  -7.5T  (cal/mol)  ilia) 

HJJti  =  -40975  ♦  7. ST  fcal/mol)  Q4bl 

All  parameters  are  now  fixed  and  a  complete  phase  diagram  can  be 
calculated.  The  parameters  are  summarized  in  Appendix  B. 


5 . 0  Calculation  of  Tl-C  Phase  Diagram 

Combination  of  the  equations  in  Appendix  B  and  the  parameters 
in  Appendix  C  permit  calculation  of  the  phase  diagram  by  standard 
methods.  The  result  is  presented  in  Fig  S. 

Figure  6  compares  the  experimental  points,  the  phase  boundnrie 
(dashed  curve)  suggested  by  Rudy  (10)  and  the  calculated  results 
Isolid  lines).  The  agreement  with  the  suggested  solidus  curve  is 
very  good  whereas  the  agreement  with  the  suggested  liquidus  curve 
is  not  so  good.  However,  the  liquidus  curve  has  not  been  established 
experimentally  to  date.  It  should  however,  be  mentioned  that  the 
appearence  of  order  or  molecular  species  in  the  melt  will  deform  the 
liquidus  curve  to  a  shape  more  like  the  one  proposed  by  Rudy.  HOI. 


6.0  Calculation  of  the  Ti-N  Phase  Diagram 

The  calculated  titanium-nitrogen  phase  diagram  is  shown 
in  Fig  7.  In  addition  Hultgrens  phase  diagram  is  displayed  using 
dashed  lines.  The  solid-phase  equilibria  are  in  a  very  good  agree¬ 
ment  and  the  fcc/L-equilibrium  is  satisfactory.  In  both  Hansen's 
and  Hultgren's  phase  diagrams,  the  latter  equilibria  was  designated 
by  dashed  lines  to  indicate  that  the  shape  of  the  curves  is  very 
uncertain.  The  most  serious  discrepancy  between  the  reported  phase 
diagram  and  the  calculated  one  is  the  temperature  of  the  three-phase 
equilibrium  L/bcc/hcp. 
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Figure  6.  Liquid-Fee  Equilibrium  in  the  Titanium-Carbon 
System.  Solid  Curves  are  calculated.  Dashed 
Curves  and  Experimental  Points  are  from 
Reference  9. 
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In  Hultgren’s  diagram  it  is  around  2293K  (2020  °C)  whereas  the 
calculation  yielded  2037k  (1764°C).  It  is  easily  demonstrated 
that  this  part  of  Hultgren's  diagram  cannot  be  self  consistent 
since  close  to  the  edges  of  a  phase-diagram  the  width  of  a  two- 
phase  field  is  approximately  given  by  AG/RT,  where  AG  is  the 
free  energy  difference  between  the  two  phases  for  the  solvent 
(i.e.  Ti  in  this  case).  At  2293K  we  have  for  the  width  of  the 
L/bcc-f ield , :  AG/RT  =  0.15,  whereas  Hultgrens  diagram  yields 
about  0.02.  A  two-phase  field  which  is  wider  than  0.02  will 
tend  to  push  the  three-phase  equilibrium  to  lower  temperature. 
Altering  the  width  of  the  two-phase  field  by  changing  solution 
phase  parameters  cannot  be  done  realistically  because  the  nitrogen 
concentration  in  the  liquid  is  too  low.  It  might  instead  be 
possible  to  change  the  description  of  the  bcc  phase  but  the  N- 
concentration  in  that  is  also  rather  low  and  non-realistic  values 
are  required  to  obtain  agreement  with  Hultgren.  In  the  latter 
case  one  also  has  to  change  the  hep  and  fcc-phases  if  one  wants 
to  keep  the  agreement  with  the  solid-phase  equilibria  in  the 
reported  phase  diagram. 

7 . 0  Calculation  of  Ternary  Sections  in  the  Ti-C-N  system 


In  general  calculation  of  the  Ti-C-N  system  requires 
description  of  the  C-N  binary  in  addition  to  the  Ti-C  and  Ti-N 
binaries  discussed  above.  However,  since  the  information  required 
to  evaluate  such  C-N  interactions  is  not  available,  they  were 
assumed  to  be  equal  to  zero  as  a  first  approximation.  On  this 
basis  the  ternary  sections  at  1600,1800  and  2500K  for  the  N-C-Ti 
system  shown  in  Figures  8,9,10  were  calculated.  The  calculation 
shows  that  the  fcc-field  extends  from  the  TiC  to  the  TiN  side  and 
that  the  hep  phase  disolves  very  little  carbon. 

8.0  Vapor  deposition  of  solid  Ti,.  -.C  N 
_ _ _ (1-x-y)  x  y 

From  a  practical  point  of  view  it  is  interesting  to 
determine  what  composition  an  input  gas  must  have  in  order  to 
yield  a  specific  composition  x,y  of  the  carbo-nitride .  The 
limiting  values  at  1800K  in  the  binary  Ti-C  system,  using  a 
mixture  of  TiCl^,  CH^  and  of  Ptot  =  1  atm,  was  calculated 
by  Bernard  et.  al.  (12).  In  principle  their  calculation  could 
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culated  Isothermal  Section  at  1600°K  in  the 
anium-Carbon-Nitrogen  System. 
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be  repeated  not  only  for  the  extreme  compositions  (corresponding  to 
equilibrium  with  graphite  and  almost  a  pure  Ti  respectively)  but  for 
a  whole  family  of  different  compositions.  Each  composition  would  be 
represented  as  a  curve  in  the  P^^ci  -  P0^  “  plot.  The  calculation 

is  quite  straight  forward.  From  the  given  solid  phase  composition 
the  C-  and  Ti-activities  can  be  calculated.  By  keeping  these  activities 
and  the  total  gas-pressure  fixed  and  solving  the  non-linear  system 
of  equations  formed  by  equilibrium  and  stoichiometric  conditions  for 
the  gas-mixture  one  obtains  P^ici  when  P°^  is  specified  or  vice 

versa.  If  a  nitrogen-bearing  component,  e.g.  NH^,  is  added  to  the 
input-gas  the  deposited  phase  will  also  contain  some  N.  The  input- 
gas  composition  required  to  obtain  a  certain  solid-phase  composition 
is  calculated  in  the  same  way  as  in  the  binary  case.  From  the  solid 
phase  compositon  the  C-,  N-  and  Ti-activities  are  calculated  and  the 
values  obtained  are  inserted  to  the  equilibrium  and  stoichiometric 
conditions.  In  this  case,  however,  each  solid  phase  composition  is 
not  represented  by  a  curve  in  the  P°icl  -  p£h  -  plot,  but  by  a 

surface  in  the  P°icl  -  P°H  -  Pj^  -  plot.  In  Figures  11-13  the 

4  4  3 

calculated  iso-activity  lines  in  the  fcc-phase  have  been  ploted  for 
each  component  at  1800K  so  that  the  activities  can  be  estimated  easily 
for  a  certain  composition  at  this  temperature.  As  expected,  the 
nitrogen  and  carbon  activities  required  to  obtain  a  solid  phase 
-without  vacancies,  i.e.  x^^  =  0.5,  are  infinitely  high. 
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Fiaure  11.  Calculated  Iso-activity  lines  for  Titanium 
in  Ti(CfN)  at  1S00K 
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Appendix  B 


For  a  ternary  regular  solution  we  apply  the  following 
equation  for  the  integral  Gibbs  energy: 

V  xA°GA  *  XB  °r’B  *  xc  °GC  *  RT  [XA  lnxA  *  xBlnxB 
*  XClnXJ  *  xAXB  LAB+  xAxC  LAC  *  XBXC  LBC  (A1) 

where  xa,xb,xc  are  tbe  no^ar  fractions  and  subjected  to  the 
condition  : 

XA*  XB  *  XC  *  1  CA2) 

GA  ,  gb,  °GC  are  the  Gibbs  energies  for  the  pure  components 
in  the  phase  under  consideration,  and  LAB,  LAC ,  Lg(~  are  inter¬ 
action-parameters.  In  a  subregular  solution  the  interaction 
parameters  exhibit  a  linear  concentration  dependence  in  the 
binary  edge-systems.  In  this  report  the  following  representation 
has  been  applied  in  the  binaries: 

LAB  '  °LAB  *  lLAB  (XA-V  (A3) 

The  same  representation  will  also  be  applied  for  a  ternary  solution. 
By  standard  methods  the  partial  quantities  are  derived: 


where 


nj"°Gj*  RT  lnxj*F'r'j 


j “A,B ,C 


nA*  XB  (1’XA}  °LAB  '['I  *  2xb  (l"xA')(xA"xB)]  l'L; 
♦  xc(l-xA)  °lAC  *£x2  ♦  2xc(1-xa)(xa -xc)JLLAc 


XBXC  LBC 


2xbxc^xb"xc.^  lbc 


CB  ’  XA(l’xB)  °LAB  U  *  2xA(1'xB)(xA~xB)J  \b  * 
*xc(1’xb)0lbc  *  Cxc  *  2xc(1“xb)(xb’xc)J  1lbc 
*  xaxc  °lac  -  2xAxc  (xa“xc)  llAC 

Ecc  "  XA(1-XC)  °LAC  [xA  •  2xA(1-xcHxa-xC]]  1Lac 

+  xb(1-xc1<5lbc  -£xb  •  2xB(l'xc1(xB'xc)J  lbc 

”XAXB  °LAB  -  2xAxBfxA'xB)  ^AB 


The  sublattice  model  takes  the  following  form  for  the  case 
of  Ti  on  one  sublattice  and  a  mixture  of  C,  N  and  Va  (=vacant  sites) 
on  the  other  sublattice  (i.e.  Ti  (C,N,Va)  ) 

cL  C 

Gm  a  yC  °GTiaCc  +  yN  °GTiaNc  +  yVa  a°GTi 

+  yCyN  LCN  +  yC  yVaLCVa 
Ti 

+  yNyVaLNVa  (A8) 


+  RT  c 


[yclny< 


yNmyN 


Va 


lnyVa] 


where  the  site  fractions,  y  can  be  calculated  from  the  ordinary 
mole  fractions  by  the  equations: 


yc  -  2  Jtc/(1‘xc  "  V 


yN  ■  2  V(1'xc  "  xn> 

C 


( A9 ) 
(A10) 


and 


yVa 


:  *  1- 


yC  “  yN 


(All) 


The  subscript  m  in  equation  A8  signifies  that  G  is  counted  per 
mole  of  formula  unit:  Ti  (C,N,Va)  .  In  order  to  compute  G  for 
a  mole  of  atoms  (or  a  gram  atom)  equation  A8  must  be  multiplied 
by  the  factor  (1/a)  (l-xN-xc)  .  The  following  values  of  a  and  c  hold: 

fee:  a= 1 ,  c= 1 
bcc:  a=l ,  c=3 

hep:  a=*2,  c  =  l  ( A 1 2 ) 

°GX.  „  and  °GT.  „  are  the  Gibbs  energy  of  formation  for  the 
* 1  C  » i„N 

a  c  a  c 

compounds.  These  compounds  may  or  may  not  be  stable.  Moreover, 
experimental  data  may  or  may  not  be  available  for  the  compounds. 


The  L-parameters  are  interaction  parameters  and  are  defined 
to  vary  with  composition  as: 


L  =  °L 
Lij  L 


1L. .  (y .  -  y, ) 
ij  i 


(A13 ) 


% 


s-v.'  VVV  VV*- 'v’v 


■  ...  1 _ ■ 
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The  partial  quantities  become: 


Gc  -  1  (°GTi  c  -  a  °GTi)  +  RT  (hyc 
C  3  c 

*  eg„ 


Gn  -  I  C  0Gt1  .n  -  a  °GT1)  *  RT  (Iny 
c  a  c 

♦  EGVI 


gT1  ■  °GTi  *  RT  5  lnYVa  *  E(:Ti 


Furthermore ; 

E(3C  i  f  ^GliaCr  "  E^Tia  Va J 
c  a  c  a  c 


>"yVa> 


-  myVa) 


Ep  1  fEp  Er  \ 

GN  c  C  GTi.N  "  GTi  VaJ 
a  c  a  c 


EGTi  =  i  Si  Va 
ii  a  a  c 


(A13 ) 


(A14) 


(A15) 


ERiiacc  -  yN  (1-yG}  °lJn  -K +  2yN(1-yc)(yc  -xnO  1l 

yVa(1"yC}  LCVa  +  CyVa  +  2yVa(1“yC)  <-yC_yVa')J  llCVa 
yNyVa  °LNVa  "  2yNyVa  (yN  '  yVa}  ^NVa 

GTiaNc  =  yC(1“yN5  °LCN  '  [yC  "  2yC(1"yN^(yC"yN0  ^CN 
+  yVa(1'yN}  °LNVa  +  fyVa  +  2yVa  ^1-yN^  (yN  "  yVaO  ^N' 


1T  Ti 
lCN 


( A16  ) 


l.Ti 

bNVa 


yCyVa  LCVa  "  2yCyVa  ^yC"yVa^  ^CVa 


i  (A17) 


EGTi -Va  3  Yc  (1’yVa5  °LCVa  “  [yC  “  2yC(  1”yVa)  ^yC‘yVa*.j|  N 
av  c 

♦  yN  (1_yVa^  LNVa  “  CyN  "  2yN  ^  1-yVa^yN"yVa-)]  *LNVa 


yCyN  °LCN  “  2yCyN  ^yC“yN^  1lCN 


(A18 ) 
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Summary  of  parameters  for  Ti-C-N  system,  (cal/mol,  conversion 
to  J/mol  by  multiplying  by  4.184). 


Liquid 


°G^i  -  °djjjp  »  4920  -  2.9T 
°GC  "  °C  =  27300  "  6*5T 


o~L  1  orgas  _ 

“n  •  7  ^  ■  ' 


2077  +  15.98T 


Subregular  solution,  Gm  given  by  Eq.  Al. 

°4 -  = 0 


LCTi  *  ‘23185  "  10 -ST,  iLcri  =  30516  -10.5T 
)lNri  =  “68482  "  7-SI =  40975  +  7.5T 


(Ref.  8) 


°4tC  -  °($?  =800 

0(^iC  "  °^TiP  "  °cf  c  "43972  +  2.528T  ' 

°4S  "  °Gti>  '  7  0(^aS  “  -79666  +  21.4ZT 

Sublattice  model,  On  given  by  Eq.  A8  a=c=l 

OjTifcc  l.Tifcc  n 
LCN  ■  lcn  =  0 

0lCVaCC  =  ~16099  +  0.36T,  1L^CC  =  -38801  +  11.03T 

^  ■  -25oo°  *  4-4T-  xCacc  *  ° 


(Ref.  8) 


5574 


b 


» 


» 


I 


t 


» 


i 


L 


hS2 


Ref  state  for  Ti 


°nhcp  -  2  °r,hcp  °rGr 
l,Ti2C  2  GTi  "  GC 


°r,?'P  -  2  °r,^P  .  1  Ofigas  .  .93490  ,  24.71T 

l  2 

sublattice  model,  Gm  given  by  Eq.  A8  a  =  2 ,  c  =  1 

o.Tihcp  l.Tihcp 
LCN  »  LCN  =0 


o. T ihcp 
LCVa 

o, T ihcp 
LNVa  = 


bcc 


'  -50000,  1Ljj^cp  =  0 
-9985  -3.44T,  1Ljj^cp 


=  0 


°r’XiC  -  °r.^P  *  1040  -  0.9T 

GTiC,  “  °G^?P  -  3  °G^r  =  -44996  ♦  0.49T 

O 

°GTiN,  ”  °GTiP  “  3/2  °fiN3S  =  -104155  ♦  31.22T 
3  *’2 

sublattice  model,  Gm  given  by  Eq.  A8  a  «  1,  c 

o.Tibcc  _  l.Tibcc  n 
CN  "  lCN  =  0 


oLTibcc 


CVa  »  -50000,  1Lj^CC  =  0 

o.Tibcc  l.Tibcc  n 
^NVa  ~  *  u 


JNVa 


Tetragonal 

°^t*N  *  2  °G*‘P  -  l  °G*a: 5  ■  -120692  ♦  43.06T 


'2 


Stoichiometric  Phase 


Gas  °G||S  -  °C^?p  a  112685  -  48.024T  +  1.808  TlnT 
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(Ref.  6) 
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CALCULATION  OF  THE  AljO-j-ZrOj-Y^  SYSTEM 
1 . 0  Introduction 

Zirconium  Dioxide  is  an  important  technological 
ceramic  because  of  its  refractory  nature  and  abrasive  prop¬ 
erties.  In  combination  with  additions  of  A^O^,  SiC,  SiO  , 

Y2°3'  Ca0  arK*  ;*'ts  ProPert ^es  can  k<2  tailored  for  a  wide 
variety  of  uses  as  a  thermal  and/or  an  oxidation  barrier  or 
a  structural  ceramic .Recent  developments  (1-3)  have  shown 
that  by  combination  with  Al2C>3,  Si^N^  and  SiC  enhanced  fracture 
toughness  can  be  developed  in  Zr02  composites.  Heuer  (4)  has 
shown  that  tough  zirconia-containing  ceramics  can  be  fabricated 
which  derive  their  usefull  strength  from  the  stress-induced 
martensitic  transformations  of  zirconia  percipitates  from  the 
tetragonal  (T)  to  monoclinic  (M)  structures.  Electron  diff¬ 
raction  and  transmission  electron  micrographic  studies  have  shown 
that  in  partially  stabilized  ZrC>2  with  tetragonal  (T)  phase 
particles  in  a  cubic  (A)  matrix  and/or  tetragonal  (T) 
structures  in  a  fine-grained  alumina  matrix  microcrack  formation 
around  crack  tips  leads  to  enhanced  toughness.  In  view  of  the 
wide  range  of  applications  for  Zr02  and  ZrC>2  composites  it 
seems  worthwhile  to  add  this  compound  to  the  exisiting  data 
base  (5,6)  described  in  Sections  II  and  ITI.  Accordingly  the 
follov/ing  analysis  of  the  Zr02~Al20^  and  Zr02~Y20^  were  carried 
out  and  coupled  with  the  available  description  of  the  y2°3~Ai'203 
shown  on  page  11  to  calculate  isothermal  sections  in  the  A12C>3- 
Zr02-Y203  system. 

\ 

2.0  Analysis  of  the  Y2°3~A^2°3'  ZrQ2-Y203 and  Zr02~Al203  System 

The  analysis  of  the  binary  systems  was  performed 
along  the  line  described  in  Section  II  using  available  phase 
diagrams  (7-11)  and  thermochemical  data  (12)  for  the  compounds 
of  interest,  in  the  course  of  this  study  it  was  noted  that  Y2C>3 
exhibits  a  high  temperature  transition  from  the  Mn203  structure 
to  a  high  temperature  form  on  heating  above  2500K  (8) . 
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Figure  6.  Calculated  l/SY^O^  -  l/5Al^0^  Phase  Diagram 
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Figure  1.  Recalculated  (1/5)  Y202- (1/5) AljO^  p^ase  Diagram 

This  fact  has  been  reflected  by  recalculating  the  YO-AO  diagram 
in  Figure  1  above  by  applying  the  new  description  for  YO  in  place 
of  that  given  in  Section  II  to  reflect  allotropism  in  YO.  Although 
this  is  a  small  change  it  will  be  employed  in  future  consideration.-. 

Tables  1  and  2  provide  the  lattice  stability,  solution  phase 
and  compound  phase  parameters  required  to  calculate  the  ZO-YO 
and  ZO-AO  phase  diagrams  which  agree  with  the  experimental  results 
available  (8,9) .  The  recalulated  YO-AO  phase  diagram  is  obtained 
by  employing  the  revised  lattice  stability  description  for  YO  shown 
on  page  132  with  the  previoulsy  published  (5)  solution  and  comp¬ 
ound  phase  parameters  shown  on  pages  6  and  7. 

3.0  Calculation  of  the  Al20.j-Zr02-Y202  System 

Figure  4  and  5  display  isothermal  sections  in  the 
A^2°3""Zr02~Y2°3  system  at  temperatures  between  2700K  and  1700K. 
These  sections  illustrate  the  ranges  of  liquid  and  solid  solution 
stability  and  the  extent  of  compound  intrusion  as  a  fuction  of 
temperature.  Expansion  of  the  data  base  by  analysis  of  other 
ZrC>2  binary  systems  sch  as  ZrC^-C^O^,  Zr02~CaO,  Zr02~MgO, 

Zr02-Si02,  Zr02-SiN4  and  Zr02-Be0  would  permit  synthesis  of  many 
technologically  important  quasi-ternary  systems. 
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TABLE  1 


Summary  of  Lattice  Stability  Parameter s}T  in  Kelvins 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms)  ,  Tin  Kelvins) 


A0= ( 1/5 )  A1203,  Z0= ( 1/3 )  Zr02,  Y0=(l/5)  Y^,  L=  Liquid 
C=  Corundum,  Y=high  temperature  (hexagonal)  ^2°3'  B“bcc 
(Mn203)  Y2°3'  A=Zr02  cu^ic#  T=Zr02  Tetragonal,  M-Zr02  Monocl 
20Z0LA- (1/3)  Zr02  (liquid) - (1/3)  Zr02  (CaF2)  Cubic 


AOAOLC=23640-10. 209T 
AOAOLA=  -  9.832T 
AOAOLT=  -10.586T 

ZOZOLA=29008-  9.832T 
ZOZOLC=  -10.209T 
ZOZOLY=17991-  8.368T 
ZOZOLB=24853“  9.205T 


ZOZOAT=1987-0 .  753T 
ZOZOTM=2008-1 . 381T 

YOYOLY=22694-  8.368T 
YOYOLB=26878-10. 042T 
YOYOLA=22615-  9.832T 
YOYOAT=-8368-  0.753T 
YOYOTM=-8368-  1.381T 


TABLE  2 

Summary  of  Solution  and  Compound  Phase  Parameters 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms) ,  T  in  Kelvins) 

LZOAO=l7573 ,  LAOZO=39748,  TZOAO=TAOZO=AZOAO=AAOZO=62760 
CAOZO=CZOAO=62760 
Solution  Phases 


L  ZOYO"  14016  +  4.184T 
LYOZO=  14016  +  4.184T 
BZOYO=  2929  +  8.368T 
BYOZO=  2929  +  8.368T 
AZOYO=-12552  +11.297T 
l.YOZO-  837  +  4.602T 


YZOYO=  2929  + 
YY0Z0=  2929  + 
TZOYO=-12552  + 
TYOZO=  837  + 
MZOYO=  4184  + 
MYOZO=  837  + 


8.368T 

8.368T 

11.297T 

4.602T 

11.297T 

4.602T 


Compound  Phase  =  P  =  (1/19)  (3Zr02*2Y2C»3)  =  ( 1/ 1 9 )  ZOgYD1Q- 

Z0o  . _ „ Y0_  r„,  Base  =A ,  Compound  Parameter  C-  31798  + 
0.474  0.526 

11.09T,  Gibbs  Energy  of  Formation  from  (1/3)  Zr02 (M)  and 
(1/5)  Y203(B),  AGf =-286-2 . 86T 


Y0=1/5Y„0  A+B  Z0*l/3  ZrO 


L / 3  ZrO 
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CALCULATION  OF  THE  HEAT  OF  FUSION  OF  EUTECTIC  FLUORIDE 


GLASSES  IN  THE  ZrF4~BaF2  AND  THE  ZrF^-LaF^  SYSTEMS 
1.0  Introduction 

The  current  analysis  of  the  ZrF^ -BaF 2  and  ZrF^-LaF^ 
systems  shown  on  page  65  is  the  first  step  in  calculation  of 
the  ternary  LaF2“ZrF4~BaF2  system  shown  on  page  69.  Very 
little  experimental  data  exist  for  these  systems.  One  method 
of  improving  the  level  of  confidence  in  the  calculation;,  and 
the  current  degree  of  knowledge  concerning  these  systems  won  id 
be  to  measure  the  heat  liberated  by  solidification  (oi  absorber 
during  fusion)  of  liquids  of  the  eutectic  composition  at  the 
eutectic  temperature.  Such  measurements  can  readily  be 
performed  by  differential  scanning  claorimetry.  The  calcula¬ 
tions  presented  below  provides  values  for  the  heat  of  fusion 
based  on  the  description  contained  in  Section  V. 

2 . 0  Calculation  of  Heat  of  Fusion 

The  heat  of  fusion,  AH,  at  the  eutectic  temperature, 
Tg,  and  the  eutectic  composition,  XE,  for  the  ZF-LF  case  is 
defined  as 

(1) 

(2) 

where  XE=  0.25,  TE=850K  and  LZFLF=-11088  J/g.at.  LLFZF= 

22595  defining  the  enthalpy  of  the  liquid.  The  last  two 
terms  are  (0.75)  (12845)  and  (0.25)  (12560),  from  page  61, 
respectively  so  that 

AH=  -500  +  9634  +  3140=  12274  J/g.at. 

since  ZFq  75Lfq  25  corresP°nds  to  (3/20)  ZrF4~(l/16)  LaF^, 
then  one  gram  atom  corresponds  to  25.1  grams  of  ZrF^  and 


AH=  HL[TE,  xe]  -  (l-xE)  H^p  -  xE  Hlf 


AH=  xE(l-xE)  t(l“XE)  LZFLF  +  X  LLFZF ] 

+  (1-xe)ahzf  +  xE  ahlf 


12.3  grams  of  LaF^  or  a  total  of  37.4  grams.  Hence  a  mixture 
containing  32.9  w/o  of  LaF^  would  exhibit  a  heat  of  fusion  of 
328  J/gm.^ on  the  basis  of  these  calculations.  Similarly  a 
calculation  of  the  heat  of  fusion  for  the  ZF-BF  case  suggests 
that  a  mixture  containing  42.7  w/o  BaF2  would  exhibit  a  heat 
of  fusion  of  281  J/gm.  near  800K. 

Measurement  of  these  heats  of  fusion  by  scanning 
calorimetry  would  be  very  valuable  in  improving  the  present 
calculations. 


CALPHAD  CALCULATION  OF  THE  EFFECT  OF  A1F3  ADDITIONS  ON 

THE  GLASS  COMPOSITIONS  OF  TERNARY  ZrF4-LaF3~BaF2  FLUORIDES 

(presented  at  the  Second  Interantional  Symposium  on  Halide 
Glasses.  Rensselaer  Polytechnic  Institute,  Troy  New  York 
12181,  USA,  2-5  August  1983,  Sponsored  by  AFOSR,  ONR, 

Corning,  GTE,  and  Spectran  Corporation.) 

1 . 0  Introduction 

The  discovery  of  a  family  of  non-oxide  glasses 
based  on  mixtures  of  ZrF^  or  HfF^  with  other  metallic  fluorides 
by  Poulain  and  coworkers  offers  great  potential  in  optical 
fiber,  windows  and  source/detector  applications.  Due  to 
limited  phase  diagram  data  available  for  binary,  ternary  and 
multicomponent  fluoride  systems  employed  to  synthesize  those 
glasses  progress  in  identifying  new  compositions  has  proceeded 
along  empirical  lines.  In  order  to  remedy  this  situation,  the 
CALPHAD  method  (1-3)  for  coupling  phase  diagram  and  thermo¬ 
chemical  data  has  been  applied  to  develop  a  data  base  covering 
metallic  fluorides.  The  objective  is  to  permit  computation  of 
multicomponent  phase  diagrams  which  can  be  used  to  identify  the 
composition  range  where  the  liquid  is  most  stable.  The  latter 
offers  opportunities  for  prediction  of  new  glasses.  Currently, 
the  data  base  covers  combinations  of  0.2  ZrF4(ZF),  0.25  LaF3 ( LF ) , 
0.333  BaF2(BF),  0.333  PbF2(PF),  0.5  NaF  (NF) ,  0.5  RbF  (RF) ,  0.5 
CsF  (CF)  and  0.5  KF  developed  along  the  lines  described  earlier 
for  III-VI,  II-VI  and  SIALON  systems  (2-3).  This  has  been  pre¬ 
sented  in  Section  V  Figure  11  on  page  69  shows  calcualtion  of 
isothermal  sections  in  the  LF-ZF-BF  system  illustrating  the 
range  of  composition  where  the  liquid  has  the  greatest  stability. 
The  compositions  agree  with  those  where  Poulain  and  cqworkers 
found  glass  formation  (4,5)  and  were  presented  in  1982  (6)  and 
recently  published  (7).  The  current  data  base  has  been  extended 
to  cover  A1F3  (AF=0. 25A1F3)  by  employing  experimental  data  pro¬ 
vided  by  Thoma  (8)  to  evaluate  AFAFLT=13389-8 . 368T  J/g.at.  and 
LAFZF=25104  J/g.at.=LZFAF,  LAFBF=4184  J/g.  at . =LBFAF  and  LAFLF= 
20920  J/g.  at.=LLFAF.  These  parameters  permit  calcuation  of  the 
range  of  maxium  stability  of  the  liquid  phase  for  ternary  LF-ZF- 
BF  glasses  with  additions  of  4 ,  9  and  15  atom  percent  of  AF . 


2 . 0  Calculation  of  the  Range  of  Liquid  Stability 


In  line  with  the  development  in  Section  V  the 
equilibrium  between  a  liquid  with  atom  fractions  of  LF, 

X2  atom  function  of  ZF  and  atom  fractions  of  BF  and  x4 
atom  fractions  of  AF;  and  solid  LF  having  the  J  structure 
solid  ZF  having  the  M  structure  and  solid  BF  having  the  K 
structure  is  given  by  equations  (1)  to  (3) ,  when  LLFBF=LBFLF ; 
LZFAF=LAFZF;  LLFAF=LAFLF  and  LBFAF=LAFBF  as  in  the  case  at 
hand ! 

LFLFLJ=-RT  In  x1  -  x^2  (x1  +x2 )  ( l-x2+x2  (x1  +x2 )  *)  LLFZF 

-X2  (x2(x1+x2)“1  -  x^  (x1-t-x2)”1LZFLF 

-x3(l-x1)LLFBF  -x4  (1-x.,)  LLFAF 
+x2x3 (x2+x3) -1  (x2LZFBF  +  x3LBFZF) 

+X2X4  LZFAF  +X3X.LAFBF  (l) 

ZFZFLM=-RT  In  X2~X1X2 (X1 +X2) ~1 ( 1-X2+X1 {x1+x2) _1)  LZFLF 

-x1(x^(x^x2)  1“X1X2)  (X.j+X2)  1LLFZF  -x4  (l-x2)  LZFAF 

-x2x3(x2+x3)  1  ( l-x2+x3  (x2+x3) -i)  LZFBF  +  X-jX^LLFBF 

+x3 (x2+x3) “1 (x2~x3 (x2+x3) “1)lbfzf+x1x4llfaf+x3x4lafbf 

(2) 

BFBFLK=-RT  In  Xj-X^l-X-j)  LLFBF  -x4  ( 1-X.j )  LAFBF 

~X3X2^X3+X2^  ^ ( 1”X3+X2 (x2+x3) -^)  LBFZF  +x^x4LLFAF 
+X2(X2+X3)"L(X3  -  x2(x2+x3)“1)  LZFBF  +  x2x4LZFAF 
+X1x2(x1+x2)”1(x1LLFZF  +  x2  LZFAF) 
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Equations  (1-3)  are  solved  at  773K,  873K  and  973K  for  x 
0.00,  0.04,  0.09  and  0.15  respectively  in  Figures  1-G  which 
show  the  tie  lines  from  the  liquid  phase  composition  to  soiid 
J  (LF) ,  solid  M  (ZF)  and  solid  K  (BF) .  Figures  7  and  8 
show  only  the  composite  liquid  phase  boundries  between  773K 
and  973K  for  AF=0.00  to  0.15  and  Figure  9  shows  the  liquid 
phase  boundries  at  773K  compared  with  limited  AF=0.04.  The 
present  calculations  for  AF=Q.Q9  and  AF=0.15  suggest  that 
at  these  levels  glass  formation  will  be  easier  and  that:  the 
center  of  the  glass  forming  region  will  shift  toward  lower 
LF  concentrations 

Note  in  Proof 

At  the  Second  International  Conference  on  Halide  Glasses, 
where  this  presentation  was  made,  it  was  noted  that  the  compound# 
ZrF^*BaF2  and  ZrF^LaF^  had  been  reported  to  be  stable  phases. 
Exact  references  were  not  provided  but  one  participant  suggests 
a  congruent  melting  point  for  ZrF4*BaF2  (ZFQ  625BF0  375^  aL 
817K.  Comparison  with  Figure  5  on  page  65  shows  that  if  such 
a  compound  were  inserted  and  a  congruent  melting  point  would 
be  required,  the  liquidus  temperature,  currently  about  1000K 
would  have  to  be  lowered.  Such  an  alteration  would  require 
small  changes  in  the  values  of  LZFBF  and  LBFZF  in  Table  2  on 
page  61.  However,  the  eutectic  temperature  of  800K  predicted 
in  Figure  5  page  65  and  the  succeeding  results  of  Section  V  and 
the  present  section  would  not  be  materially  affected.  A 
reference  to  differential  thermal  analysis  studies  of  ZrF^-BaF7 
by  Y.  Kowmoto  and  F.  Sakaguchi  submitted  for  publication  to 
the  Bulletin  of  the  Chemical  Society  of  Japan  was  reported 
but  has  not  been  aquired  to  date. 
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Figure  3  Calculated  Solid-Liquid  Equilibria  in  LF-ZF-BF 
with  AF=0 . 04 


Figure  5  Calculated  Solid-Liquid  Equilibria  in  LF-ZF-BF 
with  AF=0.09  (above)  AF=0.15  (below) 
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